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Variationsprinzip für Kernroationen 
Wilhelm Macke und Richard Lenk 
Inhaltsübersicht 


Für deformierte axialsymmetrische Kerne werden Lösungen mit einem | . 
Zusatzdrehimpuls senkrecht zur Figurenachse aus der Variationsforderung 
abgeleitet, daß die Energie einen Extremwert unter den beiden Bedingungen 
konstanter Normierung und vorgegebenen Drehimpulses annehmen soll. Die == 
exakten Lösungen dieses Variationsverfahrens sind Eigenlösungen des Hamil- 
tonoperators H. Geht man aber von einem nichtkugelsymmetrischen Modell- J 
operator H° aus, so entstehen Lösungen, die einem mit der Winkelgeschwindig- 3 
keit @ rotierenden Kern entsprechen. Die Erwartungswerte von Energie und 4 
Drehimpuls hängen in erster Näherung quadratisch bzw. linear von der Winkel- _ 
geschwindigkeit ab. Der Proportionalitätsfaktor wird als Trägheitsmoment 
des Kerns interpretiert. Der Trägheitstensor ist entartet, der Kern rotiert nur 
um Achsen senkrecht zur Figurenachse. Es wird eine allgemeine Formel für 
das Trägheitsmoment angegeben, die den Einfluß der Wechselwirkungen be- 


Die Analyse der experimentellen Daten führt für die Atomkerne mit den | 
Massenzahlen 150 < A < 190 und A > 222 zu folgenden Ergebnissen!): Die 
Kerne besitzen große innere Quadrupolmomente, sind also stark deformiert. 
Dabei bleibt die Kernform jedoch im wesentlichen rotationssymmetrisch um — 
eine ausgezeichnete Achse. Die niedrig angeregten Zustände dieser Kerne lassen _ 
sich zu Folgen 
+ 1) 
zusammenfassen. Dabei ist J der Kernspin, g das Kernträgheitsmoment, Z, 
eine von J unabhängige Konstante. (1) beschreibt mit einer Genauigkeit von _ 
wenigen Prozent die gemessenen Energieniveaus einschließlich des untersten = 
Niveaus einer solchen Rotationsbande. Am wichtigsten ist der Fall, daß der 
unterste Zustand des Rotationsspektrums (1) der Grundzustand des Kerns 
ist. Das Trägheitsmoment kann nach (1) aus den Experimenten entnommen _ 
werden, für jede Rotationsbande hat g einen bestimmten Wert. Verglichen 


1 
mit dem Trägheitsmoment g;tarr für starre Rotation gilt grob g = 5 Istarr 


Zur theoretischen Erklärung der gemessenen Trägheitsmomente existieren 
im wesentlichen folgende Versuche: a) Die hydrodynamische Theorie beschreibt 


N Alder, Bohr, Huus, Mottelson u. Winther, Rev. mod. Physics 28, 432 Chap. V | 
(1956); S. A. Moszkowski, Artikel in Handbuch der Physik (S. Flügge) Springer 
(1957) Berlin/Göttingen/Heidelberg. 


l Ann, Physik, 7. Folge, Bd. je 


4 
7 
| 
| 


2 > Annalen der Physik. 7. Folge. Band 9. 1961 TK A 


die Kernrotation als wirbelfreie Strömung einer Kernflüssigkeit, diese Methode 
liefert um einen Faktor 3—5 zu kleine Trägheitsmomente. b) Das cranking 
model?) postuliert eine Rotation für das deformierte mittlere Potential, in dem 
sich die Nukleonen unabhängig voneinander bewegen. Konkrete Rechnungen?) 
mit impulsunabhängigen Potentialen (harmonischer Oszillator und Potential- 
kasten) haben g = starr ergeben, wenn der ruhende Kern durch den Grund- 
zustand im deformierten Potential beschrieben wird. Dieses Ergebnis konnte 
von Lüders®) allgemein für große Teilchenzahlen begründet werden. Eine 
mögliche Impulsabhängigkeit des Potentials?) ändert nur die Form der Theorie, 
nicht das Ergebnis. c) Schließlich ist das cranking model ergänzt worden durch 
die Annahme, daß die Kernmaterie superfluide Eigenschaften ®) besitzt. Durch 
die Einführung einer entsprechenden Energielücke an der Fermioberfläche 
läßt sich Übereinstimmung mit den experimentellen Werten für g erreichen. 
Von dieser Möglichkeit soll jedoch im folgenden abgesehen werden. 

In Kapitel 2 werden die Eigenschaften eines Kernmodells bei deformierten 
Kernen zusammengestellt. In Kapitel 3 werden neue Lösungen aus einem 
Variationsprinzip abgeleitet, das hinsichtlich der exakten Lösungen dem Va- 
riationsprinzip zur Schrödingergleichung äquivalent ist. Dabei wird das 
Verfahren zunächst in der einfachsten Näherung durchgeführt, wobei der 
Hamiltonoperator einfach durch den Modelloperator ersetzt wird. Die Deu- 
tung des Verfahrens führt in Kapitel 4 zu der Vorstellung eines mit bestimmter 
Winkelgeschwindigkeit rotierenden Kerns. In Kapitel 5 wird der Einfluß der 
Wechselwirkungen zwischen den Nukleonen auf das Trägheitsmoment allge- 
mein untersucht, die in Kapitel 3 verwendete Näherung wird also fallenge- 
lassen. Schließlich wird in Kapitel 6 gezeigt, daß ein analoges Verfahren für 
die Kerntranslation zum richtigen Zusammenhang zwischen Energie und Ge- 
samtimpuls führt. 


2. Kernmodelle und deformierte Kerne 


Jot 
Dem exakten Problem Sy 

wird ein Modellproblem te 


H® — H’=-T+W (2) 
zugeordnet, das möglichst einfache Eigenschaften besitzen soll. Entsprechend 
H=H’+U! uUu=-V—-W (3) 


kann der Hamiltonoperator H in Modelloperator H® und Störung U® zer- 
legt werden. Beriicksichtigt man U® durch Stérungsrechnung’), so gilt 


E= E + (0, B®) mit B= 
E—H® 
2) D. R. Inglis, Physic. Rev. 96, 1059 (1954); 108, 1786 (1956). 
*3) A. Bohr u. B. R. Mottelson, Kgl. Danske Videnskab. Selskab., Mat.-fys. Medd. 
30, No.1 (1955); R. D. Amado u. K. A. Brueckner, Physic. Rev. 115, 778 (1959). 
4) G. Liiders, Z. Naturforsch. 15a, 371 (1960). 
5) D. R. Inglis, Nuclear Phys. 8, 125 (1958). 

6) V. I. Belyaev, Kgl. Danske Videnskab. Selskab., Mat.-fys. Medd. 31, Nr. 11 (1959). 
A. B. Migdal, Nuclear Phys. 18, 655 (1959). 
7) W. Macke, Ann. Physik, Leipzig (6) 20, 80 wt 
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Das Problem besteht darin, die Aufspaltung (3) — innerhalb der durch die 
Modellstruktur von H® eingeengten Möglichkeiten — so einzurichten, daß die 
Differenz E — H® in (4) möglichst klein wird. Dies läßt sich optimal jeweils 
nur für einen einzigen Zustand erreichen, etwa für den Grundzustand des 
Systems. Im Einteilchenmodell zum Beispiel sind die ®° Determinanten aus 
Einteilchenwellenfunktionen, durch diese Forderung ist die Struktur von H 
weitgehend festgelegt. Je nach der Näherung, die bei der Auswertung von (4) 
benutzt wird bzw. erforderlich ist, kommen in diesem speziellen Fall zur ge- 
nauen Bestimmung von H® zwei Methoden in Betracht: Hartree-Fock- 
Näherung und Brueckner-Theorie. In jedem Falle gilt | 
E = (®°, AIO NEA (5) 
Der exakte Hamiltonoperator H ist invariant gegen beliebige Trans- 
lationen und Drehungen, also gilt 


[H, $] =0 [H, 3] = 0. (6) 
Die Eigenlösungen ® der Energie aus (1) können daher gleichzeitig Eigen- 
funktionen des Gesamtimpulses ® sowie von Quadrat und z-Komponente des 
sesamtdrehimpulses % sein. Beim Übergang zum Modelloperator H® geht die 
Symmetrie von H zumindest teilweise verloren. Durch Einführung eines 


raumfesten mittleren Einteilchenpotentials, das in H® enthalten ist, wird auf 
jeden Fall . 
[H°, B] = 0. (7) 
Die Modellfunktionen sind nicht Eigenfunktionen von ®; außerdem ver- 
schwindet der Erwartungswert von ®. Das Verfahren liefert nur Zustände, bei 
denen der Erwartungswert des Schwerpunktes einen festen Wert hat. Dabei 
führt der Schwerpunkt möglicherweise Schwingungen aus, die im exakten 
Problem nicht auftreten und beim Modellproblem zur Erhöhung der Energie 
und zum Auftreten unrealistischer Zustände®) (spurius states) führen. Die 
freie Translation des Schwerpunkts und die zugehörige Schwerpunktsenergie 
ist dagegen im Modell nicht enthalten. 

Bei deformierten Kernen verliert H® die Kugelsymmetrie, die H besitzt. 
Bleiben — wie im folgenden angenommen wird — die Kerne axialsymmetrisch, 
dann enthält H° eine ausgezeichnete Richtung e, die Rotations- oder Figuren- 
achse des Kerns. Zerlegt man den Drehimpuls nach 


=eoce J—e X (e X FJ) =eoeF¥+Q3Y (8) 
in eine Komponente in Richtung von e und eine senkrecht zu e, wobei die 


durch „‚o““ getrennten Vektoren nicht miteinander zu verknüpfen sind und der 
Projektionstensor Q in (8) durch 
Q=I-—eoe I = Einheitstensor (9) 
if 
(2°, e 3] =0 (10) 


Die Modellfunktionen sind also nach 


definiert ist, so gelten die Vertauschungen 


8) J. P. Elliot u. T. H. R. Skyrme, Proc. Roy. Soc. A 232, 561 (1955); S. Garten- 
haus u. C. Schwartz, Physic. Rev. 108, 482 (1957). 
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Eigenfunktionen von e §. Allein aus den Vertauschungsrelationen der Dreh- 
impulskomponenten folgt dann bekanntlich, vgl. z. B. ?), daß der Erwartungs- 
wert von Q% in allen Zuständen ®Y, verschwindet: 


(Pix, Dix) = 0. (12) 
Es existieren zu H® also nur Zustände ohne Drehimpuls senkrecht zu e. (Die 
Gleichungen (11) und (12) gelten an sich stets, also auch bei kugelsymmetri- 
schem H®, nur kann dort e beliebig gewählt werden, während es hier durch 
H® selbst festgelegt wird.) Das exakte Problem enthält dagegen auch Zustände, 
bei denen ein Drehimpuls senkrecht zur Symmetrieachse des Kerns auftritt 
(Rotationszustände). Ähnlich wie bei der Translation des Schwerpunkts wer- 
den diese Kernrotationen vom Modell H® nicht geliefert. Dagegen kann auch 
hier nicht verhindert werden, daß der Kern kleine Pendelbewegungen um die 
durch H® festgelegte mittlere Achsrichtung ausführt, 4° kann also auch un- 
physikalische Zustände dieser Art mit enthalten. 


3. Variationsprinzip 


Die Schrédingergleichung (1.1) ist Variationsprinzip 

6|H — AN] = 6[(@, H ©) = 0 (1) 
äquivalent. Die Näherungslösungen ®® dieses Problems führen bei deformierten 
Kernen zu (2.11) und (2.12). Um Lösungen aufzufinden, die nicht nur einen 
Drehimpuls in Richtung von € besitzen, soll bei der Variation (1) der Erwar- 
tungswert des Drehimpulses auf einen Wert A gemäß 


(8,309) (2) 
festgelegt werden. Berücksichtigt man die Nebenbedingung (2) bei der Varia- 
tionsaufgabe (1) mit einem Lagrangeschen Parameter @, so entsteht 


—a3—AN|=0. (3) 


Die Variationsforderung (3) ist gelegentlich bereits in der Literatur!®) benutzt 
worden. Bei vollständiger Variation sind die beiden Variationsprinzipien (1) 
und (3) praktisch gleichwertig. Die Lösungen von (1) sind stets Eigenfunktionen 
vs zu einer Komponente von J, etwa zu J,, wobei allerdings die z-Richtung völlig 
beliebig ist. (Auch die Überlagerung mehrerer Lösungen mit verschiedenen 
Eigenwerten von J, gibt nur eine Eigenfunktion zu einer neuen z’-Richtung.) 


Dagegen ist bei (3) die z-Richtung durch die Richtung von @ Fuge va eh (3) be- 
sitzt nur Lésungen mit 
Im übrigen sind die ® die Lösungen (2.1) des Problems (1). Als Erwartungs- 
wert des Drehimpulses werden nach nur die Werte 


») W. Macke, Quanten, 196 und 476, Akad. Verl.-Ges. Geest & Portig K.-G., Leipzig 
(1959). 
10) E. P. Gross, Nuclear Phys. 14, 389 (1959); D. J. Thouless, Nuclear Phys. 21, 
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eh- geliefert, wobei die Richtung von X durch die von @ festgelegt ist, lol bleibt 
1gs- willkürlich wählbar. 

Bei teilweiser Variation, also beim Aufsuchen von Näherungslösungen, ist 
(12) die Nebenbedingung (2) dagegen entscheidend. Sie führt bei deformierten 
Kernen zu neuen Lösungen, die durch H® noch nicht erfaßt werden. (3) geht 


ie aus (1) hervor, indem man nach 

ırch H>H-—-o% (6) 

yet zum Hamiltonoperator den Term — @& hinzufügt. In der einfachsten 4 
Näherung, nämlich H ~ H°®, sind also Lösungen ®® zu { 
wer- 

[H°— 8° = (7) 

pa aufzusuchen. Falls ®% als kleine Störung behandelt werden kann, gilt 


»=[1 + @ (8) 
in erster Näherung der Schrödingerschen Störungsrechnung. Der Pro- 
nzip | jektionsoperator P schließt den ungestörten Ausgangszustand ° aus: 

(1) PP=0=-MP. (9) 
rten Bezeichnet man den Zustand ©° durch ), dann lauten (8) und (9) 


P,=0=<P. 
(2) I Die gestörte Funktion ist entsprechend | 
aria (8°, ® )=(1 + 3/= = 1 
(3) Erwartungswerte physikalischer Größen werden nach tn Ad 
(9°, 4 <A> + @ IE! AD + (A 3 + A&G J) (12) 
n (1) ©”) 1+ 


gebildet. Dabei ist das quadratische Glied im Zahler nur dann sinnvoll, wenn 
die in (8) weggelassenen Anteile der Wellenfunktion keine Beiträge gleicher y 


Ordnung zu A liefern. Vernachlässigt man Terme in w* und höher völlig, so 
vereinfacht sich (12) zu 
A = <A) + (A £9' + <@Ye0' A> + (13) 


(4) Daraus ergibt sich der Erwartungswert des Drehimpulses bis zur Ordnung 


ozu 
J=<H + 60 DM. (14) 
Für den Erwartungswert der Energie ist A = H® zu setzen. Wegen 
0 
(5) { H% = = Pa 0 (15) 


verschwinden die in w linearen Terme sowie entsprechende quadratische, die 
eipzig | in ( 13) fehlen und bei denen A ebenfalls direkt an > anschließt. Also gilt i 


E= He — <H®, + 3506 (16) 
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mit 
Die beiden in (14) und (16) eingeführten Tensoren stimmen also überein. 
Der Tensor @ ist offenbar symmetrisch. Um weitere Eigenschaften zu 


finden, multipliziert man ihn zunächst mit dem Einheitsvektor e. Wegen 
(2.10) und (2.11) entsteht dann 


mit (10) gilt also 


e6=0=6e, ebe = 0. (19) 
Damit bleibt die Drehimpulskomponente in Richtung der Figurenachse nach 
(20) 
ungeändert. Zerlegt man ferner @ so wie $ in (2.8) in zwei Komponenten 
@=ecea+Qa, 

so ist nach 


die Komponente von @ in Richtung von e ohne Einfluß auf Drehimpuls und 
Energie. Aus Symmetriegründen muß ferner 


ada = g = konst. ae=0 a? = 1 (23) 


für alle auf e senkrechten Einheitsvektoren a gelten, speziell also auch für 
zwei orthogonale Vektoren a’ und a’’. Aus a’ und a’”’ können aber durch Linear- 
kombination neue Vektoren konstruiert werden, für die (23) ebenfalls gelten 
muß, etwa 


1 [73 
a= (a’ +0") 
Nun ist (23) eine Forderung 
HF: 
== Oa’ +. a” Oa” + a’ Oa” + 0” =g, (25) 
die sofort auf 
seat 


a’ da” = a” da’ =0 ea’ =ea’ a’=0 (26) 


führt. Wählt man also e als eine der drei Achsen, so ist 6 stets auf Haupt- 
achsen transformiert. Wegen (19) und (23) ist der Eigenwert g einfach durch 


darstellbar. Der gesamte Tensor 6 hat die Form are ni 


wobei Q der in (2.8) und (21) verwendete Projektionstensor ist. Im Unter- 
schied zu (14) sind in (28) die beiden Vektoren § skalar verknüpft, entsprechend 
der geforderten Spurbildung. 
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4. Deutung durch Kernrotationen | 

Nach (3.6) ist bei Variation der Energie mit vorgegebenem Drehimpuls — 

H durch H — @ § zu ersetzen. Das entspricht dem Ubergang zu einem mit der 
Winkelgeschwindigkeit & rotierenden Bezugssystem, in dem der Hamilton- 


operator gerade die Form H —@ § besitzt. Diese Behauptung soll zunächst 


unter Weglassung des Spins, also 


pag werden. Dazu werden die kanonischen ee zu 
— @ & gebildet: 


(2) 
Zusammenfassung der beiden Gleichungen gibt j are 3 
(3) 


zur Kraft — 0H/ér; treten Coriolis- und Zentrifugalkraft hinzu. 
Eine quantenmechanische Überlegung führt ebenfalls zur gleichen Deutung 


des Zusatzterms — w §. Die Zeitableitungen beliebiger Erwartungswerte ent- 


halten nach 
az 
3, 4] + = 


der für drei typische Fälle untersucht werden soll. Für A = r, ergibt sich nn Dr 


einen Zusatzterm 


wobei die Vertauschbarkeit u 
[p,w x r) =0 (7) 


benutzt wurde. Entsprechend gilt fiir A = p, ib 


Für A = § entsteht schließlich 


3]- + {3 (63) - (9) 


ein Ausdruck, der sich auf Grund der Vertauschbarkeit von @ mit $ und der 
Unvertauschbarkeit der einzelnen Komponenten von % als 
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schreiben läßt. Dabei sind in 


die Drehimpulsvertauschungen zusammengefaBt und in (10) verwendet. Die 
Gleichungen (9) bis (11) gelten alle auch fiir Bahndrehimpuls 2 und Spindreh- 
impuls S einzeln. In (6), (8) und (11) zeigen die Erwartungswerte eine zusätz- 
liche Zeitabhängigkeit, die der Beobachtung in einem mit @ rotierenden 
Bezugssystem entspricht. 

Das in Kapitel 3 durchgeführte Variationsverfahren kann also wie folgt 


gedeutet werden: Der Lagrangesche Parameter @ hat die Bedeutung der 
Winkelgeschwindigkeit desrotierenden Kerns. Dessen Struktur wird untersucht 
im mitbewegten Koordinatensystem, sie wird dort beschrieben durch die 
Wellenfunktion ®® von (3.8). Die Rücktransformation auf das ruhende Koor- 
dinatensystem liefert dort eine nichtstationäre Lösung 
x 1), (12) 
wobei rt symbolisch für die Koordinaten aller Teilchen steht. Für die Bildung 
von Erwartungswerten skalarer Größen (Energie) ist die Zeitabhängigkeit in 
(12) jedoch ohne Bedeutung, da das Argument der Wellenfunktion nur Inte- 
grationsvariable ist, deren Bezeichnung beliebig ist. Bei den Erwartungswerten 
von Vektoren führt die Transformation (12) zu einer REN, vom 
Typ 
hun 
falls der I = im rotierendem System zeitunabhängig ist. Das 
gegenüber (10) geänderte Vorzeichen entspricht der Rücktransformation auf 


die Lösungen ® von = 2 


Der Effekt (13) ist im Formalismus von Kapitel 3 nicht mit enthalten, er 
tritt lediglich bei der Betrachtung nichtstationärer Lösungen (12) auf. Er kann 
jedoch benutzt werden, um die entsprechend (3.22) wirkungslose Komponente 


e © in Symmetrierichtung durch die Forderung 


(15) 
festzulegen. Winkelgeschwindigkeit und Drehimpuls haben dann gleiche Rich- 
tung. Der Erwartungswert des Gesamtdrehimpulses mit den Lésungen (12) 
wird nach (13) zeitunabhängig. 

Der in (3.28) definierte Tensor wird nun zum Trägheitstensor des Kerns. 
Seine im Projektionstensor Q zusammengefaßte Entartung bedeutet, daß der 
Kern stets um eine zur Figurenachse e senkrechte Richtung rotiert. Das Träg- 
heitsellipsoid hat die Gestalt einer Stricknadel in Richtung von e. Der Zahlen- 
wert 


ist das Tragheitsmoment um eine zu e senkrechte Achse, J, die entsprechende 
Komponente von %. 
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Die Vorstellung eines mit bestimmter Winkelgeschwindigkeit rotierenden 
Kerns entspricht einer klassischen Beschreibung der Rotationsfreiheitsgrade 
des Kerns. Das durchgeführte Verfahren liefert Rotationszustände nur bei 
deformierten Kernen, d. h., die Mindersymmetrie der Kernform gegenüber der 
Symmetrie von H ist Voraussetzung für das Zustandekommen einer kollek- 
tiven Rotationsbewegung. Anders ausgedrückt: Nur bei Mindersymmetrie 
der Kernform erscheint ein Teil der Bewegung aller Nukleonen als kollektive 
Bewegung des Gesamtsystems. 


Die Auswertung des Variationsprinzips mala 


—®3—AaN]|=0 
ist bisher in Kapitel 3 nur in der einfachsten Näherung H — H° durchgeführt 
worden. Genauer wären Lösungen ® zum Hamiltonoperator 


aufzusuchen, um das Variationsproblem (1) zu lösen. Faßt man wie in (2) U® 


und — @ & zu U(w) zusammen und behandelt diesen Term nach der Störungs- 
rechnung, so entsteht die Lösung 


B(o) = U(w) 1 — B(o)]. 
Der Erwartungswert nat 


kann ausgewertet werden und ist entsprechend (3.14) und (3.16) als 


5. Einfluß der Wechselwirkungen 


zu interpretieren. 
Der in (5) auftretende @)> kann entsprechend = 
eB\. 
= BOY + (Z| 
= 
Anz = (Bo) + (BiG + 


nach Potenzen von @ entwickelt werden. <B,> enthält 3 überhaupt nicht, 
es stellt die Differenz zwischen der Modellenergie EP fiir den ruhenden Kern und 
dem exakten Eigenwert dar. Zur Auswertung der in (7) geforderten Ableitungen 
ist zunächst die Abhängigkeit der Energienenner von w durch eine Entwick- 


lung nach ungestörten Energienennern Se 
8) 
zu beseitigen. Das Ergebnis!) ist (Tae? sale 


B(w)) = <U B 
Bio) = +B, 
1) J. Goldstone, Proc. Roy. Soc. A 289, 267 (1957). ay a 
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wobei B(w) durch ln nit 

B(w) = (U (wo) — (U(w)>) Bw] (10) 

definiert ist und der Index L in (9) bedeutet, daß bei der Auswertung von (9) 
nur über alle Verbundcluster (linked cluster) zu summieren ist. 
Der Erwartungswert (Bo) unterscheidet sich nun von 


= 


dadurch, daß in der durch Entwicklung von (11) entstehenden Reihe jeweils 
an einer Stelle UP durch — § ersetzt wird. Das liefert zunächst <— 9), zu- 


sammen mit @ also das in (6) enthaltene Glied <— @ 9). «Bt ist eine 
Summe über Terme vom Typ 
— ID)... <UN)),. (12) 


Nun besitzt U® als Differenz von H und H® die beiden Operatoren gemeinsame 
Symmetrie 


[U°, e3] = (4, eg] — (4°. (13) 
es hat also von Null verschiedene Matrixelemente nur zwischen Zuständen 
mit gleichem Eigenwert AM vone §: 1 

3 


Daher treten in (12) nur Matrixelemente von % — <)> zwischen solchen Zu- 
ständen auf, die beide zum Eigenwert hM,des Grundzustands ) gehören. Für 
diese Elemente kann aber nach FR 


(Dy, — Ov.) = (Or, Q I = 0 (15) 


3 —< einfach durch Q% ersetzt werden, dessen Elemente zwischen 
Zuständen zum gleichen Eigenwert von e &, hier speziell hM,, ebenso wie die 
in (2.12) bereits betrachteten Diagonalelemente verschwinden. Wegen (14) und 
(15) sind also alle Terme (12) Null, <@Bo> ist also mit <— @ %) identisch. 
Auch wenn in (12) % — <Q) zweifach auftritt, kann es durch Q § ersetzt werden. 
Denn zu Nichtdiagonalelementen trägt ohnehin nur Q% bei, für Diagonal- 
elemente liegt aber wegen (14) stets der Fall (15) vor. 

Der in & quadratische Term von (7) liefert nach (6) das Trägheitsmoment 
als 


2 \aao ae o L 

_ wobei (9) verwendet wurde. Der Ausdruck (16) für 6 entspricht der doppelten 


Summe über alle Verbundcluster, in denen zweimal (UP—<UP®y) durch 
Q% ersetzt wird: 


= —2¢(1 Be) e0'(1 — Byeo') QF (1—e0' 
Die ersten Terme von (17) lauten oll .ctmcbtinen 
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Das Trägheitsmoment wird also durch eine Entwicklung nach Potenzen von 
U® dargestellt, die einfachste Näherung U° = 0 wurde in Kapitel 3 bereits 
untersucht. Der allgemeine Ausdruck (17) besitzt die gleichen Entartungs- 
eigenschaften, wie sie nach (3.28) der erste Term von (18) zeigt. e ® ver- 
schwindet, da nur Q% zu (18) beiträgt. 


Um die Deutung der durch (6) definierten Größe 6 als Trägheitstensor 
zu bestätigen, soll noch gezeigt werden, daß der Erwartungswert des Dreh- 
impulses 


3 = <}> + Oo (19) 


wieder zum gleichen Ausdruck (17) für das Trägheitsmoment führt. Aus (6) 
und (19) folgt dann ebenfalls rar 


ryt 

H=H-63+63=<H> +3060 (20) 
für die Abhängigkeit der Energie von w. Die Auswertung von § mit den Wellen- 
funktionen (3) liefert 


~ (@, “1+ (Blo) B(o)y 


BO) — Blo) IP + BO) BO) (91) 


Dabei enthält (21) noch die gestörten, von w abhängigen Energienenner. Der 
Term im Nenner wird in unterster Näherung nur mit (%) im Zähler kom- 
biniert und führt dazu, daß im letzten Term des Zählers § durch Q% zu 
ersetzen ist. In höheren Näherungen werden durch den Nenner von (21) gemein- 
sam mit den Zusatztermen, die beim Übergang zu ungestörten Energienennern 
entstehen, die unlinked cluster kompensiert, wie im einzelnen nicht gezeigt 
werden soll. Das Ergebnis 


= B(w)> + <B(w) €0' 3) — <B(w) QF Bw))|, (22) 


gestattet nunmehr eine einfache Entwicklung nach Potenzen von @ 


= (3 + @[— Bod + <By £0! Q Ye BS], + ©. e. (23) 


wobei die von w unabhängigen Zusatzterme zu (%) überhaupt verschwinden. 
FaBt man die beiden Terme in der Klammer von (23) zusammen, so entsteht 
genau der Ausdruck (17), die in (23) enthaltenen komplex konjugierten Glieder 
liefern einfach einen Faktor 2. 


Das in diesem Kapitel durchgeführte Verfahren, bei dem in einer einzigen 
Störungsrechnung entsprechend (2) der Übergang vom Modelloperator H® zum 


Operator H — &% durchgeführt wird, kann am einfachsten gedeutet werden. 
wenn der erwähnte Übergang nach 


| 


in zwei Einzelschritte zerlegt wird. Dabei entsteht zunächst ein rotierender 
Modellkern, wie er in den Kapiteln 3 und 4 ausführlich untersucht wurde. Der 
zweite Schritt erlaubt eine genauere Berücksichtigung der Kernstruktur, ins- 
besondere eine Veränderung der Kernstruktur durch die Kernrotation. 
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Die Auswertung des allgemeinen Ausdrucks (16), bei dem der Einfluß der 
Wechselwirkungen zwischen den Nukleonen berücksichtigt wird, stellt ein 
schwieriges Problem dar. Rechnungen mit einfachsten Näherungsvorausset- 
zungen finden sich bei Amado und Brueckner®), Rockmore!) und in 3), 
Rechnungen für endliche Kerne unter Berücksichtigung der Wechselwirkung 
liegen nicht vor. 
6. Schwerpunktsbewegung 


Beschreibt man die Schwerpunktsbewegung durch das Variationsprinzip 
6[H —» N] = 0, (1) 


wobei ganz analog zu (3.3) » zunächst einen Lagrange schen Parameter zur 
Berücksichtigung der Nebenbedingung 


(2) 
darstellt, so lassen sich Lösungen zu (1) allgemein durch den Ansatz 
gewinnen. Mit (3) transformieren sich die Impulse nach ak ape 
=e (p, + m, ® (4) 


Im Erwartungswert von H beeinflußt der Phasenfaktor von (3) nur den kine- 
tischen Anteil, insgesamt überführt (3) das Variationsprinzip (1) in 


6[(@, Ho) — 18,0) — 02] = 0 mit (6) 
Dabei ist aber M v?/2 eine belanglose Konstante, (5) ist also das ursprüngliche 
Variationsprinzip (3.1) ohne Nebenbedingung. Die Näherungslösungen von (5) 
beschreiben einen ruhenden Kern; allgemein sei ® eine genäherte oder auch 
exakte Lösung von (5) mit der Eigenschaft 


(8,20) = 0, (6) 


daß der Erwartungswert des Impulses verschwindet. Dann führen die Lösungen 
©’ von (3) unter Verwendung von (4) auf 


BE = (0, = [BP + Mv]d)= Mv (7) 


und 


Ho) = (© ‚[# +0 B+ 0) = @, 40) 45%, (8) 


wobei M die in (5) definierte Gesamtmasse ist. Das Variationsprinzip (1) führt 
also auf Lösungen, bei denen nach 


H = (®, H®) + ® ea dew (9) 


2) R.M.Rockmore, Physic. Rev. 116, 469 (1959). rs 
13) R. Lenk, Ann. Physik 9, 79 (1961). 
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der Zusammenhang zwischen Energie und Impuls die Masse M als Trägheits- 
parameter des Systems enthält. 


So wie nach Kapitel 4 H— © % den Hamiltonoperator in einem rotieren- 
den Bezugssystem darstellt, so ist hier H — v ®der Operator in einem mit der 
Geschwindigkeit vd bewegten System. Bildet man nämlich die zu (4.4) analogen 
zusätzlichen Zeitableitungen 


A=-z 
so erhält man für = ‘ 


(10) 


und für A = p, 


pat [Be, = 0. (12) 
Der Parameter v gewinnt damit die Bedeutung der Geschwindigkeit des be- 
wegten Kerns. Das Verfahren kann so gedeutet werden, daß die Struktur des 
Kernes in einem mitbewegten Koordinatensystem untersucht wird. 


Dresden, Institut für Theoretische Physik der Technischen Universität. 


Bei der Redaktion | am 17. Juni 1961. 
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Über ein Bi-Orthogonalsystem der Symmetriecharaktere 


Inhaltsübersicht 


Es werden Matrizen A,, und B,, eingeführt, welche die Linearbeziehungen 
zwischen den Normalformen eines Symmetriecharakters bzw. die Beziehungen 
zwischen den symmetrischen und antisymmetrischen Normalformen dar- 
stellen. Diese Matrizen lassen sich aus einer dritten Matrix, der Z-Matrix 
berechnen, welche nur von den Eigenschaften der verwendeten Tableaus 
abhängt, und sich mühelos konstruieren läßt. Mit Hilfe der Matrix A,, wird 
die Reduktion auf das linear unabhängige System durchgeführt. Wenn der 
reduzierte Teil der H-Matrix die Einheitsmatrix ist, bilden die symme- 
trischen und antisymmetrischen Funktionen ein Bi-Orthogonalsystem. Für 
die Symmetriecharaktere von Systemen bis zu 6 Nukleonen sind die Bi- 
Orthogonalsysteme angegeben. 


In der Theorie der Symmetriecharaktere wird im allgemeinen die Existenz 
bestimmter Orthogonalsysteme benutzt, wobei jedoch die Methoden, solche 
Orthogonalsysteme explizit herzustellen, sehr unübersichtlich sind. Anderer- 
seits gibt es eine sehr übersichtliche Methode, ein vollständiges linear unab- 
hängiges System von Funktionen zu bilden, welches jedoch nicht zu einem 
Orthogonalsystem führt. Diese Methode benutzt die Symmetrisierungs- bzw. 
Antisymmetrisierungsoperatoren bestimmter Youngscher Tableau. 


Der Zweck der vorliegenden Arbeit ist, zu zeigen, daß die gleichzeitige 
Benutzung des symmetrischen und antisymmetrischen Funktionssystems auf 
verhältnismäßig einfache Matrizen führt, mit deren Hilfe die totale Wellen- 
funktion hergestellt werden kann. Außerdem wird gezeigt, daß unter gewissen 
Bedingungen diese beiden Systeme ein Bi-Orthogonalsystem bilden, welches 
all das leistet, was man von einem Orthogonalsystem erwarten kann. Es wird 
vermutet, daß die Bedingungen immer erfüllbar sind, wenn dies bisher auch nur 
für Systeme bis zu 6 Nukleonen demonstriert werden konnte. 


Die Bi-Orthogonalität des Systems hängt von der Auswahl der Tableau 
ab. Es läßt sich zeigen, daß die übliche Auswahl der Tableau (wachsende Zahlen 
in Zeilen und Kolonnen) nur bis zu Systemen von 4 Nukleonen bi-orthogonale 
Systeme liefert. Für 5 und 6 Nukleonen führt jedoch eine andere Auswahl 
der Tableau zum Ziel. 
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en “= er . Die Youngschen Tableau und zugehörigen Operatoren 


Wir betrachten einen Symmetriecharakter [9,, 93, ---, Jn] = {hy, Ra, 
mit dem N Rahmen. Der Charakter [3, 2] = {2, 2,1} z.B. hat den 


Rahmen 


Der zugehörige Anti-Charakter ist definiert 


Im obigen Beispiel ist der Anti-Charakter {3, 2} = [2, 


durch [Ay, hs, An) 


2,1] mit dem 


Rahmen 


| 


Eine bestimmte Zahlen 1, ‚k auf den Rahmen, sei 
ein Tableau 7. Jedem Tableau 7 entspricht ein umgeklapptes Tableau 7 des 
Anti-Charakters, welches wir auch das Anti-Tableau nennen werden. 

Ein Tableau 7 werde durch eine Permutation pr- r in das Tableau 7 


überführt. 


5 = pr'r = prr T’. (1) 

Es seien p(T), p’(7’)..., Permutationen, welche nur Zahlen in einer Zeile 
des Tableau 7’ vertauschen und ¢, q(T), ex q'(T) . . . Permutationen, welche 


nur Zahlen innerhalb einer Kolonne vertauschen. Hier sei aul al 


+ 1 wenn é,q eine gerade Permutation una a 
wi |- 1 wenn &,g eine ungerade Permutation. 
Für jedes Tableau gibt es g verschiedene Permutation p(T) und h ver- 
schiedene q(T), wo 
h=h,!h,!...,h,,! 


Offenbar ist p(7)-! auch eine Funktion, welche zu den p(7') gehört und ent- 


sprechend für g(7’). niet 
Für das Anti-Tableau 7 gilt 


Wir definieren die Operatoren P, Q: 

Q(T) = Sq(7). 


(4) 
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Des weiteren gelten die Formeln ee 
pry P(T) = P(T’)prr prrP(T)= 
prr AT) = Q(T") prr Q(T) = Q(T") prr (7) 


p(T) P(T) = P(T) p(T) = P(T) p(T) Ps(T) = Ps(T) p(T) = e, P*(7) 
q(T) Q(T) = Q(7) (7) = Q(T) @(T) = Q(T) = & O (7). 
Das Auftreten oder Nichtauftreten des Faktors ¢ in diesen Gleichungen macht 


man sich leicht klar, wenn man die identischen Permutationen in der Summe 
für P oder Q betrachtet. 


— Wir betrachten jetzt einen wichtigen Zusammenhang zwischen zwei ver- 
schiedenen Tableau 7 und 7’. Wir unterscheiden zwei Fälle: 

a) Wir sagen 7’ sei zu 7 orthogonal, wenn es mindestens ein Zahlenpaar 
i, j gibt, welches im Tableau 7 in einer Zeile steht und im Tableau 7” in einer 
Kolonne. Das heißt, es gibt ein p(7) und ein g(7”), so daß 


p(T)=—qg(T)=P,, (Definition für 7’ | T) 
wo P,, die Vertauschung von i und j bedeutet. Es folgt 
P,, AT)=— Q(T’) P,, Q(T’) = Q(T"). 

b) Es gibt kein Zahlenpaar, welches obiger Bedingung geniigt. Dann stehen 
alle Zahlen, welche in 7’ in einer Kolonne stehen, in T in verschiedenen Zeilen. 
Wir können daher den Übergang von T zu 7’ vollziehen, indem wir erst die 
Elemente in den Zeilen von T so umordnen, daß sie in die richtige Kolonne 
kommen und dann die Elemente in den Kolonnen umordnen. Es folgt, daß 
ein g(T’) und ein p(T) existiert, so daß 

pr = &4(T') p(T) (wenn 7” nicht | 7’). 
Wir definieren zwei Permutationen ¢(7'), ?(7’) durch die Gleichungen 


q(T) = prr Q(T") pr = p(T) q(T") p(T) 


(8) 


AD UM 
pr kann also auch in der Form : 


pre = P(T)4(T)= eq 
dargestellt werden. Wir betrachten jetzt die Elemente einer Zeile aus dem 
Tableau 7. Die Operation g(7) verteilt diese Elemente auf verschiedene 
Zeilen, ändert aber nicht ihre Spaltenzugehörigkeit, so daß alle Elemente 
sich in- verschiedenen Spalten befinden. Die Operation p(T) permutiert diese 
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Elemente, läßt sie also alle in verschiedenen Spalten. Es folgt,da8 7" nicht | 7. 
Damit ist der Zirkel geschlossen und es folgt, daß sämtliche Formulierungen 
identisch sind: 
T’ nicht | T identisch mit pr r = T’) p(T 
x identisch mit pr r = & p(T) q(T) 
identisch mit pr 7 =e, p(T")q(7"). ; 
Es sei bemerkt, daß die Beziehung der Orthogonalität eine einseitige Beziehung e 
ist. Wenn 7” | T folgt nicht unbedingt, daß auch 7 | 7’. Wenn beide Be 


Es seien f, und /, zwei willkürliche Funktionen, welche von den Koordinaten BLS 
der k Nukleonen 1, 2... abhängen, auf welche die Permutationsoperatoren eae 
wirken. Diese Koordinaten können z. B. die Ortskoordinaten sein, oder die 
Orts-Spin-Koordinaten, oder Spin-Ladung-Koordinaten usw. Das Matrix- 
element 


(fy | fe) = fe dey ++ de, 


ist so definiert, daß das Integral über dieselben Koordinaten zu = 
ist. Für Spin- oder Ladungskoordinaten ist das Integral natürlich durch die 
entsprechende Summe zu ersetzen. 


Wir bilden 
(Q(2") fy | 

Vertauschung der Integration über die Variablen i und j ändert TER 
Integral. Wir erhalten 


(Q(7") h | P(T) fe) Q(T") | P;; P(T) 


Im Fall (a) folgt mit (8), daß das Matrixelement verschwindet. 
Die Operatoren Q und P sind hermitisch (siehe unten). Es folgt 


(Q(2") fy | fa) = Q(2") fa) = (PCT) Q(T") hy | fa) = 0 a) 
Da dies für jedes willkürliches Paar von Funktionen gilt, folgt 
Q(T’) P(T) = P(T)Q(T')=0 wenn 7” | T. 


Es sei beachtet, daß die Vertauschungsoperatoren hermitisch wen, aber ova ae 
die anderen Permutationen. Es gilt 


(hl Pi; fe) = (Pish 
(h | Pile) =(p | fe) 


= 
wie sich leicht beweisen läßt. Hieraus folgt die Hermitizität der Operatoren 
P und Q. 


4. Die Normalformen 


angewendet auf eine Funktion F stellt eine symmetrische Normalform = % 
Tableau 7, her 
R= PQAF= 


2 ann. Physik, 7. Folge, Bd. 9 


12 
| 
ziehungen erfüllt sind, sagen wir, daß 7’ und 7 zueinander orthogonal sind. | 
b Um nur die erste Beziehung auszudrücken, sagen wir 7” ist zu 7’ orthogonal. 
ht 
me 
er- 
ar 
| 
(12) 
| 
m 
se 
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Aus dieser Form können wir durch Permutation der Elemente andere Normal- 
formen erzeugen. Wir zählen nur die Normalformen, welche rg ver- 
schieden sind +) 

wo wir zunächst willkürlich für jedes » eines der g Tableau, welches zu derselben 
symmetrischen Form führen, herausgegriffen und mit T, eens haben. 


Pv» ist definiert durch die Gleichung ” 

Aus (2) folgt 


8 8 N - 
Pv = Pou Pui Pi = Pou Pu- (15) 


Mit Hilfe der Formeln (7) erhält man die alternative Schreibweise 
= P,Q, Pprı F 
Zwischen diesen n Funktionen bestehen lineare Beziehungen der Form 


wo die j alle Elemente einer Zeile sind, und i ein Element aus einer folgenden 
Zeile. Die Anzahl der linear unabhängigen Funktionen ist gleich dem Ent- 
artungsgrad f des Charakters. Wir numerieren die Funktionen g, so, daß die 
ersten f Funktionen ein linear unabhängiges System bilden 


(16) 


(17) 


Px = Pui Pi 


Die übrigen Funktionen können dann als Linearkombinationen der 9x geschrie- 


ben werden. 
J, 


8 
Pa = x Px . 
= 
Wir führen die ee ein, wobei wir für die Matrizen große Buch- 
staben benutzen. 


(18) 


= = (93). (19) 


Der Index 1 bzw. 2 deutet an, daB die Indizes der Elemente in der Spalte 
von 1 bis f, bzw. von f + 1 bis » laufen. Entsprechend ist die Matrix 


Ay, (a,x) (20) 


1) Der Einfachheit halber sei angenommen, daß keine gleichlangen Zeilen in den 
Rahmen vorhanden sind, d.h. die Zahlen g, sind alle verschieden. 

Wenn zwei gleichlange Zeilen vorkommen, führen die Tableau, welche sich nur durch 
den Austausch der beiden Zeilen unterscheiden, zu denselben symmetrischen Normal- 
formen, so daß jeweils 2 g verschiedene Tableau zur selben Funktion 9} führen. Die allge- 
meine Theorie läßt sich jedoch sehr leicht auf diesen Fall übertragen, indem wie die x 
der Funktionen 9) in 2 identische Unterserien » = 1,2,...n/2, und»=n/2 +1, 


einteilen. Wir verlangen dann, daß sämtliche im weitesen durchgeführten Ah 
in jeder Unterserie separat durchzuführen sind, so daß keine Linear-Kombinierung 
zwischen den Unterserien erfolgt. Das heißt die im folgenden auftretenden Matrizen sind 
in reduzierter Form anzunehmen, als eine sich entlang der Diagonale wiederholende Unter- 
matrix. Dann können die beiden Unterserien formal als verschieden betrachtet werden. 
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zu verstehen. Dann gilt 

Wir führen jetzt antisymmetrische Normalformen ein durch die Definition a 


Gx = % = % Per APi= PM Pi- U) 


Diese lassen sich als Linearkombinationen der symmetrischen Funktionen dar- = 
stellen 


Haar 


Die Funktionen g; hängen davon ab, welche Tableau 7’, ausgewählt worden B a 
sind. Wir werden sehen, daß die antisymmetrischen Funktionen ebenfalls ein 
linear unabhängiges System bilden. Dann lassen sich die symmetrischen Funk- 
tionen als Linearkombinationen der antisymmetrischen darstellen und die es . 
zu B,, inverse Matrix existiert. 


onl, iss 


Wir betrachten jetzt die Matrizen are he 
i 
En = Ey = (es) (25) 
Mit (16) und (22) folgt AR 


pF) 
Wenn T, zu T, orthogonal ist, folgt aus (11) 

(wenn T, | T,). (28) 
Wenn dies nicht der Fall ist, gilt (9) FEN 


= & (wenn Tanicht | 7,). 
Also mit (15) 


Po = PUT) Ga = Gi 2) 


Denn qj ist symmetrisch für die Nukleonen, welche durch 


werden. wid 


Aus (26) folgt daher mit (12) 

Nun ist 9} antisymmetrisch für die Nukleonen, welche durch &, g(7,) per- = 
mutiert werden. Es folgt, daß iodal hug 


al- 

er- 

14) 

en. 

15) 

16) | 
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Vertauschung der Reihenfolge der Integration führt mit (14) und (22) zu dem 
Resultat 
a 8 a 8 
&q (Pra Pia 73) = | gi) 
Wir können die Ausgangsfunktion so normieren, daß e,, = 1. Dann folgt 


e,3=0 wenn T, | T, 


= 

Es sei hier bemerkt, daB wir von der Tatsache, daB Y und Yi aus derselben 
Ausgangsfunktion F hergestellt worden sind, nur zum Zweck der Normali- 
sierung Gebrauch gemacht haben. 


Wenn also 94 (j), 9,(i) aus zwei verschiedenen Funktionen F(j) und F (i) 
entstanden sind, folgt sofort 


(90) | Pi (*)) = ei. 
Aus der Definitionsgleichung (18) für die a,, folgt durch Multiplikation mit 
(31) 
x=1 


Für festes x bildet (31) ein inhomogenes lineares Gleichungssystem von f 
Gleichungen für die f Unbekannten a,,. Dieses ist lösbar, wenn die Deter- 
minante 


leal = |By| A 0. (32) 


Dies ist also eine notwendige Bedingung dafiir, daB das ausgewahlte System 
linear unabhangig ist. Es ist auch eine hinreichende Bedingung, denn wenn 
eine Linearbeziehung 


=agi=0 
A 


und dieses Gleichungssystem hat nur dann eine Lösung, wenn die Determinante 
verschwindet. 


Die Gleichung (31) kann auch in Matrizenform geschrieben und vr 


Ey, = Ay Ey 


Hatten wir erst die antisymmetrischen und dann die symmetrischen Normal- 
formen gebildet und dabei auch die Tableau 7, benutzt, so hätten wir für 
E,, die umgeklappte Matrix erhalten; dies andert nichts an dem Wert der Deter- 
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minante. Hieraus folgt, daß die Tableau 7, auch für die antisymmetrischn 
Funktionen ein linear unabhängiges System liefern. Es folgt, daß die zuder 
in (24) eingeführten Matrix B,, inverse Matrix existiert. 
Für die antisymmetrischen Funktionen folgt aus der Definition (22) und (5) 
Nun sind nach Giiekang (29) und (30) für alle Funktionen g,, welche sich in _ 
der Form &,g(T,) Px darstellen lassen, die Koeffizienten e,, =, und für alle ü 


anderen Funktionen e,,, — 0. Die Summe (34) für @ kann also auch wie flgt 
geschrieben werden 


Oder 


Elimination von ®§ durch (21) ergibt 
| = (By + By Ay) 
Vergleich mit (24) ergibt mit (33) 


By = Ey + Ey Ay, = (1 + An Ay) 


By = Eu + Ex Ex Ei} 
q 6. Die Normalformen des Anti-Charakters 


Wir behandeln die Funktionen des Anti-Charakters in ähnlicher Weise mit 
geringfügigen Änderungen, die sich nur auf das Vorzeichen der Funktionen 
beziehen. Wir führen jedoch zuerst die antisymmetrischen Funktionen ein 7 
weil diese fiir das Anti-Tableau den symmetrischen Funktionen des Tableau 
entsprechen. 


vt = P(7,) = Pe(Ty (40) 
yy = €y1 Pri vi = Eu Pop Ya: (4) 

Die Linearbeziehungen zwischen diesen Funktionen sind dieselben wie bei den 

9, (vgl. Formel (17)), nur sind die Pi; durch — P,, zu ersetzen. Das—-Zeichen — 

ist jedoch in der Definition der Daher gelten dieselben 


Relationen und es folgt | % a3 
Die symmetrischen Funktionen seien definiert durch die Gleichungen 3 12 
A 
=P (T.) (48) 


wie sich leicht beweisen läßt. 
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Wir erhalten daher die Gleichung 


(45) 


Wir nehmen an, daß die Koordinaten, von denen die und y abhängen, 
sich zu dem Gesamtsystem der Koordinaten ergänzen. Zum Beispiel könnte 
gy die Ortsfunktion und y die Spin-Ladungsfunktion sein. Die aus diesen 
gebildete Gesamt-Wellenfunktion muß in allen Nukleonen antisymmetrisch 


Die folgende Funktion erfiillt diese Bedingung { 

Vi» vr 


Elimination von Y, und ®, mit Hilfe von (21) und (42) ergibt 


ya (1 Ay Ay) Bi. (47) 


Mit (38) folgt 
| 1= & By Di = Da By (48) 
Vi Vi Pyare 


Diese Gleichung voll ausgeschrieben, lautet rer) i 


- 


Da g und yw von verschiedenen Koordinaten abhängen, folgt 
(%\ x) ~ ( Gz qi) (pi! ens == T trace (Eu Eu =]. 
(51) 
Die Funktion x ist also auf 1 normiert. 
Wenn es möglich ist, die Tableau 7, so zu wählen, daß Z,, die Einheits- 
matrix ist, so erhält man die Wellenfunktion in der einfachen Form 

1 4 a 4 
In diesem Fall bilden die g% und g, ein vollständiges Bi-Orthogonalsystem, 
für welches 
= | wa) = (53) 
Die Bedingung hierfür ist, daß die Tableau 7’, , 7, zueinander orthogonal sind 
für alle x #9. 

Wenn £,, die Einheitsmatrix ist, ist u Bedingung (32) für ein linear unab- 
hängiges System erfüllt. 


oder mit (24) und (45) j 
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Es sei hier bemerkt, daß wir außer den Funktionen Px , welche den Vektor 

(45) | © bilden, auch Funktionen Pa definieren können, welche sich zum Vektor 
®% zusammensetzen. Die Gesamtheit dieser Funktionen ist allerdings im allge- 
meinen nicht gleich der Gesamtheit aller nicht identischen Funktionen, welche 

sich aus g% herstellen lassen, wie das für die symmetrischen Normalformen der 

en, | Fall ist, da hierfür ein anderes System der Tableau 7, ausgesucht werden 


nnte | müßte. pling 
risch sh 
woraus sich die Gleichungen 
und ‘T 
ur 
46 i 
(46) 0 bas 


ergeben. 
Unter Berücksichtigung der Tatsache, daß die hier auftretenden Matrizen 


reell sind, (also E+ = E), läßt sich folgende Beziehung ableiten 
welche für die Überprüfung der Richtigkeit der Matrizen nützlich ist. 


(48) 8. Die Matrizen für die Charaktere von 5 Nukleonen 


= 


Als Beispiel seien die Charaktere für 5 Teilchen betrachtet. Der vollsymme- 
trische Charakter [5] ist trivial. 
(49) Für den Charakter [4,1] ergibt sich folgende Tabelle: 


v | Tableau | e 


(50) 


vi 


3 


=1, 
(51) 


te 


eits- 


(52) 
tem, 
(53) 
sind 
nab- 


% 


(54) 
(55) 
0) 
x | 223% 0 1 0 0 
4/13 4 5 0 0 0 I-1 


Far den Charakter [3, 2] ergibt sich, daß die übliche Auswahl der Tableau 
nicht geeignet ist. 
1 2 3 ah ‘ 
4 5 2 4 


sind nämlich nicht orthogonal zueinander, da es kein Paar gibt, welches im 
_ ersten Tableau in einer Spalte steht und im zweiten in einer Zeile. 


7 Die in der folgenden Tabelle getroffene Auswahl fiihrt jedoch zum Ziel. 


yeni. 
w=] 12 8 1 0 0 0 0 ‘nih 
2 1 0 0 
3 ri 0 1 0 0 N 
4 | : , 1 | 0 ) 1 ( a a 
( 0 ) 
2 6 | 4 
5 245 0 0 0 0 1 AA 
eh @ 
| 85 
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34 
24 
y 345 1 0 0 
21 
10 324 |-ı 1 —ı 0 
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Charakter [3, 1, 1] ergibt sich mit der üblichen Auswahl | 


| Tableau 


3 | 


| 


2 
1 
5 
2 
1 
4 
2 
1 
3 
3 
1 
2 


Ar, 


( 


1 
Aw 
Ei: 


bleau| Für den 
| 
| | 
vine 2 2 4 0 1 0 0 0 0 | 0 0 0 0 
Ziel. | 
é 2 5 0 0 1 0 0 v0 0 0 0 0 
| 
4 5 4 0 0 0 1 0 0 |—1 0 0 0 : 
{ 35 0 0 0 0 1 0 ae 0 0 
N 
| 
6 4 5 0 0 0 0 0 1| 0 0 —1 ] 
| 
‘8 35 1 0 —1 0 —1 0 0 1 0 0 
9 45 0 1 1 0 © 0 
10 4 5 0 0 0 1 1 1 —-1 —1 
| 
1 3 1 0 1 
B 1 | 3 0 1 -1 2, = 
n° | p 3 1 
Die übrigen Symmetriecharaktere sind Anti-Charaktere der obigen und daher 
im Prinzip schon eingeschlossen. 
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9. Bi-orthogonale Systeme für die Symmetriecharaktere von 6 Nukleonen D 


Der Charakter [5,1] 


v | Tableau C9 Cys 
x=1/123 4 5 1 0 0 
6 | 
=2/12346! 0 1 0 
| 5 | 
=8/12356/ 0 o 1 
| 4 
=4/1245 6] 0 0 0 
3 
=5 13 45 6] 0 0 0 
; 2 
a=6|2 3456|-1 -1 -1 
12111 
Bee. 


Der Charakter [3,2,1] 


Die Anzahl verschiedener Funktionen g, ist n = 60, der Entartungsgrad 
ist f = 16. Wir sparen uns die Mühe, die Matrix E,,, welche 44 Zeilen und 16 
Kolonnen enthält, aufzuschreiben. Die Matrix Z,, ist die Einheitsmatrix, wenn 
man folgendes System zueinander orthogonaler Tableau benutzt. 


132 312 162 364 162 142 345 365 

45 46 34 35 12 12 

6 5 5 6 6 4 
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3 4 35 16 64 62 65 45 42 62 
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K. 
Der Charakter [4,2] 


AD 
> = > o © 
—) > o > = 
o o = > = 
— > o > = 
s 8 a4 © > 
oO x > > 
o > -) — o = 
- > > = = > _ > 
~+~ 0 © + © © 
Il 
N 


27 
v | Tableau | 2 3 ey 4 & 7 yg %9 
| 
h 
12 
1| 1346 0 0 1 1 0 
| 1 5 
13 28 45 —1 —-1 0 0 —1 —ı 0 0 0. 7 
16 4 1456 0-1 -1 —ı 0 0 —1 0 
| 23 
ee | 
2 4 2 1 1 1 -1 -1 
A —1 1 1 1 —1 2 
—1 1 2 2 —1 —2 4 —ı 1 Er 
=1 —§ 1 2 4 2 
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Der Charakter [4,1,1] 
Tableau e, 1 e, 2 e, 3 e, 4 e, 5 ey 6 e, 7 Cyg e, 9 e, 10 
6 
4 
6 
4 
5 
3 
6 
3 
5 
3 
4 
2 
6 
9 
5 
2 
4 
2 
3 
1 
2 
1 
3 
1 
4 
1 
5 
1 
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1 3 1 0 —1 1 0 —1 
0 kei de 1 1 0 3 ] 1 | 0 0 
B 1 0 1 ] 3 1 0 1 0 
0 
Der Charakter [3,3] 
128 3: 0 0 0 
0 45 6 
12 4 0 1 0 0 0 
3 5 6 
0 12 8 0 0 1 0 0 
3 4 6 
13 4 0 0 0 1 0 
0 265 
_ 35% 0 0 0 0 1 
245 
2% —1 -1 —1 0 0 
345 
1 19353 0 0 —1 —1 
246 
145 0 —ı1 1 0 1 
| 236 
28 & 
l 1 5 6 1 1 0 1 0 
234 
0 4 9 1 2) 1 
= 
0 1 2 4 -1 2 4 
> 
24 
Rossendorf bei Dresden, Zentralinstitut für Kernphysik, Bereich _ 
Theorie. 
22 (ase Bei der Redaktion eingegangen am 9. Juni 1961. “a! 
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Sekundärelektronenemission von Caesium—Germanium- 
Schichten verschiedener Zusammensetzung 


Von G. Appelt 


Inhaltsübersicht 


Es werden Messungen der Ausbeute, der Temperaturabhängigkeit der Aus- 
beute sowie der Energieverteilung der Sekundärelektronen (S) an Cs—Ge- 
Schichten verschiedener Zusammensetzung durchgeführt. Die Schichten waren 
einerseits aus amorphen, andererseits aus kristallinen Ge-Schichten herge- 
stellt. Aus amorphen Ge-Schichten bildet sich durch Einwirkung von Cs-Dampf 
die Verbindung CsGe, während bei kristallinen Ge-Schichten das Cs wahr- 
scheinlich nur als Adsorptionsschicht vorliegt und das Grundgitter keine Ände- 
rung erfährt. Cs—Ge-Schichten beider Typen weisen hohe Ausbeutewerte auf. 
Es ergeben sich Anzeichen dafür, daß die hohen Ausbeuten wesentlich durch 
eine Erniedrigung der Potentialschwelle an der Oberfläche verursacht werden, 
während die Austrittstiefe der S durch die Cs-Behandlung der Schichten un- 


Cs—Sb- und Cs—Bi-Schichten als typische Vertreter der A!—BVY-Verbin- 
dungen wurden hinsichtlich ihrer Sekundäremissionseigenschaften vor kurzem 
systematisch untersucht!). Für Cs—Sb bestätigten die Messungen bezüglich der 
hohen Ausbeutewerte und der Energieverteilung der S frühere Ergebnisse 
anderer Autoren. Sie zeigten ferner, daß Cs—Bi hinsichtlich der Sekundärelek- 
tronenemission (SEE) dem System Cs—Sb außerordentlich ähnlich ist. Darüber 
hinaus ergab sich aus den Messungen ein Hinweis darauf, daß die unerwartet 
hohen Ausbeutewerte dieser Verbindungen in erster Linie auf eine niedrige 
Potentialschwelle an der Oberfläche zurückzuführen sind, eine Eigenschaft, die 
alle A'—BY-Verbindungen auszuzeichnen scheint. In dieser Hinsicht unter- 
scheiden sich die A'—BY-Verbindungen von einigen anderen Stoffgruppen mit 
hohen Ausbeutewerten, nämlich den Alkalihalogeniden und Erdalkalioxyden, 
bei denen die hohe Ausbeute in erster Linie durch eine große Austrittstiefe der 
S hervorgerufen wird. 

Angesichts dieses Ergebnisses ist es von großem Interesse, auch andere 
Stoffe mit niedriger Austrittsarbeit daraufhin zu untersuchen, ob sie einen 
hohen Ausbeutekoeffizienten besitzen und ob sich auch hier wie beiden Al—BV- 
Verbindungen ein Hinweis darauf zeigt, daß die hohen Ausbeutewerte auf die 
niedrige Austrittsarbeit zurückzuführen sind. 


1) G. Appelt u. O. Hachenberg, Ann. Physik '7) 6, 67 (1960). 
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Es ist bekannt?), daß Cs—Ge im sichtbaren Spektralbereich einen gut meß- 
baren äußeren Photoeffekt besitzt und daher eine niedrige Austrittsarbeit hat. 
Dieses AI-BIV.System wurde daher für die Untersuchungen ausgewählt, die 
sich wie bei den A'—BY-Verbindungen!) auf Messungen der Ausbeute in Ab- 
hängigkeit von der Energie der Primärelektronen, auf Messungen der Tempe- 
1- raturabhängigkeit der Ausbeute und auf die Ermittlung der Energieverteilung 
der S erstreckten, wobei die Schichten jeweils bei verschiedenen Zusammen- 


setzungen gemessen werden sollten. al. 


2. Herstellung der Schichten 


Die Cs—Ge-Schichten wurden technologisch in der gleichen Weise wie die 
Cs—Sb- und Cs—Bi-Schichten!) hergestellt. Weitere Einzelheiten darüber fin- 
den sich in®). Auf einen Objektträger wurde zunächst eine Ge-Schicht !) auf- 
us- | gedampft, welche man dann Cs-Dampf;) 


Ge- aussetzte. Die Temperatur während dieses 

ren | Formierprozesses betrug etwa 140 °C. Es a 

ge- | ist bekannt, daß man beim Aufdampfen kristallin 
npf | von Ge je nach der Temperatur des Ob- §& @ te ST og I) 
jektträgers hoch- oder niederohmige Ge- 

de- ? Schichten erhält. Nach Untersuchungen = “ 

| von G. Haas*) und S. A. Semiletow’) ~% 


rch | handelt es sich dabei um amorphes und BR 


-7 
len, |} kristallines Ge. Die Grenztemperatur, 0 


un- f oberhalb der die Aufdampfschichten kri- 077 N \y / 
stallin werden, beträgt 370 °C. Für die vor- 


| 


liegenden Untersuchungen wurde bei der JE 
Präparation der Schichten sowohl von E lem... 
amorphen Ge-Schichten als auch von kri- Abb. 1. Verlauf des Leitwertes während 


N der Formierung von amorphen und 
bin- | stallinen ausgegangen. kristallinen Ge-Schichten in Cs- 
zem Die Kontrolle des Formierprozesses Dampf 

der | erfolgte wie bei den früheren Untersu- 

isse | chungen!) durch Messung der Leitfähigkeit. Abb. 1 zeigt den zeitlichen Ver- 
lek- | lauf des Leitwertes während der Formierung für amorphe und kristalline Ge- 
iber | Schichten. Beide Kurven haben qualitativ den gleichen Verlauf: nach einem 
rtet | Anstieg des Leitwertes folgte in Minimum (I) und darauf wieder ein monotoner 
rige | Anstieg bis zu einem Leitwert (II), der offenbar einem Gleichgewichtszustand 
‚die | zwischen dem Cs-Dampf und der Schicht entspricht. Eine weitere Fortdauer 
iter- § der Formierung ruft dann keine Leitwertsänderungen mehr hervor. Während 


Jen 2) N. Schaetti u. W. Baumgartner, Helv. Phys. Acta 28, 524 (1950). 

a ; 8) G. Appelt, Diplomarbeit, TH Dresden 1958. 

si 4) Als Ausgangsmaterial wurde einkristallines n-leitendes Ge benutzt. Es war spek- 
tralrein bis auf Spuren von Ca, die aber vermutlich aus den bei der Analyse benutzten 
dere | Gefäßen stammten. 

inen 5) Das Cs lag in metallischer Form in Ampullen vor. Die Spektralanalyse zeigte 
BY Spuren einer Reihe von Verunreinigungen (Al, B, Ca, Cu, Fe, Mg, Na, Si, Sn, Ti), die nach 
~~, # Erfahrungen von C. Kunze (Vak.-Technik 8, 168 ( 1959)) zum großen Teil dem Glas ent- 
f die | stammen. Infolge der Destillation während des Formierprozesses wirkt Cs sicher in sehr 
reiner Form auf die Ge-Schichten ein. 


6) G. Haas, Physic. Rev. 72, 174 (1947). Wee ts A 
S. A. Semiletow, Kristallografija 1, 542 (1956). wr 
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sich bei kristallinen Ge-Schichten der Leitwert nur um wenig mehr als eine 
Zehnerpotenz ändert, erstrecken sich die Leitwertsänderungen amorpher Ge- 
Schichten über fünf Zehnerpotenzen. Auch ist das Leitwertsminimum hier 
wesentlich stärker ausgeprägt als bei den kristallinen Ge-Schichten. 

Ein Leitfähigkeitsminimum läßt bei einem Zweistoffsystem immer eine 
stöchiometrische Verbindung vermuten, und es soll deshalb im folgenden ge- 
prüft werden, ob die Minima in den Kurven der Abb. 1 einer stöchiometrischen 
Zusammensetzung zugeordnet werden können oder nicht. 

Das System Cs—Ge wurde bereits 1948 von E. Hohmann’) im Rahmen 
einer zusammenfassenden Bearbeitung der Alkaligermanide und -silicide in 
Form kompakten Materials untersucht. E. Hohmann stellte eine stabile 
Verbindung mit der Formel CsGe fest und vermutete eine Zwischenverbindung 


20 + 7 
Sur 
3 = 
| if \ Ge (aniston) 
02 a i \ Cs-Ge (I) 
Cs-Ge (I) 


4% 05 0607 mn 6 18 19 20 212 
Wellenlänge A /44——= 


Abb. 2. Spektrale Durchlässigkeit von Ge- und Cs—Ge-Schichten auf Grundlage amorphen 
und kristallinen Germaniums 


mit der Zusammensetzung CsGe,. Die Formel CsGe, ist inzwischen von R. 
Schafer?) im Jahr 1957 auf Grund eines Analogieschlussesals sehr wahrschein- 
lich bestätigt worden. Schäfer konnte ferner nachweisen, daß nur eine 
Zwischenverbindung vorhanden ist und diese sich bei ungefähr 500 °C zersetzt. 
Ihre Bindungsenergie ist offenbar sehr gering, denn bei der thermischen Zer- 
setzung von CsGe bildeten sich beim Abdestillieren von Cs gleichzeitig CsGe, 
und reines Ge, wie röntgenographische Untersuchungen zeigten. Ferner hat 
CsGe nach R. Schäfer kubische Struktur mitz = 32 mol/Zelle und der Gitter- 
konstanten a = 13,64 + 0,02kX und eine Dichte o = 4,28 g/em? bei 25 °C. 
Um nun zu prüfen, ob sich in den Leitwertsminima eine der existierenden 
Verbindungen ausbildet, wurden Messungen der optischen Durchlässigkeit 
durchgeführt, und zwar vor Beginn des Formierprozesses an reinem amorphen 
bzw. kristallinen Ge, dann in den dem Leitfähigkeitsminimum (I) und etwa dem 
Endzustand (II) entsprechenden Zusammensetzungen. Nach Abb. 2 findet man 
für amorphes Ge nach der Formierung eine Durchlässigkeitskurve mit ver- 
änderter Absorptionskante bei A = 0,7 u, während sich die optische Durch- 
lässigkeit bei kristallinem Germanium durch das Formieren nicht ändert. 


8) E. Hohmann, Z. anorg. allg. Chem. 257, 113 (1948) 
®) R. Schäfer, Diss. Münster 1957. 
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Daraus kann man schließen, daß das Leitwertsminimum bei Schichten, 
die aus amorphen Ge-Schichten erzeugt wurden, einer stöchiometrischen 
Cs—Ge- Vi ve entspricht. Höchstwahrscheinlich handelt es sich den Er- 


amorphen Ge- Schichten um die Verbindung CsGe, während CsGe, wegen seiner . 
geringen Stabilität nicht auftritt. Für die Verbindung CsGe spricht ferner das 
Ergebnis einer in?) durchgeführten Abschätzung des Quellfaktors (dem Ver- 


hältnis der Dicke der ausformierten Schicht zur Dicke der ursprünglichen 
reinen Ge-Schicht) aus den Durchlässigkeitsmessungen der Abb. 2. Es ergibt S 
sich ein Zahlenwert des Quellfaktors, der demjenigen sehr ähnlich ist, der sich 
aus dem Quotienten der Molvolumina M/o von CsGe und Ge errechnen läßt. 
Eine CsGe-Schicht ist übrigens im Aussehen von einer Ge-Schicht kaum zu a 7 
unterscheiden. Im Vergleich mit der in der Durchsicht kräftig braun gefärbten we 
Ge-Schicht erscheint CsGe jedoch etwas rotbraun. mt 


Bei den formierten kristallinen Ge-Schichten dagegen tritt im Leitfahig- 
keitsminimum keine stöchiometrische Verbindung auf. Wir haben es hier 
sicher mit einer Ge-Schicht zu tun, an deren Oberflache Cs adsorbiert ist. 

Auf Grund des Verlaufes der Leitfähigkeitskurven wird nun der Versuch _ 
unternommen, Schichten mit annähernd reproduzierbaren Eigenschaften - 
durch zu erzeugen, daß man den Formierprozeß an ausgezeichneten Punkten 
der Leitfähigkeitskurven unterbricht. Für die SEE-Messungen wurden folgende  —s_— 
Punkte ausgewählt: 


1. reine Ge-Schicht pret PP Ge (amorph) Ge (kristallin) 

2. Formierzustand I CsGe Ge + Cs 
(Leitwertsminimum) 

3. Formierzustand II CsGe + Cs Ge + Cs 


(Endwert des Leitwertes) 


Wahrscheinlich ist nach dem eingeschlagenen Verfahren nur bei amorphen 
Ge-Schichten im Leitfähigkeitsminimum, also beim Auftreten der stöchiometri- 


Ge- Schichten in beiden Formierzuständen brauchen .gleiche Leitfähigkeite. ve 
eigenschaften nicht unbedingt auch gleiche Zusammensetzungen zu bedingen. — 
Diesem Umstand muß bei der Beurteilung der Versuchsergebnisse Rechnung 


getragen werden. a 


3. Die Versuchsréhre ae 


Diefür die SEE-Messungen an den Cs—Ge-Schichten verwendeten Versuchs- => 
röhren waren von der gleichen Bauart wie die zur Untersuchung der Cs—Sb- 
und Cs—Bi-Schichten benutzten Röhren!). Der einzige Unterschied bestand 
im zusätzlichen Einbau eines kleinen Ofens, mit dem der Objektträger während 
der Ge-Bedampfung zur Erzielung kristalliner Schichten geheizt werden 
konnte. Der technologische Ablauf des Pumpprozesses war dem bereits in!) a 
beschriebenen ähnlich. Nach der Ge-Bedampfung konnten die Versuchsréhren _ 
bei einem Druck von etwa 5 - 10-7? Torr von der Pumpe abgeschmolzen werden. 


4. Meßmethoden und Durchführung der Messungen 


Die Meßmethodik war im Prinzip die gleiche wie in!). Die Schaltung erfuhr 
jedoch im Hinblick auf die Messungen der Energieverteilung der S 


3 Ann. Physik, 7. Folge, Bd. 9 
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eine entscheidende Verbesserung insofern, als die Differentiation der Gegen- 
spannungskurve nunmehr auf elektrischem Wege durchgeführt wurde. Das 
neue Verfahren brachte eine erhebliche Steigerung der Meßgenauigkeit. Vor 
allem wurden die beim rechnerischen oder graphischen Differenzieren der 
Gegenspannungskurve entstehenden zufälligen Fehler vermieden. Dem Ver- 
fahren liegt das Prinzip zugrunde, der zwischen Probe und Auffangkäfig liegen- 
den Gleichspannung U, noch eine Wechselspannung U, hinreichend kleiner 
Amplitude in Reihe zu schalten. Dann muß auch der über den Auffangkäfig 
fließende Strom eine Wechselkomponente enthalten, die der Steigung der 
Gegenspannungskurve im Punkte U, und damit der Energieverteilung der 8 
proportional ist. Die Amplitude der Wechselspannung U, bestimmt die Energie- 
auflösung. 

Abb. 3 zeigt die verbesserte Versuchsschaltung. Gegenüber der in!) benutzten 
Schaltung ist in die Zuleitung zur Probe zusätzlich zur Gleichspannung U, die 
Wechselspannung U, eingeschaltet, und zwar ist die dem oberen Abgriff des 
Potentiometers R, entsprechende Spannung als U, wirksam. In die Zuleitung 
zum Auffangkäfig ist der Widerstand R, eingefügt. Die an ihm auftretende 
Wechselspannung wird verstärkt und an einem Röhrenvoltmeter RV zur An- 
zeige gebracht. Der über die Kapazität C,, zwischen Probe und Auffangkäfig 
fließende Wechselstrom wird dadurch unwirksam gemacht, daß C, , zusammen 
mit C, und R, eine Wheatstonesche Brücke bildet, bei deren Abgleich der 
im Nullzweig liegende Widerstand R, stromlos ist. 

Zur Messung wurde zunächst ohne Strahlstrom bei beliebiger RC-Genera- 
tor-Spannung das Potentiometer R, auf Minimum des Ausgangssignals ein- 
gestellt und danach die Ausgangsspannung des Generators so verändert, bis 
U, den vorgeschriebenen Wert hatte. Darauf wurde bei möglichst konstantem 
Strahlstrom die Spannung U, in Abhängigkeit von U,, gemessen. Da der Strahl. 
strom i, während der Dauer einer Messung nicht hinreichend konstant gehalten 
werden konnte, wurden die Meßwerte U, rechnerisch auf einen konstanten 
Wert des Strahlstromes bezogen. Bei den Messungen betrug i, = 2: 10-7 A, 
U,err = 100 mV bei f = 36 Hz und der auf U, bezogene hi a ew 1mV. 
SchlieBlich wurde der MeBvorgang noch automatisiert, es wurde die Gleich- 
spannung U, durch ein motorbetriebenes Potentiometer zeitproportional ver- 
ändert und die Anzeige des Röhrenvoltmeters mit einem Kompensations- 
schreiber registriert. Während der damit erreichten erheblich kürzeren Meß- 
dauer blieb der Strahlstrom hinreichend konstant. Bei sonst unveränderten 
Parametern betrug die Meßfrequenz f = 130 Hz. Dadurch konnten Schwe- 
bungserscheinungen zwischen Meßfrequenz und Wechselrichterfrequenz des 
Kompensationsschreibers (50 Hz) vermieden werden. Einige Kurven der Ener- 
gieverteilung wurden nach diesem Verfahren gewonnen. 

Die gemessene Funktion U, = f(U,) muß noch zwei meßtechnisch begrün- 
deten Korrekturen unterworfen werden, ehe sie der Energieverteilung der $ 
proportional ist. Zunächst muß der Abszissenwert U, um den Spannungsabfall 
I: R, korrigiert werden (J, = Gleichstromkomponente des Auffangkäfig- 
Darüber hinaus ist die W des 


dem differentiellen Innenwiderstand der Röhre nur dann streng 


proportional, solange (so oh ‘ani > R,ist. Anderenfalls müssen auch die Ordinaten- 
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werte U, korrigiert werden. Unter den Bedingungen der vorliegenden Bektiing ; 
braucht die Abszissenkorrektur nicht durchgeführt zu werden, denn sie liegt _ 
höchstens in der Größenordnung von 0,1 V. Die Ordinatenkorrektur dagegen 
macht sich bemerkbar, allerdings nur bei größeren Ausgangssignalen U,. Bei _ 
U, < 0,4 V beträgt die Korrektur weniger als 2,5%, bei U, = 1V,2V bzw. 3V 
beträgt sie z. B. 6, 13 bzw. 20°, und zwar sind die der wahren Energieverteilung 
entsprechenden Ordinatenwerte um die angegebenen Prozentzahlen größer als 
die gemessenen. Die unkorrigierte Funktion U, = f(U,) ist also der Energie- 
verteilung der S ähnlich bis auf die Maxima, die etwas verflacht wiedergegeben 
werden. Bei den vorliegenden Messungen wurde die Ordinatenkorrektur bei Phy 
Ausgangsspannungen U, > 0,4 V durchgeführt. San 
Es sei noch bemerkt, daß eine ähnliche Schaltung zum unmittelbaren N 
Differenzieren der Gegenspannungskurve vor einigen Jahren von N. B. Gornij _ 
und L. M. Rachowitsch!P) veröffentlicht worden ist. Der grundsätzliche _ 
Unterschied zu der hier verwendeten Anordnung besteht in der Art der Be- 
seitigung von Störeinflüssen. Im Gegensatz zu der vorliegenden Anordnung 
mit vollständiger Abschirmung von Versuchsröhre und Meßleitungen arbeiteten - 
Gornij und Rachowitsch mit einer gegenphasigen Störkompensation, mit 
der sich der Storpegel höchstwahrscheinlich nicht in so großem Maße senken 7 5 
läßt. Daraus folgen einige Nachteile hinsichtlich Energieauflösung und Ab- DR a 
veichungen zwischen wahrer Energieverteilung und gemessener Funktion re 


U, = f(U,). Die letzteren lassen sich allerdings durch die bereits geschilderten 4: oe 
Korrekturen beseitigen. 


Mit der Schaltung in Abb. 3 können auch die Messungen der Ausbeute 


durchgeführt werden. Um bei diesen Messungen jedoch den Einfluß der Ter- 
22 
Heizbarferie 
| Batterie für 
RC-Generator Wehnelt- 
ZSL-AUSIONG sponnung 
Batterie für 
Platten- 
lie hts Abb. 3. Schaltung der Versuchsanordnung 


tiärelektronen zu unterdrücken, muß das Netz gegenüber dem Auffangkäfig 
eine negative Vorspannung erhalten. Die Verbindung zwischen Netz und Auf- 7 > 
fangkäfig in Abb. 3 wurde daher an der mit einem Kreuz gezeichneten Stelle in 


0) N. B. Gornij u. L. M. Rachowitsch, Shurn.exp. teor. Fis. (russ.) 26, 454 (1954). - = 
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getrennt und eine Batterie von 12,5 V zwischengeschaltet (Minuspol an der 
Netzelektrode). Beziiglich der Gleichstréme ist die Schaltung dann mit der 
in!) gezeigten Versuchsschaltung vollkommen identisch. Als Ausbeute berech. 
neten wir den Ausdruck 


Hierin bedeuten J}, J, den über die Probe fließenden Strom bei positivem 
bzw. negativem Hilfspotential (50 V) der Probe gegenüber dem Auffangkäfig, 
i, den Strom der echten S mit E < 50 eV, i, den Strom der reflektierten und 
rückdiffundierten Primärelektronen (P) mit Rn >50eV, i, den 
primären Strahlstrom, der bei den Ausbeutemessungen & 6- 10-8 A betrug 
Das entspricht bei dem verwendeten Strahlsystem einer primären Strahl. 


stromdichte von S 2 - 10-* A/em?. 


Zur Messung der Temperaturabhängigkeit der Ausbeute wurde 
die Versuchsröhre in einem Thermostaten aufgestellt, dessen Temperatur 
zwischen — 70 °C und 100 °C beliebig eingestellt werden konnte. Es wurde 
wie auch in!) die Ausbeute bei einer konstanten P-Energie £8 oberhalb des 
Ausbeutemaximums in Abhängigkeit von der Temperatur gemessen. 


5. Meßergebnisse 
a) Messungen der Ausbeute 


Die Ausbeutekurven für 
Cs—Ge-Schichten, die ausge- 
hend von amorphem Germa- 


nium prapariert wurden, sind 


& 


A 


Pie] in Abb. 4dargestellt. Amorphes 
ud DS Ge besitzt unseren Messungen 
sy ag zufolge eine Maximalausbeute 
on *| von 1,1, die bei EI = 600 

Fo — 700 eV erreicht wird. Für 

20 37 ~ CsGe im Formierzustand I be- 
en Cae trug Ömax bei allen gemessenen 


Proben etwa 4,5 (Kurven 2a 
bis c), während im Formierzu- 


I la stand II (CsGe + Cs) bei den 
f Fe einzelnen Proben stark abwei- 


RNR 


10 
0 ‘ed chende Ausbeutekurven mit 
1 2 3 36 Maximalausbeuten zwischen) 
Energie der Primärelektronen Ep °/ key —— und 6,9 erhalten wurden 
_ Abb. 4. SEE-Ausbeute von Cs—Ge-Schichten ver.  (Kurven3a—e). Die überein 
schiedener Zusammensetzung auf Grundlage amor- stimmenden Ausbeutekurven 


phen Germaniums. Kurven1: Ge(amorph);Kurven im Formierzustand I können 
CsGe Kur I) 3: +Cs dahingehend gedeutet werden, 

ormierzustand II); Kurve 4:nichtgenau bekannte 
Zusammensetzung in der Nähe des Leitwertmini- daß ur diesem Zustand die er 
mums. (Die Buchstaben a—e kennzeichnen 5 ver-  bile stöchiometrische Verbin- 
schiedene Proben) dung CsGe vorliegt. Im For- 
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mierzustand II dagegen haben wir CsGe-Schichten mit Cs-Überschuß, der bei 
den einzelnen Proben sehr verschieden sein kann. Während der Formierung 
durchläuft die Austrittsarbeit sicher ein Minimum, bei dem dann die SEE- > 
Ausbeute maximale Werte annimmt. Man erhält eine um so höhere Ausbeute, ae ; 
je naher der betreffende Zustand demjenigen mit minimaler Austrittsarbeit 
ist. Für dieses Verhalten spricht Kurve 4e mit der sehr hohen Maximalaus- — 3 
beute von 7,8, die an einer Cs—Ge-Probe von nicht genau bekannter Zusam- _ 
mensetzung in der Nähe des Leitwertsminimums (wahrscheinlich zwischen a 
Formierzustand I und II) gemessen wurde. Das Ausbeutemaximum aller mit 
Cs behandelten Schichten in Abb. 4 liegt bei 600 bis 800 eV P-Energie und _ 
damit also bei gleichen Werten ES „ wie bei den reinen Ge-Schichten. ER 
Abb. 5 zeigt die Ausbeutekurven der Schichten, die ausgehend von kristelli- DEI 
nem Germanium präpariert wurden. Für kristallines Germanium erhielten ir | 
Maximalausbeuten von 1,65 und 1,83 bei E2,, = 700 bzw. 1400 eV (Kurven a Be 


la u. b). Nach dem Formieren mit Cs a5 

betrug die Ausbeute der Schichten ' 8 m 

bmax 4— 8, wobei in den Formier- | 

zustanden I und II stark unterschied- 2 

liche Ausbeutekurven erhalten wurden. & 6 

Dies ist nicht verwunderlich, denn bei = 2» 

diesen Schichten brauchen, wie schon % R 
erwähnt, ähnliche Leitfähigkeitseigen- Ja 

schaften keinesfalls ähnliche SEE-Eigen- 2% 

schaften zu bedingen. Die Lage des Aus- 

beutemaximums blieb bei jeder Probe 

ungefähr bestehen mit Ausnahme der 2 ja i) 
Probe a, bei der sich das Maximum von ‘coe Sees oe 
Kurve 3a gegenüber dem der Kurven la Ib | 
und 2a zu etwas niedrigeren P-Energien 


Vergleicht man die Ergebnisse der Energie der Primörelektronen Ep /keVv— = 
Ausbeutemessungen mit denen anderer Abb. 5. SEE-Ausbeute von Cs—Ge- Pe 
Autoren, so findet man nur zum Teil eine 
Übereinstimmung. Für Ge erhielten H. maniums. es, Ge (kristallin); 
Gobrecht und F. Speer!!) dmax = 1,25 Kurven2: Ge + Cs(Formierzustand 
bei Z9,, = 500eV, ferner J.B. Johnson Kurven3: Ge+Cs(Formierzustand II). 
und K. G.McKay®) dmax = 1,15 bei (a und b verschiedene 
400 eV. Aus neuerer Zeit sind Messungen 
von W. G. Bolschow und W. K. Seleznew!3) sowie von W. G. Bolschow a bis 
u. Mitarb.) bekannt mit dmax = 1,05 bei ES, = 500 eV. Diese Zahlenwerte _ 
stimmen mit unseren Meßergebnissen für amorphe Ge-Schichten recht ong 
überein, wenn man bedenkt, daß die genannten Meßergebnisse anderer Autoren | co ty 
z.T. an kompaktem Material und auch mit etwas unterschiedlichen MeB- 
schaltungen gewonnen wurden. W.G. Bolschow u. Mitarb.) geben auch 


Er 


1) H. Gobrecht u. F. Speer, Z. Physik 135, 602 (1953). 
12) J. B. Johnson u. K. G. McKay, Physic. Rev. 98, 668 (1954). 
13) W. G. Bolschow u. W. K. Seleznew, Shurnal techn. Fis. (russ.) 26, 1657 (1956). — 
14) W. G. Bolschow, L. N. Dobrezow, A. A. Scharinow, T. W. Kratschino u. 
M. K. Repnikowa, Fisika twerd. tela 1, 1768 (1959). Si 
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eine Ausbeutekurve für Cs—Ge an mit 
Ömax = 1,8 bei ZZ» = 350 eV, offenbar an 
kompaktem Material gemessen, das nach 
kurzer Glühung bei 800 °C Cs-Dampfen bei 
einer Temperatur von 150 °C ausgesetzt 
wurde. Dieser Ausbeutewert liegt erheblich 
unter den von uns gemessenen Zahlenwerten. 


b) Messungen der Temperaturabhängigkeit 
der Ausbeute 

Die Messungen der Temperaturabhängig- 
keit der Ausbeute brachten zutage, daß sich 
Cs—Ge-Schichten in den verschiedenen Zu- 
ständen instabiler verhalten als Cs—Sb- und 
auch Cs—Bi-Schichten!). Die Meßpunkte 
streuten im allgemeinen so stark, daß eine 
Temperaturabhängigkeit der Ausbeute nicht 
mit Sicherheit nachgewiesen werden konnte. 


Ss 


energieEp > Epm Nur nach längerem Elektronenbeschuß 

(Zp = 1500 eV; i,» 107 A, 

t~ 1h) konnte an einer 

9 : | CsGe-Schicht die in Abb. 6 

Ni ‘ | dargestellte Kurve gemessen 

IN | werden. Es zeigt sich eine 

25 „il nur sehr schwache Abnahme 

der Ausbeute mit der Tem- 
N peratur. 

c) Messungen 


der Energieverteilung 
der Sekundärelektronen 


Die Messungen der Ener- 
gieverteilung der S sind in 
den Abb. 7 und 8 dargestellt. 


Energieverteilung (aut ip normiert) 


Die MeBkurven in Abb. 7 gel- 
ten fiir Cs—Ge-Schichten, die 
auf Grund von amorphem 


un 


Gegenspannung ——= 


Abb. 7. Energieverteilung der Sekundärelektronen von 
Cs—Ge-Schichten auf Grundlage amorphen Germani- 
ums, auf i, normiert. Kurven 1: Ge (amorph); Kur- 
ven 2: CsGe (Formierzustand I); Kurven 3: CsGe + Cs 
(Formierzustand II); Kurve 4: Schmiegungskurve (Er- 
klärung siehe Text); Kurve 5: Kurve 4 aus Abb. 8. 
(a—d kennzeichnen verschiedene Proben) ty 


Germanium präpariert wur: 
den. Sie wurden mit Hilfe 
des Kompensationsschreibers 
registriert und dann umge- 
zeichnet. In Abb. 8 sind die 
Kurven für die Schichten auf 
Grundlage kristallinen Ger- 
maniums dargestellt, die 
punktweise gemessen worden 
sind. Sämtliche Kurven sind 
auf konstanten Strahlstrom 
normiert und direkt in 
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Abhängigkeit von der Gegenspannung U, dargestellt. Das allmähliche Aus- Sr 


laufen der Kurven in Richtung negativer Spannungen U,, insbesonders die 


kleinen Nebenmaxima der Kurven la und 1bin Abb. 8, werden durch Ter- 


tiärelektronen verursacht und sind nicht reell. Man kann daher die in Rich- _ 
tung der negativen Abszissenachse liegenden Kurvenflanken durch eine 
Gerade annähern. Dem Schnittpunkt dieser Geraden mit der Abszissenachse _ 
ist die S-Energie E, = 0 zuzuordnen. Der durch die Geraden festgelegte Abszis- __ 


senabschnitt entspricht dem Kontaktpotential, der Differenz der Austritts- 


arbeiten von Auffangkäfig 
und Cs—Ge-Schicht. Unter 
der Voraussetzung, daß sich 
die Austrittsarbeit des Auf- 
fangkäfigs (Graphit) durch 
die Cs-Dämpfe nicht verän- 
dert, wird die Verschiebung 
der Kurven in Richtung ne- 
gativer Spannungen durch 
eine Verkleinerung der Aus- 
trittsarbeit der Cs—Ge- 
Schicht verursacht. 

Sieht man von den Kur- 
ven 2b und 3c in Abb. 7 
zunächst einmal ab, so haben 
die gemessenen Energiever- 
teilungskurven (Abb. 7 u. 8) 
folgendes Ergebnis. Das Kur- 
venmaximum liegt um so 
höher, je mehr der Punkt 
E, = 0in Richtung negativer 
Spannungen U, verschoben 
ist, d.h. also, je mehr sich die 
Austrittsarbeit der Schicht 
verringert hat. Da die maxi- 
malen Verschiebungen der 
Kurven in den Abb. 7 und 8 
je etwa 3 V betragen, muß 
sich die Austrittsarbeit bei 
beiden Schichttypen von 
4,8eV (Germanium)'5) 
auf etwa =1,8eV ver- 
ringert haben. Aus der lang- 
welligen Grenze des äußeren 


50 
25 2b 


Energieverteilung (aut ip normiert) rel Einheiten 


10 

5 

-5 0 5 10 15 2 


Gegenspannung Up/V —= 


Abb. 8. Energieverteilung der Sekundärelektronen von 
Cs—Ge-Schichten auf Grundlage kristallinen Germa- 
niums, auf i, normiert. Kurven 1: Ge (kristallin); 
Kurven 2: Ge + Cs (Formierzustand I); Kurven 3; 
Ge + Cs (Formierzustand II); Kurve 4: Schmiegungs- 


kurve (Erklärung siehe Text). (a und b kennzeichnen IE 


verschiedene Proben) 


Photoeffektes beider Schichttypen?) folgt auch eine Austrittsarbeit von & 1,9eV, 
ein Zahlenwert, der mit = 1,8eV gut übereinstimmt. 


Ferner zeigen die Energieverteilungskurven bei allen mit Cs behandelten 


Ge-Schichten bei Spannungen U, oberhalb des Maximums einen nahezu über- 
einstimmenden Verlauf. Dieses Verhalten ist durch die Hüllkurven 4 in Abb. 7 


») J. A. Dillon, Jr. u. H. E. Farnsworth, J. appl. Phys. 28, 174 (1957). 
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40 Annalen der Physik. 7. Folge. Band 9. 161 
und 8 dargestellt, die aus später ersichtlichen Gründen als quadratische Hyper- f Maxi 
beln gezeichnet sind. Bei Cs—Ge-Schichten auf Grundlage kristallinen Germa- f Schic 
niums hat die Hüllkurve jedoch einen anderen Verlauf als bei den Schichten | runge 
auf Grundlage amorphen Germaniums. Zum Vergleich ist die Kurve 4 aus tinuie 
Abb. 8 in Abb. 7 als Kurve 5 eingezeichnet. Dies wird sicher dadurch verur- f ven ı 
sacht, daß es sich bei den Kurven in Abb. 7 um Schichten mit der Gitterstruktur § sich d 
des CsGe, in Abb. 8 jedoch um Schichten mit Ge-Gitter handelt. Die Energie- | phen 
verteilungskurven für Ge lassen bei Spannungen U, oberhalb des Maximums 
auf Grund ihres Verlaufes eine Zuordnung zu irgendeiner Hüllkurve nur schwer 
erkennen. Aus Abb. 7 geht jedoch deutlich hervor, daß die Kurven 1a—c bei 
Spannungen U,, > 4 V einen anderen Verlauf haben als die Schar der Kurven B 
2—3. In Abb. 8 dagegen laufen die Kurven 1a und b bei U, > 10 V indie Schar | ausge 
der Kurven 2 und 3 hinein. Auch dieses Verhalten scheint darauf hinzudeuten, | 1. die 
daß bei kristallinem Ge nach dem Formieren das Gitter erhalten bleibt, bei | Obert 
amorphem Ge dagegen nicht. der P 

Die Kurven 2b und 3c in Abb. 7 fallen aus dem Rahmen der allgemeinen f den 4 
Ergebnisse etwas heraus. In Abb. 2b weist die negativen Spannungen U, zu- f zelprı 
gekehrte Kurvenflanke eine gewisse Feinstruktur infolge einer erhöhten Ter- § stanz 
tiärelektronenzahl auf. Dadurch liegt der Punkt E, = 0, der durch die Flanken- # zung, 
tangente bestimmt wird, zu weit im Negativen. Bei Kurve 3c hat sich offenbar weise 
auch die Austrittsarbeit gewir 
des Auffangkäfigs etwas Z 
verändert. Hinweise dar- (ig be 
auf, daß sich die Aus- f hoher 
trittsarbeit von Graphit # Einze 
bei sehr langer Ein- 
wirkung von Cs-Dämp- 
fen ändert, finden sich 
bei H. Findeis!®). den I 
In Abb. 9und 10sind | i8¢ di 


Energieverteilung (auf Mox.normiert) 
S 
wn 


20 die Energieverteilungs- stark 

Energie der Sekundärelektronen Es /V ——= kurven in der allgemein Disk 

Abb. 9. Energieverteilung der Sekundärelektronen von üblichen Form darge- D 
Cs—Ge-Schichten auf Grundlage amorphen Germaniums, stellt, nämlich auf Maxi- f Grun 


i iert. i ie in Abb. 7 - 
auf Maximum normiert. (Bezeichnungen wie in ) bzgl. 


des Kontaktpotentials | nach 


~ 


korrigiert. Bei amor-f maß 
pi phem Ge (Abb. 9) liegt | Dafü 
das Maximum bei (sieh: 


2,2eV. Die Halbwerts- ten § 
breite der Kurve beträgt f zufül 


etwa 7eV. CsGe-Schich- 
|  Überschuß zeigen eine} 


S 


Energieverteilung (auf Max .normiert) 
2 
N ™ 


20 recht einheitliche Kur- 


5 10 15 h it . mech 
Energie der Sekunddrelextronen E, ——= 
Abb. 10. Energieverteilung der Sekundärelektronen von 1) H. Findeis, Di = 
Cs—Ge-Schichten auf Grundlage kristallinen Germaniums, plomarbeit, Humboldt-Uni- ” 
auf Maximum normiert. (Bezeichnungen wie in Abb. 8) versität Berlin 1961. = 
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Maximum bei 0,5 — 0,7eV und einer Halbwertsbreite von etwa 1,7 eV. Bei den 
Schichten auf Grundlage kristallinen Germaniums (Abb. 10) sind die Verände- 
rungen nicht so stark. Das Maximum rückt von 1,3eV bei Ge anscheinend kon- 
tinuierlich bis auf 0,5eV herunter, während sich die Halbwertsbreite der Kur- 
ven von 4,5eV bis 2eV verringert. Auch in diesen Darstellungen offenbart 
sich der unterschiedliche Charakter der Cs—Ge-Schichten auf Grundlage amor- 
phen und kristallinen Germaniums. 


6. Diskussion der MeBergebnisse = | 


Bei der Diskussion der MeBergebnisse soll von der bekannten Vorstellung 
ausgegangen werden, daß die SEE in drei Einzelprozesse zu unterteilen ist: 
1. die Anregung der S durch die P, 2. die Diffussion der angeregten S zur 
Oberfläche, 3. der Austritt der S aus dem festen Körper durch Überwindung 
der Potentialschwelle an der Oberfläche. Obwohl der Messung stets nur die in 
den Außenraum austretenden S zugänglich sind, die bereits an allen drei Ein- 
zelprozessen teilgenommen haben, und obwohl ferner die untersuchten Sub- 
stanzen auf Grund der Herstellungsmethode hinsichtlich ihrer Zusammenset- 
zung, ihrer Homogenität und ihres Gitteraufbaues große Unsicherheiten auf- 
weisen, lassen sich doch einige allgemeine Aussagen über die Einzelprozesse 
gewinnen. 


Zunächst haben die Ausbeutemessungen gezeigt, daß es sich bei den mit 
Cs behandelten Ge-Schichten — amorphen und kristallinen — um Stoffe mit 
hoher SEE-Ausbeute handelt. Es steht nun zur Diskussion, welchem der drei 
Einzelprozesse die Erhöhung der Ausbeute hauptsächlich zuzuschreiben ist. 


Die P-Energie E%,,, bei welcher das Ausbeutemaximum auftritt, liefert 
einen Anhaltspunkt über die mittlere Austrittstiefe x, der S und damit über 
den Diffusionsprozeß. Beträgt ES beispielsweise nur einige hundert eV, dann 
ist die mittlere Austrittstiefe gering und die Diffusion der S unterliegt einer 
starken Bremsung durch verschiedene Elementarprozesse. (Eine ausführliche 
Diskussion der möglichen Vorgänge findet sich in!”).) 


Dagegen liegt das Ausbeutemaximum bei höheren Energien E}„, wenn auf 
Grundeinergeringen Bremsung dieSnoch aus größeren Tiefen in den freien Raum 
austreten können. Da sich bei allen Schichten die Lage des Ausbeutemaximums 
nach der Cs-Bahandlung praktisch nicht verändert hat, wird sich auch das Aus- 
maß der Bremsung und damit die Austrittstiefe der S nicht verändert haben. 
Dafür spricht auch die nur schwache Temperaturabhängigkeit der Ausbeute 
(siehe dazu!)). Daraus folgt, daß die Erhöhung der Ausbeute der mit Cs behandel- 
ten Schichten sicher nicht auf eine Änderung des Diffusionsprozesses zurück- 
zuführen ist. 


Die Energieverteilungsmessungen in der in Abb. 7 und 8 dargestellten 
Form geben über die Verhältnisse an der Potentialschwelle Auskunft. Aus der 
klassischen Theorie der SEE von E. M. Baroody'®) wie auch aus der quanten- 
mechanischen Theorie 1?) folgt für die Energieverteilung der im Innern eines 


17) O. Hachenberg u. W. Brauer, Advances in Electronics 11, 413 (1959). 


18) E.M. Baroody, Physic. Rev. 78, 780 (1950). ee 7 | 
19) H. W. Streitwolf, Ann. Physik (7) 8, 183 (1959). tel inde une _ 
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Kristalls angeregten S ein Ausdruck 

(E = Energie der 8, E, = Fermienergie). 


Die in allen Tiefen 0 < x < x, mit dieser Energieverteilung angeregten S 
diffundieren zur Oberfläche, an welcher unmittelbar innerhalb der Potential- F 
schwelle eine Energieverteilung @(E — E,) vorliegt, die der Energieabhängig- f der 
keit der Anregung N;) sehr ähnlich sein wird. Nach Durchlaufen der Potential- 
schwelle verringert sich die Energie sämtlicher in den freien Raum austretenden 
S um den Betrag der Austrittsarbeit der Schicht gy, d. h. die Funktion 
G(E — E, —@) stellt im Bereiche E, = E — E, — y > 0 qualitativ die Energie- 
verteilung der S im Außenraum dar. 


Der gemeinsame Verlauf der Kurven in Abb. 7 und 8 bei Spannungen U, 
oberhalb der Maxima kann nun so gedeutet werden, daß in diesem Bereich 
die MeBkurven dem Verlauf der Energieverteilung @(E — E,) im Innern un- 
mittelbar an der Oberfläche folgen. Man kann sich leicht überlegen, daß dieses 
Verhalten von der Austrittsarbeit der Schicht nicht beeinflußt wird. Insbeson- 
ders lassen sich die Energieverteilungskurven in Abb. 7 und 8 durch die als 
quadratische Hyperbeln gezeichneten Schmiegungskurven (jeweils Kurve 4) 
gut annähern, woraus auch experimentell der ähnliche Verlauf der Energie- 
verteilung G am inneren Rand der Potentialschwelle und der Energieabhängig- 
keit der Anregung £,;) ersichtlich ist. Von dieser Funktion G wird dann beim 
Überwinden der Potentialschwelle gleichsam ein der Austrittsarbeit der Schicht 
entsprechendes Stück abgeschnitten. Abb. 7 und besonders Abb. 8 geben das 
dem beschriebenen Modell entsprechende Verhalten sehr deutlich wieder. Es 
folgt daraus auch, daß es bei einer Verringerung der Austrittsarbeit einer größe- 
ren Zahl von S ermöglicht wird, in den Außenraum zu treten, d. h. die Ausbeute 
steigt an. 


Eine grobe quantitative Abschätzung der Ausbeuteerhöhung allein in- 
folge der Verringerung der Austrittsarbeit führt zu einem ähnlichen Ergebnis 
wie bei Cs—Sb und Cs—Bi!), nämlich, daß die nur durch Einfluß der Verringe- 
rung der Austrittsarbeit verursachte Ausbeuteerhöhung in der Größenordnung 
der bei den Schichten gemessenen Ausbeuteerhöhungen liegt. Dies folgt auch 
schon aus der starken Ähnlichkeit der auf i, normierten Energieverteilungs- 
kurven von Cs—Ge und denen von Cs—Sb und Cs—Bi!). Die Erniedrigung der 
Potentialschwelle ist also auch bei den Cs—Ge-Schichten eine wesentliche 
Ursache für die hohen Ausbeutewerte dieser Schichten. 


Über den Prozeß der Anregung der S geben die Messungen keinen direkten 
Hinweis. Sofern die bisherigen Überlegungen überhaupt eine quantitative 
Abschätzung zulassen, so kann man daraus schließen, daß der Anregungs- 
prozeß wahrscheinlich an den hohen Ausbeutekoeffizienten nur einen geringen 
Anteil hat, da der Diffusionsprozeß keine große Änderung erfährt und die Ver- 
ringerung der Austrittsarbeit schon allein eine den Meßergebnissen ent- 
sprechende Erhöhung der Ausbeute bedingt. 


Die Cs—Ge-Schichten stellen somit Stoffe dar, bei denen Anzeichen dafür 
vorhanden sind, daß die hohen Ausbeutewerte durch eine Erniedrigung der 
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Potentialschwelle an der Oberfläche hervorgerufen werden. Insofern bestehen 
hinsichtlich der SEE große Ähnlichkeiten mit den A!—BY-Verbindungen. 


Herrn Professor Dr. O. Hachenberg danke ich verbindlichst für die An- 
regung zu dieser Arbeit sowie für zahlreiche fördernde Diskussionen. 


Berlin-Adlershof, Heinrich-Hertz-Institut für Schwingungsforschung 
der Deutschen Akademie der Wissenschaften. 


Bei der Redaktion eingegangen am 16. Mai 1961. 
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_ Zur Temperaturabhangigkeit der Sekundärelektronen- 


emission von Ionenkristallen 
Von Bernhard Petzel Ani Pie 


In Ergänzung und Erweiterung einer bereits erschienenen Arbeit!) wurde 
die Sekundäremissionsausbeute 6 an KCl-Einkristallen mit einer verbesserten 
Meßanordnung bei hohen Primärelektronenenergien im Temperaturbereich 
zwischen —180 und +300 °C gemessen. Die beobachtete Temperaturabhän- 
gigkeit von 6 stimmt bei tiefen Temperaturen gut mit dem theoretisch zu er- 
wartenden Verlauf?) überein. Die Aufnahme von Gegenspannungskurven bei 
unterschiedlichen Temperaturen zeigte, daß die Energieverteilung der Sekun- 
därelektronen bei hoher Energie der Primärelektronen geringfügig von der 
Temperatur abhängt. 

1. Einleitung 


Messungen der Sekundärelektronenaubeute 6 an KCl- und KBr-Einkri- 
stallen im Temperaturbereich zwischen Zimmertemperatur und 300 °C!) hatten 
ergeben, daß die Ausbeute 6 bei hohen Energien der Primärelektronen mit 
zunehmender Temperatur etwas steiler abfällt, als es die Dekkersche Theorie?) 
bei alleiniger Wechselwirkung der inneren Sekundärelektronen mit den Gitter- 
schwingungen erwarten läßt. Es erschien daher sinnvoll, die Messungen in das 
Gebiet tiefer Temperaturen auszudehnen, um hieraus gegebenenfalls weitere 
Aufschlüsse über den Diffusionsprozeß der angeregten inneren Sekundärelek- 
tronen zu gewinnen. Messungen der Sekundärelektronenemission an Ionen- 
kristallen bei tiefen Temperaturen liegen bislang nur an aufgedampften dün- 
nen Schichten vor’). 

Neben den Ausbeutemessungen sollte durch Messung der Energieverteilung 
der Sekundärelektronen bei unterschiedlichen Temperaturen die aus der Dek- 
kerschen Theorie folgende Aussage nachgeprüft werden, daß die Energiever- 
teilung von der Temperatur nicht abhängt. 


2. Meßanordnung 


Die Meßanordnung war gegenüber der in der früheren Arbeit!) verwendeten 
verbessert worden. Um die Aufladungen auf der Kristallloberfläche so klein 
wie möglich zu halten, wurde als Vorverstärker ein rauscharmer ballistischer 
Impulsverstarker *) ®) mit einer Ladungsempfindlichkeit von 2 -10-% As bei 
einem Signal- Rausch-Verhältnis 40:1 verwendet. 


B. Petzel, Ann. Physik (7 ). 6, 55 (1960). 

A. J. Dekker, Physica 21, 29 (1955). 

M. W. Gomojunowa, Sh. techn. Fis. 26, 2399 (1958). 

N. R. Whetten u. A. B. Laponsky, Physic. Rev. 107, 1521 (1957). 
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Petzel, in Vorbereitung. 
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Die Meßröhre®) für die Ausbeutemessungen bei tiefen Temperaturen besaß 
zwei Kühlfinger. Der eine von beiden diente als Probenträger. Er bestand aus 
einem Kupfertopf, der über ein Kovarzwischenstück mit dem Glasmantel der 
Meßröhre verschmolzen war. Zum Kühlen der Probe wurde flüssige Luft 
in den Kupfertopf eingefüllt, während zur Erwärmung ein kleiner Heizkörper 
eingeführt werden konnte. Der zweite Kühlfinger war seitlich angebracht und 

os ermöglichte eine Kühlung des Kollektors. Er wurde stets zuerst mit flüssiger 
Luft gefüllt. Die Temperatur auf der Probenoberfläche wurde mit einem 
Kupfer-Konstantan-Thermoelement gemessen. 

Die Meßröhre für die Energieverteilungsmessungen besaß keine Kühl- 
möglichkeit, die Messungen blieben auf Temperaturen oberhalb der Zimmer- 
temperatur beschränkt. Innerhalb des Kollektors befand sich eine kugelförmige 
Netzelektrode zur Unterdrückung der Tertiärelektronen. Zu diesem Zweck 

irde $  erhielt der Kollektor ein gegen die Netzelektrode um 10 V erhöhtes Potential. 
ten | Der gegenseitige Abstand der Maschendrähte des Netzes war klein gegen den 
ich Abstand zwischen Netz und Kollektor. Hierdurch sollte die Verminderung 
an- | des Energie-Auflösungsvermögens der Gegenfeldanordnung infolge des Durch- 
a griffs durch die Netzmaschen klein gehalten werden. 
bei Beide Meßröhren besaßen eine lösbare, hochvakuumdichte Verbindung mit 
are einem Dichtungsring aus 0,1 mm dicker Kupferfolie, der zwischen zwei kraf- 
der tigen, aus rostfreiem Stahl bestehenden Flanschen mit entsprechend bearbeite- 
ter Oberfläche zusammengepreßt wurde’). Jede Meßröhre wurde vor dem 
i Einsetzen der Probe evakuiert und 30 Stunden bei 420 °C ausgeheizt. Durch 
kri- Lösen der Flanschverbindung wurde die Meßröhre geöffnet und ein frisch 
‚ten gespaltener KCl-Einkristall eingesetzt. Die Messungen erfolgten an laufender 
mit Pumpe bei Drucken um oder unter 10-7 Torr. Vor Beginn der Messungen wurde 
ie?) jeder Kristall auf etwa 300 °C erwärmt, ein weiteres Ausheizen der Meßröhre 
wurde nicht vorgenommen. 
lek- Die Messungen erfolgten 
mit Einzelimpulsen in der [ 
ün- bereits beschriebenen 
Weise!). Für die Ausbeute- 7 
messungen wurden imTarget- | 
ek- kreis mit dem ballistischen 
Verstärker an Stelle der Strö- st j 
me die entsprechenden La- 4t 
dungen gemessen. Die Wieder- 
lein potentials zwischen den Mes- 80 90700 30 
pe‘ leitung. Dieses Verfahren ge- Ab». 1. Ausbeute ö in Abhängigkeit von der absoluten 
bei Temperatur 7 für zwei Primärelektronenenergien, ge 
währleistet,daßkeineinneren messen an einem KCl-Einkristall. Gestrichelte Kurve 
Raumladungen unterhalbder berechnet nach der Formel 
Oberfläche zurückbleiben. 6 ~ {2[exp(hv/k T)— 17 + 

6) Herrn M. Müller und Herrn A. Werner möchte ich auch an dieser Stelle fiir die 
sorgfältige Ausführung der glastechnischen und feinmechanischen Arbeiten danken. _ 

7) R. Brymmer u. W. Steckelmacher, J. Sci. Instrum. 36, 278 (1959). iti 


- | 
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Die Pod wiesen bei Verwendung einer Meßröhre mit Fiaudthivedbinidlune ausg 
höhere Ausbeutewerte auf als bei Benutzung einer allseitig verschmolzenen ände 
Meßröhre. Im Maximum der Ausbeutekurve ergaben sich an KCl-Einkristallen } vers 
Werte bis ömax ~ 25. Die charakteristische Temperaturabhängigkeit der Aus- f pera 
beute bei hohen Primärelektronenenergien war jedoch nicht verändert. Abb. 1 
enthält die an einem KCl-Einkristall gemessenen Ausbeutewerte in Abhän- | App 
gigkeit von der absoluten Temperatur für zwei Primärelektronenenergien. f Sek 


Ausbeute und Temperatur sind Beo 
logarithmisch aufgetragen. Die Ten 
Energie der Primärelektronen f der 
u, x Ts 35% ner liegt mit 6 bzw. 8 keV genügend Ene 

5 106 weit oberhalb der zum Ausbeute- 
E Ep = 250ev ff 104 maximum gehörenden Energie, so ] 
daß die Eindringtiefe der Primär- f Arb 
102 elektronen groß im Vergleich zur bei 
er MR 08 -6 4 2 0 +2 tronen ist. In die Abb. 1 ist der der 1 
Spannung zwischen Probe u. Netzelektrode(V) auf Grund der Theorie von Dek- | oat, 

zer?) si Terl: 

Abb. % An bei einer 7 dee 
Primärelektronenenergie von 250 eV gemessene Ausbeute eingetragen. Besonders : 
Gegenspannungskurven. Die Kurven wurden bei Temperaturen unterhalb der WR 


seitlich verschoben zwecks Korrektur des Kon- Zimmertemperatur ist eine gute 


ss taktpotentials Übereinstimmung zwischen den 
gemessenen und berechneten Wer- 
| : pr ten festzustellen. Bei höheren 
Arey | Temperaturen ist die Überein- 
4 kh we aes } = stimmung nicht so gut. Es wurde 
106 bereits vermutet!), daß in diesem 
} Ep = 6 kev Bereich die Grundvoraussetzung 
dr 19* der Dekkerschen Theorie, daß 
die inneren Sekundärelektronen 
I nur mit den Gitterschwingungen 
wechselwirken, nicht mehr voll 
Spannung zwischen Probe u.Netzelektrode(V) erfüllt ist. 


In Abb. 2 und 3 sind je zwei 
Abb. 3. An einem KCl-Einkristall bei einer hei 35 °C und 200 °C aufgenom- 
Primärelektronenenergie von 6 keV gemessene mene Gegenspannungskurven für 
Gegenspannungskurven. Die Kurven wurden 


seitlich verschoben zwecks Korrektur des Kon- (ie Primärelektronenenergien 250 
taktpotentials und 6000eV dargestellt. q, ist 


die von Netz wor Kollektor bei 
der betreffenden Gegenspannung aufgenommene Ladung, qkx die entspre- 
chende Ladung bei einigen Volt Absaugspannung. Bei einer Energie von 250 eV 
ist die Eindringtiefe der Primärelektronen beträchtlich kleiner als die Reich- 
weite der Sekundärelektronen und demzufolge die Ausbeute von der Tempe- 
ratur unabhängig!). Auch die Energieverteilung der Sekundärelektronen wird 
in diesem Fall durch die Temperatur erwartungsgemäß nicht beeinflußt (Abb. 2). 
Bei hoher Primärelektronenenergie sollte nach Dekker?) die Energieverteilung 
gleichfalls von der Temperatur nicht abhängen. Alle Proben, an denen von uns 
die entsprechenden Messungen durchgeführt wurden, zeigten eine schwach 
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ausgeprägte, aber eindeutig außerhalb der Meßungenauigkeit liegende Ver- 
änderung der Gegenspannungskurve bei Änderung der Temperatur. Nach Abb. 3 
verschiebt sich das Maximum der Energieverteilung bei Zunahme der Tem- 
peratur geringfügig in Richtung höherer Energiewerte®). 

Wir glauben nicht, daß die in der Dekkerschen Theorie verwendeten 
Approximationen (unendlich ausgedehntes Medium, räumlich gleichförmige 
Sekundärelektronenproduktion) für die Abweichung zwischen Theorie und 
Beobachtung verantwortlich sind. Wir halten es für möglich, daß das für höhere 
Temperaturen vermutete Auftreten anderer Wechselwirkungsprozesse neben 
der Bremsung der Sekundärelektronen durch die Phononen eine Änderung der 
Energieverteilung verursachen kann. 


Herrn Prof. Dr. Hachenberg danke ich für die Anregung zu dieser 
Arbeit, Herrn J. Griebel für seine Hilfe beim Aufbau der Meßapparatur und 
bei der Durchführung der Messungen. 


8) Schulman u. Dementjew [Sh. techn. Fis. 25, 2256 (1955)] haben bei Messungen 
der Energieverteilung an einem NaCl-Einkristall bei 1200 eV Primärenergie keine Tempe- 
raturabhängigkeit der Energieverteilung festgestellt. 


Berlin-Adlershof, Heinrich-Hertz-Institut für Schwingungsforschung 
der Deutschen Akademie der Wissenschaften. 


ik 
Bei der Redaktion eingegangen am 16. Mai 1961. ar 
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Röntgen-K- und L,,,-Niveaubreiten 


Von Armin Meisel und Wadim Nefedow  — 


7 Inhaltsübersicht 


2 Unter Benutzung von experimentellen Werten fiir die K-Niv eaubreiten 

oe und Ka,-Halbwertsbreiten ergab sich, daß man die Abhängigkeit der K- und 

__ [qyy-Niveaubreiten von der Ordnungszahl Z in der Form (Z + P)* gut wieder- 
geben kann. Vergleiche mit dem Experiment und mit anderen Berechnungen 
von Niveaubreiten wurden angestellt. 


~~ 


Nach quantentheoretischen Vorstellungen setzen sich die Breiten der 
Röntgenemissionslinien additiv aus den Breiten der an den betreffenden 
Elektronenübergängen beteiligten Energieniveaus zusammen. Bei Annahme 
einer Multiplettstruktur der Linien!) sollte jedoch der Beitrag eines bestimmten 
Niveaus zur Breite verschiedener Linien nicht konstant sein. Man kann z. B. 
von einem Beitrag des Zisr-Niveaus zur Ka,-Linienbreite sprechen, wenn man 
sich die Ka,-Halbwertsbreite linear aus den K- und Lin-Niveaubreiten zu- 
sammengesetzt denkt; der Beitrag von Ly zur La,-Breite wird sich aber 
im allgemeinen davon unterscheiden, sofern man wieder die La,-Linienbreite 
als Summe der Zyr- und My-Niveaubreiten betrachtet. Dies geht schon dar- 
aus hervor, daß das Lyy-Niveau beim Entstehen der La,-Linie ionisiert ist, 
was beim Entstehen der Ka,-Linie nicht der Fall ist. Wendet man dagegen an- 
dere Additionsschemen als die lineare Zusammensetzung an, so kann der Bei- 
trag eines Niveaus zu den Halbwertsbreiten verschiedener Linien durchaus 
gleich sein. In diesem Sinne wollen wir im folgenden versuchen, die Beiträge 
der K- und Zur-Niveaus zu den Ka,-Linienbreiten zu bestimmen. 

Von der Multiplettstruktur müssen auch die Absorptionsspektren beein- 
flußt werden. Was man gewöhnlich als kontinuierliche Absorption bezeichnet, 
sollte eine endliche Anzahl von unaufgelösten bzw. nicht auflösbaren einzelnen 
Linien sein. Ferner stellen die gemessenen Niveaubreiten (Tab. 1) ein Maß für 
die Beiträge der entsprechenden Niveaus zu den Linienbreiten dar. Der Grad 
der Zuverlässigkeit der in dieser Tabelle angeführten Werte ist jedoch recht 
unterschiedlich. Bereits Blochin und Ssatschenko?) erwähnten, daß die 
von Beeman und Friedman?) angegebenen Niveaubreiten für die Über- 
gangselemente nicht richtig sein können; da die unaufgefüllte 3d-Schale an 


1) 8. folgende Arbeit. 

2) M.A. Blochin u. W. P. Ssatschenko, Nachr. Akad. Wiss. UdSSR, phys. Ser. 
21, 1343 (1957). 

3) W. W. Beeman u. H. Friedman, Physic. Rev. 56, 392 (1939), 
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der Absorption teilnimmt, ist die Anwendung der Methode von Richtmyer, 
Barnes und Ramberg’) nicht gerechtfertigt. Bei den Daten für Mg und Al®) 
handelt es sich um korrigierte Werte, bei denen die Korrekturen etwa 75% 
der gemessenen Größen (1,5 bzw. 1,25 eV) betragen. Die K-Niveaubreite für 
Ag kam dadurch zustande, daß man ,,somewhat ; 
arbitrarily‘‘ den Wert fiir Zırı zu 2eV schätzte und Tabelle 1 3 
dann mit Hilfe der Ka,-Linienbreitenmessungen von Experimentell bestimmte _ 


Allison®) die K-Niveaubreite bestimmte”). Der _ X-Niveaubreiten in eV, 
nach der Zusammenstellung 


Wert für Au stellt auch keinen Meßwert dar; rela- yon Blochin und Ssa- 
tivistische Rechnungen?) lieferten für die Au K- tschenko?) 
Niveaubreite 47 eV und konnten außerdem die DM 05 
Lß,, ß,-Intensitätsverhältnisse richtig wiedergeben. 13 AL 0,35 
Man sieht, daB das experimentelle Material sehr 18 Ar 0,58 
unvollständig und unsicher ist. Strenggenommen 19K 0,73 
lagen keine experimentellen Beweise fiir ein Z4- = = LB f 
Gesetz der K-Niveaubreiten vor. — Von Zinn?) 28 Ni | oo 
wurden für die K-Niveaubreiten Werte veröffent- 29 Cu au 
licht, die um mehr als das Dreifache die gemessenen 30Zn a 
Ka,, a,-Halbwertsbreiten übertreffen. Da die dort = ome 1,90 
angeführte Methode völlig unbegründet ist, kann 47 Ag 8 i 
man den Z4-Verlauf der MeBwerte nicht als einen 74W 37 
experimentellen Beweis ansehen. — Die Verwendung 79 Au 62 5 
gemessener Fluoreszenzausbeuten bildet zwar eine ” 
gute Stütze für die Schlußfolgerung, daß für Z> 30 ein Z4-Gesetz gelten 
sollte?); sie leidet aber daran, daß zwei zweifelhafte Meßwerte für die K- N 
Niveaubreiten (Ag und Au) benutzt wurden und nur drei annehmbare Wette —_— 
(Ga, Kr und W) übrig blieben. ee. 
Die in der vorliegenden Arbeit angewandte Methode zur Bestimmung der 
K- und Zir-Niveaubreiten besteht in folgendem: 
1. Weil für die Ka,-Halbwertsbreiten einAEx., = A (Z + P)*-Gesetz gilt”), 
kann man ein ähnliches Gesetz auch für die K- und L[yy-Niveaubeitrage ‘a 
den Ka,-Linienbreiten erwarten (in Analogie zu den Moseleyschen Gesetzen 
für die Energie der ausgestrahlten Photonen und die Energiewerte der einzelnen 
Terme). 
2. Für ein (Z + £)*-Gesetz spricht auch die Gleichung für die Fluoreszenz- i ee 
ausbeuten Wx nach Burhop!!): 
( Wx a4+bZ+ mit a = — 0,044 a) 
(für alle Z b= 0,0346 ° 
c = — 1,35 - 10% 
4) F.K. Richtmyer, S. W. Barnes u. E. Ramberg, Physic. Rev. 46, 843 (1934), 
5) J. H. Munier, J. A. Bearden u. C. H. Shaw, Physic. Rev. 58, 537 (1940). a 
6) S.K. Allison, Physic. Rev. 44, 63 (1933). Ps; 
?) L. G. Parratt, Physic. Rev. 54, 99 (1938). oa 
H. W. S. Massey u. E.H.S.Burhop, Proc. Cambridge philos. Soc. 82, 454 
JOD). 
*) W. H. Zinn, Physic. Rev. 46, 659 (1934). Ku 


10) A. Meisel u. W. Nefedow, Z. Chem. 1 (1961), H. 11 (im Druck). 
u) E.H.S.Burhop, J. Physique Radium 16, 625 (1955). N 
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Bei Vernachlässigung des dritten Gliedes erhält man daraus: 


Wy \1/4 = b4(Z + d)* 
b(Z d) oder Wx N . 
Da andererseits die Fluoreszenzausbeute definiert ist als d 


(/", = Strahlungsbeitrag, I’, = Auger-Beitrag, J’, = gesamte Niveaubreite), 
sollte für die K-Niveaubreite gelten: 

AE, 1+ + dy. (2) 
Strenggenommen müßte man statt der /'-Werte die entsprechenden Über- 
gangswahrscheinlichkeiten setzen; bei Annahme einer Multiplettstruktur der 
Linien müßten weiterhin die entsprechenden Energiebreiten für jede einzelne 
Multiplettkomponente getrennt berechnet werden. Man kann jedoch die ge- 
messenen W „-Werte als über alle Multipletts gemittelt betrachten, wie es oft 
bei Messungen der L-Fluoreszenzausbeuten (ohne L;-, Lyy-, 
geschieht. 

3. Es bestehen zwei voneinander unabhängige Möglichkeiten zur Über- 
prüfung der Richtigkeit unserer Ansätze: 

a) Bei gegebener Form der Gleichung kann man die Interpolation zwischen 
zwei Punkten (in diesem Falle Kr und W) als einfache Festlegung betrachten, 
die K-Niveaubreiten bestimmen und damit die Ziır-Beiträge zu den Ka,-Linien- 
breiten aus den Ka,-Halbwertsbreiten für 36 < Z < 56 berechnen. Diese Bei- 
träge sollten auch die gleiche Form der Z-Abhängigkeit haben: 

AE,,,, = + MM. 

b) Durch Extrapolation dieser Werte kann man die Ly;-Niveaubreiten 
fiir solche Elemente bestimmen, fiir die experimentelle Daten vorliegen, und 
die Ubereinstimmung priifen. 

Auf Grund der unter 1 und 2 angefiihrten Punkte setzen wir fiir die K- 
Niveaubreiten (dabei soll A natürlich nicht als irgendeine Abschirmungskon- 
stante betrachtet werden): 

= g(Z + h). 
Unter Benutzung der fiir Gleichung (1) angegebenen Zahlenwerte ergibt sich 
aus (2) 
(AEx)'/* ~ 0,0346 (Z — 1,27). 

Bei Beriieksichtigung der beiden experimentellen Werte von Kr und W erhält 
man eine sehr ähnliche Gleichung: 

(ABz)K = 0,0355 (Z— 1,2) für Z> 36. (3) 
Wie schon aus Dimensionsgründen hervorgeht, ist die zahlenmäßige Überein- 
stimmung dabei vollkommen zufällig. Gleichung (3) bietet die Möglichkeit, 
verschiedene AEx-Werte miteinander zu vergleichen. Absolute Werte lassen 
sich auf diese Weise wegen fehlender genauer experimenteller Daten nur an- 
nähernd berechnen. 

Nach (3) wurden für die Elemente 37 < Z < 50 die K-Niveaubreiten und 
unter Benutzung der Ka,-Linienbreiten von Gokhale!) die Zyr-Niveau- 


2) B.G.Gokhale, Ann. Physique 7, 852 (1952), 
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breiten ermittelt. Man kann auf diese Weise die erhaltenen Zyrr-Werte "2 
folgender Beziehung gut wiedergeben (Abb. 1): eats 
= 7,15 - (Z + 128,5) für Z> 36. (4) 


Tab. 2 zeigt die befriedigende Übereinstimmung der von Blochin und | 
Ssatschenko?) mit den von uns berechneten K-Niveaubreiten; die Unter- _ 
schiede in den Zirr-Niveaubreiten sind aus- 
schließlich dadurch bedingt, daß für ihre 
Bestimmung von Blochin und Ssa- 
tschenko linear korrigierte Ka,-Halbwerts- 
breiten verschiedener Autoren benutzt wur- 


(aE, 


ist 


den. Es kommt aber gar nicht so sehr auf die 72} m 3 = 
absoluten Werte an, da die Genauigkeit der (DE, (1285) 
K-Niveauwerte ohnehin klein ist, als viel- Pf 
mehr auf den allgemeinen Gang, und dieser 77 f f ei 

ist in der genannten Arbeit schlecht ausge- 36 40 “4 48 


prägt. Die Unterschiede der Werte in den Abb.1. L,,;-Niveaubreiten in eV 2 
Spalten 4 und 5 liegen im Meßfehlerbereich als lineare Funktion von (Z + P)* 
der Untersuchung Gok hales?%). 


Tabelle 2 
K- und L,,,-Niveaubreiten in eV 
AE, | AE 
Element nach (3) in ?) aus nach (4) in ?) 


berechnet berechnet berechnet berechnet 


37 Rb 2,46 1,94 1,96 
38 Sr 2,13 201 | 1,68 
39 Y 2,06 2,06 
40 Zr 3,40 | 3.56 2,12 2,11 1,85 
41 Nb 3,77 | | 22 2,16 2090 
42 Mo 4,16 410 | 2,28 2,21 1,76 a 
43 Te 4,58 | | 2,26 in 
44 Ru 5,03 | 2,15 2,31 1,94 = 
45 Rh 5,52 | 5,38 | 2,36 2,36 1,92 
46 Pd 6,05 | | 239 2,42 1,99 
47 Ag 6,60 6,28 | 2,54 2,47 2,32 
48 Cd 7,20 | | 2,69 2,53 
49 In 7,83 | | 2,60 2,59 
50 Sn 2,65 

In Tab. 3 sind die nach (4) Tabelle 3 


Berechnete und gemessene 


bereehneten Zim-Niveaubreiten in eV 


mit experimentellen Werten ver- 
glichen ; die Übereinstimmung ist Element | „nach (4) | experimentell 


qualitativ recht befriedigend. Da | 
die Genauigkeit. mit der sich K- 


und Lyy-Beitrage zur Ka,-Halb- | +05 
wertsbreite bestimmen lassen, 78 Pt 3.6 + 
nieht allzu groß ist, kann man  —79 Au 4,4 
diese Werte auch als Beiträge der = 82 Pb 5.1 
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nn genannten Niveaus zur Breite anderer Linien betrachten; sie können ferner zur 
Orientierung über den absoluten Wert sowie zur Wiedergabe der Abhängigkeit 
von Z dienen. 
Das gleiche Verfahren wurde auf die Elemente mit 10 <Z< 36 ange. 
_ wandt. Zur Bestimmung der Z-Abhängigkeit der K-Niveaubreiten dienten 
_ die in Tab. 1 angeführten Werte für Mg, Al, Ar, K, Zn und Ga; sie ergaben die 
Gleichung 
vg (AEx)i¢ = 0,0166 (Z + 37,1) für Z< 36. (5) 
die Ermittlung der AE,, „-Formel benutzten wir unsere Messungen der 


as Halbwertsbreiten von Cu, Zn, Ga, Ge und Se?) — die entsprechenden 
Werte der Übergangselemente kann man aus verschiedenen Gründen nicht 


dazu verwenden — und erhielten 


(AE — 0,0248(Z + 11,6) für 8<Z< 36. (6) 


Die Genauigkeit der Gleichung (5) kann man aus Tab. 4 abschätzen. Eine 
Möglichkeit der indirekten Überprüfung von (6) enthält die Arbeit der Verfasser 
über die Multiplettstruktur der Röntgenemissionslinien!). Von den dort an- 

geführten drei verschiede- 


Tabelle 4 nen Wegen zur Berech- 

q Vergleich von K-Niveaubreiten in eV nung der „ungestörten" 

Element | ‚nach (5) in ?) experimentell 

| berechnet berechnet | (Tab. 1) jergangselemente be 

E | | steht einer in der Anwen- 

j 12 Mg | 0,44 0,44 0,5 dung der Gleichung (6): 

13 Al | 0,48 0,48 0,35 die gute Übereinstimmung 

u 30 Zn | 1.54 1,40 15 spricht fiir die Brauch- 

> 31 Ga 1,63 1,53 1,9 barkeit der hier abgelei- 
Bus 36 Kr 2,17 2,37 2,07 teten Beziehung. 

13) D. Coster u. H. de Lang, Physica 3, 282 (1936). 

14) J. = Bearden u. T. M. Snyder, Physic. Rev. 59, 162 (1941). 


15) R. E. Shrader, Physic. het. 49, 644 (1936). 


Leipzig, Physikalisch-Chemisches Institut der Karl-Marx-Universität. 
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Kia aes Bei der Redaktion eingegangen am 19. Mai 1961. 
| 


A thie wh joe 


Beit 


D 
pletts 
- (AE), 
AE 
bedet 
Linie 
an de 
breite 
häng 
„ung 
die n 
Z 
Rönt 
B 
renz 
den 
trons 
same 
nach 
hen, 
Diffe 
besti 
Wert 
dene 
Diffe 
hera: 
+4 
gren: 
Unsi 
und 


zur 


gkeit 


n die sii zur Multiplettstruktur von Réntgenemissionslinien 
Armin Meisel und Wadim Nefedow 


Die Arbeit enthält einige empirische und theoretische Belege für die Multi- 
(6)§ plettstruktur der Röntgenspektren. U. a. wurde gezeigt, daß die Gleichung 


Eine} (42) mn = = AE,, + AE,, wo (AELZ),,,, die Linienbreite (Halbwertsbreite) und 
„, AE,, die "Breiten der am Elektronenübergang beteiligten Energieniveaus 
ao bedeuten, für Röntgenlinien nicht gilt; daß die Cr Kx,-Linie von Cr,O, aufge- 


2 ede.) spalten ist und daß eine weitgehende Analogie zwischen Dublettabständen und 
| Linienbreiten besteht. Für das Verhältnis von Dublettaufspaltung eines der 


roch an der Entstehung der betreffenden Linie beteiligten Niveaus zur Halbwerts- 

adil breite dieser Linie ergab sich über einen großen Z-Bereich eine lineare Ab- 
hd hingigkeit von der Ordnungszahl. Mit Hilfe einer solchen Beziehung konnten 

swell „ungestörte‘‘ Kx,-Halbwertsbreiten der Übergangselemente berechnet werden, 

| (6) die mit auf zwei anderen Wegen erhaltenen Werten gut übereinstimmen. 

mung 

auch- Zunächst wird gezeigt, daß die Annahme einer Multiplettstruktur der 


gelei.) Röntgenspektren keiner der bisher bekannten Tatsachen widerspricht. 
Bekanntlich kann man die Energie eines ausgestrahlten Photons als Diffe- 
| renz zweier Terme berechnen. Betrachtet man zwei Emissionslinien, die durch 


den Ubergang eines Elek- 


trons von einem gemein- Tabelle 1 
samen Niveau auf zwei be- Differenzen der Au-Niveaubreiten AE aa AE ie ineV, 
nachbarte Niveaus entste- nach Richtmyer, Barnes und Ramberg?) 


tät. 


hen, so ist es möglich, die 
Differenz zweier Terme zu | 
bestimmen. Ein solcher 05 | 36 075 | 1,1 
Wert läßt sich auf verschie- 
denen Wegen berechnen: So kann man z. B. zur Bestimmung der Zyr—Zumn- 
Differenz die Linienpaare K (a, — 04), L(y und L(y; 
heranziehen. Man bekommt eine recht gute Übereinstimmung bis auf 
+ 4 eV 2)3), was — relativ gesehen — klein ist und oft innerhalb der Fehler- 
grenze liegt. Im Vergleich mit den absoluten Halbwertsbreiten ist jedoch diese 
Unsicherheit sehr beachtlich und entspricht etwa dem Unterschied der Zır- 
und Ian-Niveaubreiten, die man aus den Linienbreiten erhält (vgl. Tab. 1). 


Li—n | — By Be—rYs | 


0,2 


| > 


K.] S. W. Barnes u. E. Ramberg, Physic. Rev. 46, 843 
*) M. Siegbahn, Spektroskopie der Röntgenstrahlen, Berlin 1931. 
3) E. Ingelstam, Nova Acta Regiae Soc. Sci. Upsaliensis (TV) 10, Nr. 5 (1937). 


n der 
»nden 
nicht 
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Die weitere Aufspaltung der Terme kann auch so stattfinden, daB der 
Schwerpunkt des Multipletts unbeeinflußt bleibt; das trifft z. B. für die Auf. 
spaltung im Magnetfeld zu. Nach Sommerfeld‘) hängt die durchschnittliche 
Energie einer Elektronenkonfiguration nicht von der Spin-Bahn-Wechsel. 
wirkung ab, wenn jedes Niveau mit dem entsprechenden statistischen Gewicht 
2J +1 (J = Gesamtdrehimpulsquantenzahl) genommen wird. Dagegen 
können die Auswahlregeln unter Umständen den Schwerpunkt eines Multi. 


l 
Knie flee 2 
Abb. 1. Aufspaltung der 2 p®- bzw. 2 p!-Konfiguration a) im 
starken Magnetfeld, b) ohne Magnetfeld 
J 
2 
3 TER 
wh ann T 
J=0000111 1-1-1441 
LB, iM; 
Ka, 170041 
?L 
73 
| ‘4 
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i Abb. 2. Mögliches Schema zur Multiplettstruktur 
Zusammenwirkung der Elektronen innerhalb einer Schale mit gleicher innerer 
Quantenzahl stattfindet. Bei Beriicksichtigung des Zusammenwirkens von 


Elektronen mit solchen in unvollständig besetzten Schalen müßte eine Auf- 
spaltung stattfinden. Im Falle der J—J-Bindung (schwere Atome) kann man 


4) A. Semmertold, Atombau und Spektrallinien, Braunschweig 1951. N fe 
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von 4 Elektronen auf dem Ly;-Niveau und 2 auf Zy sprechen. Allen Erwartun- 
gen nach sollten also die Röntgenlinien aus vielen dicht nebeneinander liegen- 
den Multiplettlinien bestehen. Z. B. sollte sich Ka, aus 2 Linien zusammen- 
setzen, weil im starken Magnetfeld die p®-Konfiguration in 3 Niveaus aufge- 
spalten wird, von denen sich 2 in einem schwachen Magnetfeld zusammen- _ 
ziehen 5) (Abb. 1). Die Aufhebung der Entartung kann auf verschiedenen Wegen 
erfolgt sein. 

Setzt man — wie die Gründe auch dafür sein mögen — vollständige Auf- __ 
spaltung voraus, so sollte man bei Berücksichtigung der Auswahlregeln den _ 
einzelnen Linien die aus Abb. 2 ersichtlichen Strukturen zuschreiben können. 
Das aufgestellte Schema läßt sich auch für andere Linien in analoger Weise 
vervollständigen. Ganz allgemein kann man erwarten, daß die Übergänge mit 
1J = 0 mit kleinerer Wahrscheinlichkeit erfolgen als die mit AJ = + 1. Bei 
der Deutung der Intensitätsverhältnisse der Linien läßt sich möglicherweise 
hieran anknüpfen. 

In der Literatur findet man keine sicheren Angaben über die Aufspaltung 
von Linien. Vegard®) wies zum ersten Male auf die richtigen Grundlagen zur 
Berechnung der Halbwertsbreiten von Röntgenemissionslinien hin; danach 
muß sich die beobachtete Linienbreite aus den natürlichen Linienbreiten der 
einzelnen Multipletts und ihrem Abstand zusammensetzen. Schmudski’) 
ging von den gleichen Voraussetzungen aus und schlug folgende halbempirische ° 
Formel zur Bestimmung der Linienbreiten vor: re 


+6)2 2 5 wily fale 
(nm p+ 1)? 4,50, +1 halam 
nis 27 Mn, + + 5g—10) 


ch ’ tg + + 4q—10 

P, und P, geben die Anzahl der gekoppelten Elektronenpaare oder der unge- 
paarten Elektronen der Anfangs- und Endniveaus wieder; n, ! und J bedeuten 
die Quantenzahlen; g und gq werden nicht erklärt. Die Formel gibt nur den 
allgemeinen Verlauf der Z-Abhängigkeit der Linienbreiten in grober Annähe- 
rung wieder. Die auf Grund dieser Beziehung vorgenommene Berechnung 
von Energieniveaus ist deshalb zweifelhaft. In der genannten Arbeit wurde 
auch von einer Aufspaltung der Zn Kx,-Linie berichtet, was aber — nach 
unserer Kenntnis — bisher weder bestätigt noch widerlegt wurde. 

Die von Krasnikow®) beobachtete Aufspaltung der Cr Ka,,,-Linien von 
Cr,0, in 2 Komponenten der gleichen Intensität erklärte Wainstein®) durch 
die Spin-Wechselwirkung der Zy- und Zyrr-Elektronen mit den Valenzelek- 
tronen; danach sollte man beim Cr,O, nur dann eine Aufspaltung erwarten, 
wenn in dieser Verbindung zwei verschiedene Valenzzustände des Cr vorliegen. 
Setzt man jedoch nicht alle Zin- bzw. Ly-Elektronen energetisch gleich, so 


läßt sich die Aufspaltung ohne zusätzliche Annahmen erklären. — Die Auf- 
5) M. A. Blochin, Physik der Röntgenstrahlen, Berlin 1957. \ 
L. Vegard, Z. Physik 21, 6 (1920). 
?) A.S. Schmudski, Lonitomasch 15, 273 (1949). VEN Vy 


8) A. J. Krasnikow, Ber. Akad. Wiss. UdSSR 49, 346 (1945). 
’) E. Je. Wainstein, Die Röntgenspektren der Atome in chemischen Verbindungen 
und Legierungen, Moskau 1950. 
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_ spaltung selbst können wir bestätigen (Abb. 3), und auch die von Skinner, 
Bullen und Johnston™) angegebene Form der Cr Zysr-Bande von Cr,0, 
(Abb. 4) zeigt gewisse Aufspaltungserscheinungen. Ebenso kann man nach den 


genannten Autoren bei einigen anderen Oxy. 
den eine Aufspaltung der Zim-Bande er- 
warten. 

Die Hauptfrage besteht darin, ob die 
Aufspaltung eine nennenswerte Bedeutung 
hat. Man soll also feststellen, ob die ein- 
zelnen Niveaubreiten oder die Abstände 
zwischen ihnen den Hauptanteil zu den 
Linienbreiten im üblichen Sinne liefern. Es 
ist auch sehr wahrscheinlich, daß beide Tat- 
sachen einen meßbaren Einfluß ausüben. Die 
Verhältnisse können ferner verschieden sein 
bei schweren und leichten Elementen sowie 
bei unterschiedlichen Niveaus, ähnlich wie 
es bei der Dublettaufspaltung der Fall ist. 
Wenn die Niveaubreiten einen großen Anteil 
an den Linienbreiten haben, sollte man auch 
einen meßbaren Einfluß von Auger-Uber. 
gängen erwarten. Der Zusammenhang zwi- 
schen Linienbreiten und Auger-Übergängen 
wurde von Cooper”) diskutiert. Da die an 
den Halbwertsbreiten angebrachten Korrek- 
turen nicht gerechtfertigt sind, sollte man 
jedoch mit den unkorrigierten Werten rech- 
nen; unter diesen Bedingungen lassen sich 
aber die in der genannten Arbeit gezogenen 
Schlußfolgerungen nicht bestätigen. 

Im folgenden sind die empirischen Tat- 
sachen, die die Multiplettstruktur der Rönt- 
genemissionslinien wahrscheinlich machen, 
; 

1. Man kann zeigen, daß es unmöglich 


; ist, jedem Energieniveau eine solche Breite 
7. “3 IE, zuzuschreiben, daß sich alle Linienbreiten linear zusammensetzen lassen; 


für einfache Linien war dies theoretisch zu erwarten2). 


- (AE), = Halbwertsbreite der Linie, die dem Übergang E,, 
Bekanntlich wird die Niveaubreite durch 


IE 
a 


10) H. W. B. Skinner, T. G. Bullen u. J. E. Johnston, Philos. Mag. [7] 45, 1070 
Physic. Rev. 61, 234 (1942). 


(1954). 
u) J.N. Cooper, 


V. WeiBkopf u. E. Wigner, 


2. 68, 54 (1930. 


off ai 
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wiedergegeben, wobei w,, die Wahrscheinlichkeit des Überganges En > E, 
darstellt. Damit haben wir für eine Linie1*) 
2 7 p q 

Eine theoretische Berechnung für das K-Niveau von Au ergab AE = 67,15 eV, 
wobei die Strahlungsbreite AE, = 66,38 und die Auger-Breite AZ, = 0,772eV 
betragen!#). Die experimentelle Linienbreite für das Au K-Niveau soll nach 
Richtmyer und Barnes") 


sicher unter 60eV liegen; die Tabelle 2 
theoretische Berechnung Berechnete und experimentell bestimmte Linien- 
kann also höchstens die rich- breiten in eV 


tige Größenordnung wieder- ner berechnet! experimentell 
geben. Es besteht jedoch für gang nach!) | nach!$) | nach'!”) 
einzelne Niveaubreiten eine AuMx, M,N,,, 9 
„experimentelle“ Bestim- 


mungsmöglichkeit, die theo- 6 M 9,6 
retisch von Richtmyer, w 
1 
Barnes und Ramberg’) 98 8 
23, 20 


begriindet wurde. Danach 
ermittelt man _,,experimen- 
tell‘‘ die Zir- und Zim-Niveaubreiten, und zwar letztere mit besonderer Ge- 
nauigkeit. Die Gültigkeit der Gleichung (1) vorausgesetzt, läßt sich nun aus 
der Kx,-Linienbreite die K-Niveaubreite berechnen. Für die Bestimmung von 
Lu kann man und heranziehen, weil die Diffe- 
renz der Breiten dieser Linien gleich der 


Differenz von Ly; und Ly; ist. Mit Hilfe Tabelle 3 
der Lx,-Linienbreite läßt sich auch Myy Ag-Niveaubreiten in eV, nach Parratt'®) 
berechnen. Bei Benutzung verschiede- | berechnet | experimentell 
ner L-Linien kann man weitere Niveau- - 
breiten bestimmen. L, | 5,3 10? 

Man soll viele Unsicherheiten dieser Ly | 2,2 6? 
Methode beachten: So ergeben z. B. die i | 2 (2) 


Differenzen für Lim und Zir, aus Ab- 
sorptionskanten experimentell be- 
stimmt, — 0,7 eV, aus verschiedenen Linienbreiten ermittelt, jedoch positive 
Werte (Tab. 1). Abgesehen von den Differenzen, zu deren Bestimmung schwer 
vermeßbare Linienbreiten verwendet wurden, ist die Schwankung der ange- 
führten Werte trotzdem sehr hoch. Man kann nun Gleichung (1) mit Hilfe 
der von Richtmyer, Barnes und Ramberg!) berechneten Au-Niveau- 
breiten überprüfen (Tab. 2). Eine ähnliche Bestimmung wurde von Parratt!®) 
für Ag durchgeführt. Er setzte ganz willkürlich 2 eV für die Zirr-Niveaubreite 
und berechnete daraus die übrigen Niveaubreiten. Wie groß die Differenzen 
zwischen den experimentell gefundenen und berechneten Werten sind, zeigt 
die Tab. 3. 

13) A. Bril, Physica 13, 481 (1947). * 

14) E.G. Ramberg u. F. K. Richtmyer, Physic. Rev. 51, 913 (1937), 

15) F.K. Richtmyer u. 8. W. Barnes, Physic. Rev. 48, 1049 (1933). 


16) K. Molin, Dissertation Uppsala 1927. een A 
17) L. G. Parratt, Physic. Rev. 50, 598 (1936). 


18) L. G. Parratt, Physic. Rev. 54, 99 (1938). W 
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Tab. 4 gibt einen Vergleich einiger von Kiessig!?) experimentell bestimm- c) 


ter Ag-Linienbreiten mit den entsprechenden Werten, die mit Hilfe der von f fache 
Parratt}8) angegebenen Niveaubreiten berechnet wurden. Wie man sieht, ist | spaltı 
die Übereinstimmung von Theorie und Experiment unbefriedigend*). Halb 


Eine auf Grund der Gleichung (1) durchgeführte Berechnung der Fe Zun- ren 
I Niveaubreite lieferte eine negative Größe); auch für die anderen Ubergangs- re 
elemente ergaben sich negative 


Tabelle 4 Sat Werte. Hier aber darf man nicht oh 
Ag-Linienbreiten in eV Gleichung (1) anwenden, weil die f . 
Üb — Linien der Übergangselemente } 
ergang | berechnet | experimentell’*) geno! 
eine groBe Asymmetrie aufweisen st 
08-13 | 3,8 und von einer Dispersionsform 
My 1.6 | 3,7 auch nicht annähernd gesprochen 
24 15 werden kann. 
| 2. Es besteht eine weitgehende 
MN, 3,0 8,6 Analogie von Spindublettaufspal- 
tung und Linienbreite: 
a) Der Abstand der Spindubletts bleibt im Längenmaß etwa konstant, 
ändert sich aber im Energiemaß sehr stark mit zunehmender Ordnungszahl. 
Das gleiche — wenn auch nicht in solehem Maße — gilt für die Linienbreiten. 
So vergrößert sich z. B.der Ka,, &-Spindublettabstand von 22 Ti bis 74 W im 
Längenmaß um das 1,3fache, im Energiemaß aber um das 224fache; für die | yo / 
Linienbreiten findet man im Längenmaß eine Verkleinerung um das 6fache | „nd 
und im Energiemaß um das 33,5fache. Bei großem Z ist diese Erscheinung die \ 
noch stärker ausgeprägt. Es sei auch fund 
NEM gly) bemerkt, daB die Breite der Multipletts A 
40 im Längemaß bei leichten Elementen 2 
m sehr viel zur Spindublettaufspaltung | . d 
beiträgt. Die Verkleinerung der Län- 
‘pt genmaßverhältnisse bei den Multiplett- 
40 Pr abständen wird deshalb wesentlich klei- 8 
ner sein, alses bei den Halbwertsbreiten 
2 der Fall ist. Mit 
b) Bekanntlich ist der regelmäßige tört 
6 n & 8 2 H W Z Verlauf der Ka, &%-Spindublettab- 
_ Abb. 5. Ka-Dublettabstände in X stände bei den Ubergangselementen n - 
3 gestört (1) undAbb. 5). Das gleiche trifft om 
für die Ka,, x.-Halbwertsbreiten zu #4), wobei auch das Maximum der Störung ed 
bei beiden Erscheinungen übereinstimmt. Die Änderung des Ka,, «,-Spin- Mö | 
dublettabstandes läßt sich nicht nur auf eine Änderung von Form und Breite a 
der K&,, «,-Linien zurückführen. Man kann vielmehr behaupten, daß Spin- ter“ 
dublettabstände sowie Form und Breite der Linien unabhängig voneinander u 
beeinflußt werden und daher von ähnlicher Herkunft sind. liefe, 
*) Zur Möglichkeit der Wiedergabe der Halbwertsbreite einer Linie durch die Breiten es 
der beteiligten Energieniveaus siehe die vorangehende Arbeit der Verfasser. ” 
1%) H. Kiessig, Z. Physik 95, 555 (1935). “ 
20) L. Obert u. J. A. Bearden, Physic. Rev. 54, 1000 (1938). 
ss #1) A. Meisel u. W. Nefedow, Z. Chem. (1961), H. 11 (im Druck). 1, 5 


< 


4 


‘von | facher empirischer Zusammenhang. Für das Verhältnis der Kx,,&,-Dublettauf- 


t, ist | spaltung zu den Ka,- bzw. Ka,- Via 
Halbwertsbreiten im Energiemaß 
L erhielten wir eine lineare Abhän- yin ahh 
a gigkeit von der Ordnungszahl ? | 
atid (Abb. 6). Man könnte ohne weiteres 
nicht | gesamten Bereich 200<Z<74 Mr 
die | eine einzige Gerade ziehen. Es wurde *Q62(2-157) 
endl jedoch eine Teilung bei Z = 36 vor- st 
genommen, damit eine bessere Uber- 
form | instimmung auch mit genaueren RH 6 2 


chen %-Halbwertsbreiten für 37 <Z< 50 standes zur Linienbreite von Ka, (©) bzw. 


spal entnommen; unsere eigenen Messungen?!) lieferten die Kx,, .-Linienbreiten 
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c) Zwischen Spindublettaufspaltung und Halbwertsbreite besteht ein ein- 


Meßwerten möglich wird. Die Ka, Abb. 6. Verhältnis des K«a-Spindublettab- 
wurden der Arbeit Gokhales?) Ka, (@) 


für 22 <Z< 34. Nach der Methode der linearen Ausgleichung erhält man: 
AE (Lyy— Lin) 


— 0,62 (215,7) für <Z< 36 (2) 
AE,,, 

4 nm _ 9.515 (Z — 10,4) fir Z> 36, (3) 
(9) 


wo AE (Ly — Lyrr) die Energiedifferenz zwischen den Z,;- und Ly;;-Termen 
und AF x,, die Ka,- Halbwertsbreite bedeuten. Die Gleichung (3) liefert fiir 


die W Ka,-Linienbreite 41 eV; experimentell wurden 43%) und 42 eV*4) ge- 
funden. 


Aus Abb. 6 erkennt man, daß die Linienbreiten der Ubergangselemente 


vom regelmäßigen Verlauf ziemlich stark abweichen. Diese Abweichungen 
sind bei Ka, größer als bei Ka, und erreichen ein Maximum bei Z = 27—28; 
das Maximum für Kx, liegt bei Z = 25-26. Man kann nun aus Abb. 5 die 


„ungestörten‘“ Werte 
des Spindublettabstan- Tabelle 5 
des erhalten (Tab. 5). Ka-Dublettabstande 
Mit Hilfe dieser „unge- AE (Ly, — „ungestört“ Lim) 
stérten‘‘ Spindublettab- Element | un experim., 
stände und Beziehung (2) | _inX | inky in Ry 
lassen sich ‚„‚ungestörte‘ 20 Ca | 3,28 0,26 | 0,26 ; 
Kx,-Halbwertsbreiten 21 Sc 3,36 0,33 033 a 
berechnen. Zwei andere 28 Ti 3,42 0,41 eo 
Möglichkeit rB 23 \ | 3,48 0,51 0,55 i 
Möglichkeiten zur Be- 24Cr 3.54 0,62 0,68 
stimmung „ungestör- 25 Mn 3,59 0,74 0,82 

ter“ Ka,-Linienbreiten 26 Fe 3,64 | 0,88 0,96 

der Ubergangselemente 37 Co 1m 

dis 28 Ni 3,75 | 1,24 1,28 

h die ichungen 29 Cu 3,80 1,46 1,47 
(I 3) sowie (II 5) mit 30 Zn 3,86 1,71 1,71 


22) B. G. Gokhale, Ann. Physique 7, 852 (1952). 

3) F.K. Richtmyer u. S. W. Barnes, Physic. Rev. 46, 352 (1934). 
*) B. B. Watson, W. J. West, D. A. Lind u. J. W.M. DuMond, Physic. Rev. 
505 (1948). 
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(11 6)*). Alle diese Wege sind voneinander unabhängig mit der Einschrän- 
kung, daß in allen drei Fällen für die Extrapolation die Kx,-Halbwerts. 
breiten von Cu, Zn, Ga, Ge und Se benutzt wurden. Tab. 6 enthält außer 
den auf verschiedene Weise berechneten ‚ungestörten‘‘ Kx,-Halbwertsbreiten 
die K- und Z,rı-Niveaubreiten in dem Sinne, wie es in unserer vorangehenden 
Arbeit über die Niveaubreiten erläutert wurde. Die gute Übereinstimmung 
der angeführten Werte bestätigt zugleich die Brauchbarkeit von Gleichung 
(II 6). 
Tabelle 6 
Auf verschiedene Weise berechnete „ungestörte‘“‘ Kx,-Halbwertsbreiten sowie K- und 
L,ır-Niveaubreiten in eV 


„ungestört‘‘ AE, AE, 
Z | | | mach nach (IT 5) 116 
5 | nach (2) nach (I: | und (II 6) nach ( ) 
21 | 1,36 1,28 1,30 0,87 0,43 
22 | 1,43 1,40 1,41 0,93 0,48 
23 1,51 1,52 1,53 | 0,99 0,54 
tig 1,74 1,13 0,68 
At 26 1,87 1,96 | 1,96 | 1,20 ; 0,76 
27 2,04 2120 2,12 | 1,28 0,84 
# 28 2,22 2300 2,29 1,36 * 0,93 
29 2,40 248 | 2,48 1,45 1,03 


{ Die lineare Abhängigkeit des Verhältnisses von Spindublettabstand zu 
Halbwertsbreite von Z läßt sich auch bei anderen Linien nachweisen: Als 


JE (MgMg) Ergebnis der graphischen Mittelung der von 

Ingelstam?) für 42<Z< 78 angegebenen 
u Wr 04% KB,-Halbwertsbreiten (Abb. 7) erhielten wir 
40 25% die Gleichung 


15 

1E(M,,— M 


202 Beim Vergleich von (2), (3) und (4) fallt auf, 
45 535 60 1 75 daß die Koeffizienten von Z einander sehr 
Abb. 7. Verhältnis des Kß,,,-Spin- ähnlich sind. Man muß jedoch beachten, daß 
dublettabstandes zur Linienbreite im Falle (4) die Messungen eine andere Ge- 
von Kp, 
nauigkeit besitzen. 

Eine Auswertung der von Williams?) erhaltenen Ergebnisse liefert im 

Bereich 74 =< Z < 92 die beiden Gleichungen: 


*) (13) bedeutet Gleichung (3) aus der Arbeit der Verfasser über Linienbreiten®! 
(IL 5) Gleichung (5) aus der Arbeit über Niveaubreiten. ZH 
5) J.H. Williams, Physic. Rev. 37, 1431 (1931); 44, 146 (1933). (Hel; OMe 
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Aus obigen Beispielen sollte man erwarten, daß Gleichung (6) die richtige 
ist. Die Extrapolation auf Z = 47 ergibt nach (5) AEagr., = 2,68, nach (6) 
2,25eV; nach Davis und Purks?) beträgt der nicht korrigierte Wert 2,72. 
während der korrigierte bei 2,34 eV liegt. Da die Messungen von Williams?) 


nicht korrigiert sind, ziehen wir den unkor- 
rigierten Wert vor und entscheiden uns somit La) 
fir die Gleichung (5), die in Abb. 8 gra- tay ag 
phisch dargestellt ist. 20 
Man kann auch mit Hilfe der unkorri- af ER 
gierten Werte von Cooper!!) und der Mes- 220 FE He 25) 
sungen von Bril!2) zeigen, daß 
200 
AE 
Da die Genauigkeit der Lf,-Linienbreiten 6 & 6 WO 


jedoch nicht groß ist, lassen sich keine ge- Abb. 8. Verhältnis des Dyy yyy-Du- 
nauen Angaben über Aund B machen. Der blettabstandes zur Linienbreite von 


Wert fiir A entspricht etwa dem aus (2), (3) La, 
und (4), der Wert für B etwa dem aus (5). u 

Für die oben angeführten Gleichungen ist noch keine theoretische Er- 7 
klärung vorhanden. Es ist bemerkenswert, daß man nach (2), (3), (4) und (7) 7 
Differenzen von Termen bilden muB, die eigentlich nicht viel zu den Linien- 7 
breiten beitragen sollten. Es handelt sich jedoch in diesen Fällen um Linien, q 
die ein ,,Singulett‘‘-Niveau im üblichen Sinne als Endzustand haben, nämlich 7 
K für Kx,, 6, und L, für LB,. Wegen der Unvollständigkeit der erhaltenen Er- f 
gebnisse wird von einer Verallgemeinerung und einer weiteren Besprechung u 
der Gesetzmäßigkeiten abgesehen. u 


3. Die Abhängigkeit der Form und Breite der Linien von den Anregungs- 
bedingungen?”) läßt sich zwar auch ohne Annahme der Multiplettstruktur h 
deuten; sie zeigt jedoch, daß die Niveaubreiten sehr empfindlich gegenüber | 
der Anregungsart der Röntgenlinien sind, was die Multiplettstruktur wahr- 7 
scheinlich macht. 


2) B. Davis u. H. Purks, Physic. Rev. 34, 181 (1929). id 7 
7) A. Meisel u. W. Nefedow, Exp. Techn. Physik (im Druck). = 


m ‘Bei der Redaktion eingegangen am 19. Mai 1961. a u 
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er 2) Helmut Scheffers, Ann. Physik 18, 29 (1956). 


| Kar Determination of the Mean Free Path of 


Conduction Electrons) 


: By R. Kleinberg and F.J. Blatt oneh 


obra el y 
Abstract 


An earlier paper by Scheffers [Ann. Physik 33, 29 (1956)] is criticized. 
It is shown that Scheffers result for the magnetoresistance ratio of a metal 
is based on incorrect boundary conditions, and that the agreement with experi- 
ment in the case of the noble metals is fortuitous. + How idan 

wry 


Scheffers?) suggested a method for calculating the number of free elec- 
trons, N, and the electronic mean free path, J, in a metal which has an assumed 
spherical energy surface. He started with the classical equation of motion?) 

mr+m yr=elE+2 x H| 
E = (E,. 0, 0), H = (0,0, H,), J, = Nev; =c E, 


he: 
and showed that 


9 
(2) 
From this he obtained 


with 


Combining this result with the well known expression 


one obtains 
SAT 1 H (00/02 0 0) xd (5) 
> e CO, V Aeley oF, 
-5(7) N28, (6) 


bs a Using thd Kohler dingueme for copper, silver, and gold, Scheffers 
deduced numerical values for N and o,l, which are in very good agreement 
with those obtained from anomalous skin effect measurements. 


1) Supported in part by the Office of Scientific Research, US Air Force, and the 
Office of Ordnance Research, US Army. 


3) The notation is that adopted by Scheffers. ey. 
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It is now firmly established that the Fermi surface in copper touches the 
Brillouin Zone at the hexagonal zone faces‘). Thus, it deviates considerably 
from the spherical shape assumed by Scheffers. It is most probable that the 


Fermi surfaces in silver and gold resemble closely that in copper 5). Moreover, it is 
also well known that standard free electron theory predicts a vanishing magne- 
toresistance under the assumption of a spherical energy surface®). Indeed, a 
theoretical account of the observed magnetoresistance of copper demands a 
pronounced deviation from sphericity’). It is, therefore, curious that Scheffers 
should have arrived at a non-vanishing result for 4o/oq - 0, and most surprising 
that his result is in good agreement with experimental observation. 


We wish to point out that the non-vanishing result obtained by Scheffers 
is a direct consequence of his choice of incorrect boundary conditions. The 
numerical agreement is thus entirely fortuitous. 

In the normal experimental arrangement for measurement of magneto- 


resistance the boundary conditions imposed are ER Ku = a 
J = (Jz, 0,0) H = (0, 0, H,). 


This situation leads to a Hall field £, in addition to the externally imposed 
field E,, i.e. E= (Ez, E,,0). The resistivity is defined by = 


where @ is the Hall angle. In the range of applicability of the Boltzmann 
transport treatment @ is small and proportional to H. 


In contrast, Scheffers imposes the boundary condition E = (E,,0, 0). 
It follows that, for H, + 0, J must have a transverse component J, as well as 


longitudinal component J, 
- 


wird 

of) — ose? (9) 

In this way, even though Ao/oy - » = 0, Scheffers obtains a nonvanishing 

magnetoresistance ratio 


Thus, we have 


wah H? ) H? (10) 


where R 1/Nec is the Hall constant. 


4) A.B. Pippard, Rep. Progress in Physics, Vol. XXIII, p.176 (The Physical 
Society, London 1960). 

5) P. G. Klemens, Handbuch der Physik, Vol. XIV, p. 242. (Springer Verlag, 1956). 

6) A. H. Wilson, The Theory of Metals (Cambridge University Press, 1953). 

) M. G. Priestley, Philos. Mag. 5, 111 (1960). 

8) „Dabei ist nämlich angenommen, daß Strom und konstante Feldstarke E die 
t-Richtung haben und das Magnetfeld H = H, senkrecht dazu wirkt.“ This “assumption” 
presupposes a physical impossibility. Near the end of his paper Scheffers reproduces the 
standard treatment of the Hall effect and here this assumption is disregarded. 
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Of all the monovalent metals sodium probably has a Fermi surface most 
nearly spherical. Using Scheffers expression, Eq. (5), and the experimental 
data®) we obtain values for N ranging between 10 and 20 - 1022 cem®. This is 
to be compared with N = 2.3 - 1022 cm? found by Babiskin and Sieben- 
mann?®) and N = 2.5- 10% cm? from the Hall coefficient®). Evidently the 
agreement becomes worse just where Scheffers should expect it to be best. 
The rather good agreement in the cases of copper, silver and gold is clearly 
fortuitous. 


®) D.K.C. MacDonald, Philos. Mag. 2, 97 (1957); Proc. physic. Soc. 68 A, 290 (1950). 
4) Babiskin and Siebenmann, Physic. Rev. 107, 1249 (1957). 


East Lansing, Michigan (USA), Physics Department, Michigan State 
University. 


ee Bei der Redaktion eingegangen am 10. April 1961. 
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early Über eine Näherungslösung 
der Einsteinschen Feldgleichungen, _ 

welche ein quellenfreies darstellt 

Inhaltsübersicht 


Es wird in zweiter Näherung eine Lösung der Einsteinschen Feldglei- 
chungen R,, = 0 angegeben, die folgende Eigenschaften hat: Ryyx, = 0: 
Singularitätenfreiheit; g“” — n"” für r— oo; positive Gesamtenergie. Die 
physikalische Bedeutung einer derartigen Lösung wird diskutiert. as 2 


Die physikalische Bedeutung des Problems 


Die vorliegenden Untersuchungen beschäftigen sich mit der Möglichkeit, 
ein quellenfreies Gravitationsfeld zu konstruieren. Um die physikalische Be- 
deutung derartiger Betrachtungen einzusehen, orientieren wir uns zunächst 
an der speziell-relativistischen Feldtheorie. 


Hier erscheinen die einzelnen Felder als Beschreibung selbständig existie- 
render Materiearten; nach dem Vorbilde der phänomenologischen Materie 
haben sie demgemäß ihre eigene Energie, ihren eigenen Impuls, Drehimpuls 
usw., treten sie miteinander in Wechselwirkung, sind sie miteinander gekoppelt. 
Eine eindrucksvolle Rechtfertigung dieser Auffassung ergibt die quantisierte 
Form der Feldtheorien, welche die jeweils zu einem Felde gehörigen, vonein- 
ander verschiedenen und beobachtbaren Teilchen liefert. Freilich wird durch 
die Möglichkeit der Teilchenumwandlungen die Selbständigkeit der quanti- 
sierten Felder relativiert. Eine solche Relativierung tritt aber auch schon in 
der klassischen Formulierung insofern auf, als es im Falle wechselwirkender 
Felder im allgemeinen nicht möglich ist, z. B. die Gesamtenergie auf die 
Partner ohne Willkür zu verteilen. 


Aus diesem Grunde wird die Konzeption von Feldern, die relativ selbständig 
existierende Materiearten beschreiben, ihre logisch stärkste Stütze in der 
Möglichkeit finden, isolierte freie Felder zu konstruieren, zumindest unter 
idealisierten Verhältnissen. Die Idealisationen, um die es sich hier handelt, 
werden verschieden zu beurteilen sein, wie wir am Beispiel eines Feldes mit 
elektrischer Ladung und eines Maxwellschen Feldes erläutern wollen. Gehen 
wirnämlich von dem Satz aus, daß jedes elektrisch geladene System von einem 
Maxwellschen Feld begleitet ist, so wird es im Prinzip überhaupt kein freies 
elektrisch geladenes Feld, z. B. ein Diracsches Feld geben. Die Annahme eines 
solchen hat daher prinzipiell den Charakter einer Näherung, die sich auf die 


5 Ann. Physik, 7. Folge, Bd. 9 


Ki 


_ Hiermit ist gesagt, daß die Lösung der Einsteinschen Feldgleichungen über- 


PER os begleitenden Maxwellschen Feldes bezieht. Hingegen schließt 


der oben erwähnte Satz keineswegs die Möglichkeit eines freien elektromagne. 
tischen Feldes prinzipiell aus. Unter diesem Gesichtspunkt hat die Konstruk.§ weit e 
tion physikalisch sinnvoller quellenfreier Maxwellscher Felder eine besonderf muß « 
Bedeutung für die oben geschilderte Konzeption, die der klassischen speziell. überh 
3 relativistischen Feldtheorie zugrunde liegt. D. 
Es erhebt sich nun die Frage, ob man diese Konzeption auch auf die allge. Es 
meine Relativitätstheorie, auf das Gravitationsfeld übertragen kann. Beschreibt welch 
das Gravitationsfeld in dem dargelegten Sinne eine beseuidere, selbständiger und I 
Existenz fähige Materieart mit eigener Energie und ähnlichen Qualitäten,f den E 
„Gravitationsmaterie‘ ? Oder ist das Gravitationsfeld nur eine Eigenschaftfi weise 
aller Materie? Die entsprechende quantentheoretische Frage lautet: Gibt zu ve 
es Gravitonen, die sich eventuell in andere Elementarteilchen umwandehf zieler 
können ? Di 
. Wir beschränken uns auf den klassischen Aspekt. Dabei stellen wir hier durch 
auch die Besonderheit des Gravitationsfeldes, für die Metrik unseres Raum.§ sein 1 
_ Zeitkontinuums bestimmend zu sein, nicht zur Diskussion. Eine andere Beson- 
derheit des Gravitationsfeldes ist seine Universalität, insofern jedes ander 
_ materielle System Quelle eines Gravitationsfeldes ist. Diese Besonderheitf Es we 
können wir beurteilen, indem wir die Masse (genauer den Energie-Impuls.# Die € 
tensor) irgendeines Feldes als die Quelle eines Gravitationsfeldes in Parallele ein k 
setzen mit der Ladung (genauer dem Viererstrom) irgendeines elektrisch ge-# über 
= = _ladenen Feldes als Quelle eines Maxwellschen Feldes. Die bezüglich der letzt-# Diese 
_ genannten Fakten angestellten Überlegungen übertragen sich dann auf die$ Grin 
_ jetzt interessierenden Verhältnisse: Die Konzeption eines jeden freien, speziell A 
-relativistischen Feldes hat nun einen — durch die Schwäche der Gravitations-# os kei 


_ effekte gerechtfertigten — Näherungscharakter, während die Möglichkeit einesff entwi 
physikalisch sinnvollen, quellenfreien Gravitationsfeldes noch nicht prinzipiellf gibt e 
_ ausgeschlossen ist. Die Frage nach seiner Existenz dürfte für das Problem wiez 
der physikalischen Natur des Gravitationsfeldes, für das Problem der ‚‚Gravi-# diese: 
tationsmaterie‘‘ von grundsätzlicher Bedeutung sein. Nähe 
Wir geben nun vier Eigenschaften an, die bei Zugrundelegung der Ein-§ wird 
steinschen Theorie (ohne kosmologisches Glied) ein derartiges Feld haben? füllt, 


müßte. In dis 
A. Ruvxa = 0. 


haupt ein echtes Gravitationsfeld beschreibt. 
B. Abwesenheit felderzeugender Quellen. > 

Singularitäten der Gravitationspotentiale müßten nämlich vom physikalischen 8) 

Standpunkt her als Quellen des Feldes interpretiert werden. » 

C. Erfüllung der Grenzbedingung u 
Juv > Nur fiir r> oo de la 


bei Verwendung pseudokartesischer Koordinaten (7,, ist der metrische Tensor f und 1 
des Minkowskiraumes in diesen Koordinaten). Diese Grenzbedingung wird | all sc! 
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man zu wählen haben, wenn man sinnvoll von isolierten physikalischen 
Systemen sprechen kann, d. h. von solchen, deren Wechselwirkung mit anderen 
weit entfernten Systemen vernachlässigt werden darf. Eine solche Möglichkeit 
muß offenbar bestehen, wenn wir den Begriff quellenfreier Gravitationsfelder 
überhaupt sinnvoll benutzen können. 

D. Positive Gesamtenergie. 

Es gibt zwar eine Anzahl Lösungen der Einsteinschen Feldgleichungen, 
welche den Bedingungen A und B geniigen")?))*)5), sie erfüllen aber nicht C 
und D. In Anbetracht der Schwierigkeiten, die der Aufstellung einer exakten, 
den Bedingungen A bis D genügenden Lösung entgegenstehen, bzw. dem Be- 
weise der Nichtexistenz einer solchen Lösung, ist es gerechtfertigt, zunächst 
zu versuchen, mit Hilfe geeigneter Näherungsverfahren Teilergebnisse zu er- 
zielen. 


Die Entscheidung über das zu verwendende Näherungsverfahren wird 
durch die Bemerkung erleichtert, daß wegen C das Feld bei großen r schwach 
sein muß und somit dort die folgende Entwicklung gestatten soll: 

ger = + + fees, (1.1) 
Es werden hierbei die kontravarianten Dichten g“” als Feldvariable verwendet. 
Die Grenzbedingung C erfordert das Verschwinden aller ı9“” fiir r — oo. e ist 
ein kleiner Entwicklungsparameter. Ferner wird angenommen, daß das Feld 
überallschwach sei, d.h. daß die Entwicklung (1.1) im ganzen Raum gelte. 
Diese Annahme wird nicht aus physikalischen, sondern aus mathematischen 
Gründen gemacht. 


Auf dieser Basis gelang es Papapetrou®)?)8)9)10) vor allem zu zeigen, daß 
es keine Gravitationsfelder gibt, welche A bis D erfüllen und deren Fourier- 
entwicklung bezüglich der Zeit ¢ einen diskreten Anteil haben. Insbesondere 
gibt es keine zeitlich periodischen Felder dieser Art. Der Satz von Lichnero- 
wiez!!), nach dem keine derartigen stationären Lösungen existieren, ist in 
dieser Feststellung als Spezialfall eingeschlossen!*). Obwohl mit Hilfe eines 
Näherungsverfahrens hergeleitet, handelt es sich um exakte Ergebnisse. Es 
wird nämlich zunächst bewiesen, daß ,g“” die Grenzbedingung nur dann er- 
füllt, wenn der zeitlich periodische Anteil von ,g“” wegtransformierbar ist. 
In diesem Fall kann aber der Beweisgang für ,q”” und ,g“” wiederholt werden, 
und schließlich auf die Eliminierbarkeit des gesamten zeitlich periodischen An- 
teils von g“” geschlossen werden. Wie man sieht, besteht die Idee des Beweis- 


1) A. Einstein u. N. Rosen, J. Franklin Inst. 228, 43 (1937). 

*) N. Rosen, Bull. Res. Coun. Israel 3, 328 (1954). 

3) H. Bondi, Nature 179, 1072 (1957). 

4) W. B. Bonnor, J. Math. Mech. 6, 203 (1957). 

5) J. Weber u. J. A. Wheeler, Rev. Mod. Physics 29, 509 (1957). 

6) A. Papapetrou, Ann. Physik 20, 399 (1957). 

7) A. Papapetrou, Ann. Physik 1, 186 (1958). 

8) A. Papapetrou, Ann. Physik 2, 87 (1958). 

°) A. Papapetrou u. H. Treder, Ann. Physik 8, 360 (1959). 

10) A. Papapetrou u. H. Treder, Ann. Physik 6, 311 (1960). 

u) A. Lichnerowicz, C. R. Acad. Sci. Paris 222, 432 (1946); Théories relativistes 
de la Gravitation et de l’Electromagnetisme (Masson Cie, Paris 1955). 

12) Die erwähnten Theoreme konnten in den zitierten Arbeiten von Papapetrou 
und Treder auch für wichtige Spezialfälle bewiesen werden, in denen das Feld nicht über- 
all schwach ist. 
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ganges in der Aufdeckung grundsatzlicher Schwierigkeiten bei der Berechnung 
der höheren Näherung, wenn eine annehmbare niedere Näherung voraus 
gesetzt wird. 


Es war nun naheliegend, diese Idee auch für den Fall eines Feldes, das den 


_ Bedingungen A bis D gehorchen soll und bezüglich ¢ ein rein kontinuierliche 
_ Spektrum hat, auszuprobieren. Da die Untersuchungen Papapetrous hier. 


über keinen unmittelbaren Aufschluß geben, wurde, um eine Orientierung zı 


gewinnen, ein spezielles singularitätenfreies, die Grenzbedingung erfüllende al * 
ı9“” angenommen, welches keinen diskreten Anteil im Fourierspektrum be. Pr” 
züglich ¢ hat. Sodann wurde untersucht, ob sich ein ebensolches ,g“” ergibt. 
Hierüber berichten die folgenden Darlegungen, deren Ergebnis hier vorwegfi Wie d 
genommen sei: drati 

Die Antwort ist positiv. Es wird gezeigt, daß sich in der zweiten Näherung gi 
. ein Gravitationsfeld ergibt, welches alle Bedingungen A bis D erfüllt. Selbst. rg 
2 verständlich läßt dieses Ergebnis die Frage offen, ob es nicht durch höher 2 4) 
m Näherungen wieder zerstört wird. Außerdem dürfte das berechnete Feld in et di 
j der Natur nicht realisiert sein. Die Bedeutung unseres Ergebnisses wolle man 
. vielmehr daran ermessen, daß die oben erwähnte Frage nach der Natur de 
 Gravitationsfeldes in den meisten Fällen auf der Basis der ersten Näherun 
diskutiert wurde!®), und daß — wie gesagt — ein ‚g“” mit diskretem Fourier- 
 anteil bezüglich ¢ bereits zu grundsätzlichen Schwierigkeiten für ,g“” führt 

schrei 

4 prepare § 2. Entwicklung der Rechenmethode Hi 

Wir benutzen in unseren Rechnungen die de Dondersche Kon 
einfül 
= 9 1) 
= 0, Ferne 
die für jede Stufe des Näherungsverfahrens erfüllt sein soll: kan y 
v v ay ean 
19” ee = (2.2) Man « 
: 3 Hiermit erhält man für ‚g“” die d’Alembertsche Gleichung er 
Oa" 02” Mit 
. Ihre allgemeine, unseren Anforderungen genügende Lösung kann in der Form 
7 e') Pi (cos dw Ar”) (w) = Fi do ,y"" (24) 

a pals mor „were 

. a 13) Siehe etwa D. Iwanenkowu. A. Sokolow, Klassische Feldtheorie, Akademie-Ver- mit d 
lag, Berlin 1953; $ 56. 7 
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wobei I,,,; die Besselfunktion der Ordnung n + 1/2 ist. Die Realität von 
‚9 wird durch einige Forderungen an die Gewichtsfunktionen A“’} in (2.4) 
gesichert, während durch weitere Beziehungen zwischen diesen Größen die 
Koordinatenbeziehung bereits in der Form ; 
== 0 (2.7) x 

erfüllt werden kann. Die Lichtgeschwindigkeit wurde gleich 1 gesetzt. 

Mit (2.1), (2.2) erhält man für ,g“” eine inhomogene d’Alembertsche 
Gleichung 

D = (2.8) | 
Wie diese Schreibweise andeuten soll, ist die Inhomogenität 3“” eine qua- 
dratische Form in den igen der ersten Näherung. Ihrekon- © 
krete Gestalt findet man 1. c.*) Formel (3.10), (3.13).14). 

Da nun die Ableitungen von Wellenfunktionen wiederum Wellenfunktionen 
sind, haben die Ableitungen von ,g”” dieselbe Gestalt wie die rechte Seite von 
(2.4). Es werden sich lediglich die ändern. 
wird man für %#” 


dey dQ eito+Dt Cx Q) Ye (wr) (2 r) 


Hier kann man die Entwicklung => 


1aten- Pi PL =a, +42. +: ‘+a, +.. 
(2.1) einführen. ” ist in derselben Weise a oder ungeradzahlig, wie k + K. Cae 5 
Ferner gilt die Beziehung 


n>K-—kfürk<K. 
(2) Man erhält sodann mit neuen Gewichtsfunktionen Dy’,™ 


+00 +00 


mn 


schreibt sich 


(wr) pR(Qr) dw dQ. 

Mit 

Form y,(or) y„(Qr) = r) 


Ger = Leime Pm ff da dQ = [ [dw dQ jer. (2.14) 
(2,5) | Unserer von (2.8) liegt nun folgender Gedanke zugrunde: Wir 
machen den Ansatz Er 
+00 +00 2 
o.Ve- Die dort auftretenden Größen v“” bzw. w”” sind identisch mit ,g”” und.g” nicht 
mit den in dieser Arbeit so bezeichneten Größen zu verwechseln. 7 
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7 
Oey"? = res 


d.h. wir vertauschen die im mit Fourier. 
integration. 


Im folgenden benutzen wir die Abkiirzung 


Die Differentialgleichung (2.16) wird nun durch vie 


= etme (w, 2, r) 


m,n 
gelöst, wenn die ®#’” der gewöhnlichen linearen Differentialgleichung 


genügen. Ihre allgemeine Lösung ist 


r 
= (— 1)" -y, (x 7) jr r? r) dr 


Aus der Theorie de durch (2.6) definierten y-Funktionen entnimmt man nun, 
daß sie folgende Darstellungen gestatten: 


sin 

h=k+1,k>0 
h= k x 0 

wobei n; wohlbestimmte numerische Koeffizienten sind, und 


= cy o* + +---; 


diese Potenzrechenentwicklung ist fiir nichtnegative k überall, für negative k 
überall mit Ausnahme von o = 0 konvergent, wo sich ein Pol der Ordnung 
|k| befindet. 


Wegen der letzten Tatsache müssen wir in (2.20) für die Konstante B den 
Wert Null wählen: 
B=0. (2.23) 


Man kann ferner zeigen, daß die Wahl A = oo in (2.20) mit der Entscheidung 
für die halbe Summe von retardierter und avancierter Lösung der Ausgangs- 
gleichung (2.8) identisch ist. Tatsächlich wird aber aus folgendem Grunde 
diese Lösung im allgemeinen nicht existieren. Geht man nämlich mit (2.21) 
in (2.13) ein, so erkennt man, daß die Funktionen @#’”* in der Form | 


endl. Zahl sin (ar) sin (Br) 

v [2 9% 
darstellbar sind. Setzt man diesen Ausdruck in das erste Integral von (2.20) 
und benutzt für Y_„,1, wiederum die Darstellung (2.21), so erhält man 
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Hierbei wird das erste Integral der rechten Seite für A — co logarithmisch 
divergieren, während das zweite konvergiert. Somit verbietet sich dieser 
Grenzübergang, wenn nicht zufälligerweise 7“*™ verschwindet. 


Aus diesem Grunde addieren wir zu (2.20) (mit B = 0) eine Lösung der zu 
(2.19) gehörigen homogenen Gleichung, welche so gewählt wird, daß nunmehr 
der Grenzübergang A — oo ausgeführt werden kann. Die so entstehende Lö- 
sung von (2.20) soll lauten: 


r 
1/a(@, 2) 


A 
= lim (er: + (-1*ra y (ar) f N (2.25) 


In Zukunft werden wir den Stern * fortlassen. TE 
Diese Wahl der Lösung ist nur sinnvoll, wenn bei ihr die Koordinaten- 
bedingung (2.2) erfüllt ist. Es läßt sich zeigen, daß das bereits in der Form 


‚r = 0 (2.26) 
der Fall ist, so daß die Einführung der fourierfrequenzabhängigen Größe 
o in (2.25) statthaft ist. 

Wir skizzieren den Beweisgang. Zunächst betrachtet man neben (2.16) die 
Differentialgleichung 
Dt = cite Pu em, 
m,n 


und löst sie nach genau derselben Integrationsvorschrift:  __ 
y= = Pre at Hr» (2.28 a) 
r 
Hr» fim, (—1)"« y,(ar) f em. r? (x r) dr 


A 


(2.28b) 
= P_(n+1) J n n f4 ry’ Aa? 


l/o 
wobei 5", die n“”, in (2.25) entsprechende Größe ist. Zweitens läßt sich für 


jedes u, welches dieselbe Struktur wie j“” hat, also insbesondere für r — oo min- 
destens wie 1/r? verschwindet, durch explizite Berechnung zeigen, daß aus 


bei strikter Beachtung der Integrationsvorschrift die Beziehung a iA 


Yr = Wy 


folgt. (Wir flechten ein, daß unter denselben Voraussetzungen dann aus 
Dw= Div die Beziehung w = v folgt.) 
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Es gilt also 
| 
SchlieBlich zeigt man explizites daß 


Wegen des identischen Verschwindens muß nach (2.27) A” = 0, insbesondere 
auch 7”; = 0 sein. Sodann verschwindet nach (2.28a, b) auch y“? identisch 
und es folgt die Koordinatenbedingung (2.26). 

Die explizite Auswertung des Ausdrucks (2.25) gestaltet sich bei schul- 
gerechter Handhabung äußerst langwierig. Daher ist es von Nutzen, ein ein- 
faches indirektes Verfahren anzugeben. 

Man verfolge die notwendigen Rechnungen qualitativ unter Zugrunde- 
legung von Formeln über die y-Funktion, welche genaueren Aufschluß über 
deren Verhalten als (2,21), (2,22) und folglich auch über die Struktur von 
(2,24) geben. 

Im Endergebnis wird man erkennen, daß die Funktionen ©”, die folgende 
Struktur haben, wenn fourierfrequenzabhängige Koeffi- 
zienten sind und 

1 fiir gerade n 


Dem — sin(a r+ 5) + + -- ) 
+ sin (Br + + emp (2,29) 


7 uvm | vm 
+ (a r) 
— (— (ar) [Si 2ar—Si2or—Si 


Es handelt sich in dieser Formel um endliche Summen. Ci bzw. Si bezeichnen 
den Integralkosinus bzw. den Integralsinus. Zweckmäßigerweise wählt man 


är die Größe 
zunächst für die Größe o den Ausdruck BRETTEN...‘ 


wobei log y die Eulersche Konstante ist. Bei dieser Wahl wird nämlich mit 


Ci2ar+lgor— Cilor— C2or— C2Rr. (2.32) 


Man hat nun mit (2.29), (2.32) in (2.19) einzugehen und erhält zur Bestimmung 
der noch unbekannten Koeffizienten ein System linearer Gleichungen. Dieses 
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wird besonders einfach, wenn man bemerkt, daß £*”% sofort zu bestimmen ist. 


Es ist einfach proportional r””„ und kann u zur Inhomogenität des 
linearen Gleichungssystems gezählt werden. Wenn wir in der folgenden Formel a 
die Indizes u,v, m, n unterdrücken, erhalten wir mit — 


(- 1) 2 x L, = 
(2—N) L, —(—1)"40 Ly | 
(2.33) (6— N) L, +(-1)"6axZ, = 
400% 


(2 — N) 1, —(— 1)" = 

6— NL, = In 


Diese Gleichungssysteme brechen nach endlich vielen Schritten ab. In Anbe- 
tracht ihrer auBerordentlich einfachen Struktur kann man abschlieBend sagen, 
daß die Berechnung von yy” praktisch auf die Bestimmung der Größen 2 0 


kn? 
I fa welche rechts in (2.33) stehen, hinauläuft. 
§ 3. Eine spezielle Wahl für ar 


ichnet W eine singularitätenfreie Wellenfunktion 0 W = 0, so ist 


0 We 


Ta 


eine die Kae erfüllende Lösung von (2.3). Wir wählen 
+ co 


do yo(w r) ei®t 
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= 


wobei A eine positive REN ist, Es empfiehlt sich, die Größe von we y 
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ys — — 3 t PV 
cos @ ‘ Yo (@ 7) 

— — sin wt (wr), 


während die übrigen v“” gleich Null zu setzen sind. Das hiermit berechnete 
ı9“” verschwindet (bei endlichen ¢) für r > co wie 1/r4. 


Um zu zeigen, daß der Riemannsche Krümmungstensor nicht identisch 
verschwindet, betrachten wir R,,,, in der ersten Näherug. 


Es ist in unserem Falle 

a = 


Fir den von den Winkeln unabhangigen Anteil, dessen Betrachtung ausreicht, 
erhalten wir mit (3.4) 


co 

= cos wt yp (w 1) dw = 0. % 
Da sich R,;., von ,R,,,, durch Terme mit höheren Potenzen von e unter- 
scheidet, wird man schließen, daß auch der gesamte Krümmungstensor nicht 
identisch verschwindet. 


Da wir in diesem Paragraphen im Unterschied zum vorigen die reelle 
Schreibweise benutzen, schreiben wir mit 


ay"” = wh e-Ax (3.5) 
an Stelle von (2.15) i 


co co 


Es ergibt sich 


wer’ = wi (3.7) 


(2) + “| 


= (mp) Pm (at) (@, Q, 7) 
bo de 
Your) Q) = — wih (@, — 2) (3.9) 


uy 


gilt, so daB es geniigt, wey, zu berechnen, was auf Grund des im vorigen Para- 
graphen geschilderten Verfahrens geschehen kann. Dabei entsprechen die 
Funktionen Poe den obigen ®#’”. Ferner wird von (2.32) Gebrauch gemacht, 
d.h. von (2.30), welches jetzt durch zwei Formeln oj.) = ‘yee 
0(~) ara zu ersetzen ist, die implizit in die Berechnung von wf.) 
bzw. wis) eingehen. 


Zunächst zeigt man, daß ,q”” für endliche r singularitätenfrei ist, d.h, 
daß das Integral (3.6) unter dieser Voraussetzung nirgends divergiert. 


Eine Divergenz wegen des unendlichen Integrationsbereiches wird durch 
den Faktor e~-4*im Integranden verhindert. Man braucht also nur noch nach- 
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zuweisen, daß der Integrand selber keine Singularitäten aufweist. Das ge- 
schieht, indem man ihn mit Hilfe der Formeln (2.20), (2.13), (2.22) analysiert. 
Wir verzichten darauf, auf Einzelheiten einzugehen, und bemerken nur, daß 
(3.4)| wir durch unsere Wahl von 0.) bzw. oi) einen möglichen Term, der sich in 
der Umgebung von r = 0 wie log r verhält, ausgeschlossen haben. . 


inete § 4. Das asymptotische Verhalten von ,g“” 

tied Die Abschätzung des asymptotischen Verhaltens von ,g“” erfordert eine 
detaillierte Analyse der Fourierintegrale (3.6) für große r. Da ein in w und 
Q antisymmetrischer Anteil von w”” keinen Beitrag zur (3.6) liefert, kann 
w#’ von vornherein symmetrisch in @ und 22 angenommen werden. Dann 
kann die Integration statt über die positive m — 2 Viertelebene auch zweimal 
über die Achtelebene wm > 2 > 0 geführt werden. Dementsprechend schreiben 

R wir gemäß (3.6), (3.7 

ient, | wir gemäß (3.6), (3.7) 


= / dx 4* dp win + dp J da wis, 
Die Analyse dieser Ausdrücke wurde unter der Voraussetzung endlicher t, 
ıter-  d.h.t <r, durchgeführt. Außerdem wurden numerische Koeffizienten syste- 
jicht | matisch außer acht gelassen. Daher begnügen wir uns mit der Feststellung, 
daß die untersuchten Terme ‚nicht schlechter‘ als gewisse Funktionen von 
oelle | 7 für r > co verschwinden, wobei ,,gut*‘ ein möglichst rasches Verschwinden 
bedeutet. 


Alle Terme, die besser als 1/r verschwinden, werden als für das asympto- 
tische Verhalten von ,g“” unerheblich angesehen. Da die notwendigen Uber- 
legungen sehr langwierig sind, müssen wir uns an dieser Stelle darauf beschrän- 
(3.6) ken, die Ergebnisse mitzuteilen. 


Der erste Term rechts in (4.1) erweist sich als für das asympto- 
3.1) tische Verhalten von ,g”” unerheblich. 


Es bleiben also gemäß (3.7) bis (3.9) Ausdrücke der Gestalt 
COS 


3.9) | zu diskutieren. Hierbei hat Finn die Gestalt von (2.29), (2.32), wobei aber 
wegen (3.9) 2 durch — 2 ersetzt (und somit x und # vertauscht) zu denken ist. 


ara- 
die Eine genauere Analyse der in (2.29) auftretenden fourierfrequenzab- 
cht hängigen Koeffizienten /{’" und Li" zeigt dann: 

’ 
We Die l#’” enthaltenden Terme von Ft?” sind für das asympto- 
kn (B)n 


„ | tische Verhalten von (4.2) unerheblich. Die Z#’”” enthaltenden 
Terme von Féyn tragen zu (4.2) mit Termen bei, die sich nicht 
schlechter als 1/r für r— oo verhalten. 

Den größten Aufwand erfordert naturgemäß die Analyse der Beiträge zu 
(4.2), die von den Termen in (2.29), (2.32) herrühren, welche höhere Trans- 
rch zendente enthalten. Durch zweckmäßige Umformungen wird man erkennen, 
‚ch- daß diese Beiträge die nachstehend aufgeführte Struktur haben, wobei Sym- 
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bole folgender Bedeutung verwendet werden: 


p=0...9+3 
m=0...v+3—u 
h=1...? 


oA cos co 
sin (# R) cos 


é sin (nr) 4% dy. 


Diese Beiträge erweisen sich als für das asymptotische Verhalten 
von ,g“” unerheblich. 


4 sin (8 R+ 5) oo 
Diese Beiträge zu ‚g“’ verhalten sich asymptotisch nicht schlech - 
ter als 1/r. 

n 

e-4B | sin 4) e A ndy (— 1)3-” e468 f r e4: de 
0 2B 


Diese Beiträge zu ,9“ verhalten sich ebenfalls asymptotisch nicht 
schlechter als 1/r. 


ef f cos 4 r — 1 + ( - 1)3-” 


= 


é Py: > 
Das asymptotische Verhalten dieser Beiträge zu ,g“ ist nicht 
schlechter als 
a log r/r + b/r 
(a,b = const), wobei a nur dann von Null verschieden ist, wenn in IV» = 3 ge- 
nommen wird. 

Das Auftreten von Termen, die sich asymptotisch wie log r/r verhalten, ist 
ein physikalisch höchst unerwünschter Umstand. Man muß insbesondere 
befürchten, daß solche Terme die Konsistenz der nächsten Näherung zerstören 
werden. Daher liegt es nahe zu versuchen, durch Addition geeigneter Lösungen 
der homogenen Wellengleichung zu ,g”’ die fragwürdigen Terme zu neutrali- 
sieren. Um mit der Erfüllung der Koordinatenbedingung in Einklang zu 
bleiben, geschieht das am besten, indem man über die Größe o,;, neu verfügt. 
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Wir schreiben daher 


(8) neu = 0*(&, P). 
Wenn man 


wählt, kann man zeigen, daß nunmehr asymptotische Terme log r/r nicht mehr 
auftreten werden. Ferner erweist sich, daß die Singularitätenfreiheit unserer 
Lösung nicht zerstört wird. 

Unser Ergebnis ist also: Esexistiertinzweiter Näherungeinenicht- 
triviale singularitätenfreie Lösung der Feldgleichungen R,,= 0 
mit dem asy mptotischen Verhalten für große r und endliche t 


. "> 

wobei U#’eine Funktion von’ d und o allein ist. (Während unser Ab- 

schätzungsverfahren nur angibt, daß das asymptotische Verhalten nicht 

schlechter ist als in (4.3) angegeben, haben wir durch explizite Berechnung 

besonders einfache Terme, die durch andere nicht neutralisiert werden können, 

nachgewiesen, daß das asymptotische Verhalten auch nicht besser sein kann.) 


Da für das hier betrachtete Feld T} — 0 ist, wird die Gesamtenergie durch 
E =f (Ti+ ti) dr = fth@e (5.1) 


gegeben. Dieser Ausdruck soll in der niedrigsten Näherung in ¢ berechnet 


(4.3) 


werden, wobei die Einsteinsche Definition von t, zugrunde gelegt wird. 

Man findet 

2x 1 1 ¢ 


wobei x die Gravitationskonstante, ‚9 = 7.819”, und ,G,, die erste Näherung 
von ist, welches durch g*” = 6; definiert ist. 
Vergleicht man diesen Ausdruck mit 1. c.6) Formel (3.10), so findet man 
unter Benutzung der Koordinatenbedingung und J,g?=0 


ati= DO räumliche Divergenz. (5.3) 
2 


Da die räumliche Divergenz wegen des raschen Verschwindens von ,g*° keinen 
Beitrag zum Integral (5.1) liefert, kann man schreiben: 


7 
- Asg**) Bx = — fs + 


% 


r>o 


do = sin 8 dd dq. 
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Es wurde die Koordinatenbedingung in der Form = — benutzt. 
Auf dieselbe Formel gelangt man mit Hilfe der Freudschen Superpotentiale. 
Mit Hilfe der in dieser Arbeit geschilderten Methoden berechnet man 

3 


E=z At 


In unserer Näherung hat das Feld also eine positive Gesamtenergie. 


Herrn Prof. Dr. A. Papapetrou danke ich herzlich fiir die Anregung 
dieser Arbeit sowie für sein förderndes Interesse bei der Durchführung der- 
selben. 


Greifswald, Institut für theoretische Physik der Universität. 


un ss Bei der Redaktion eingegangen am 26. Juni 1961. 
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Inhaltsiibersicht 


Ausgehend von einem in einer vorangegangenen Arbeit ausfiihrlich dis- 
kutierten Variationsprinzip werden Näherungslösungen aufgesucht, die einen 
rotierenden Kern modellmäßig, aber im Rahmen des gewählten Modells mit 
optimaler Genauigkeit beschreiben. Als Modellnäherungen werden Hartree- 
Fock-Verfahren und Brueckner-Theorie untersucht. Für das Kernträg- 
heitsmoment werden dabei Formeln abgeleitet, die eine Verallgemeinerung der 
bekannten Formel des cranking.model darstellen. Die Rechnungen weisen in 
Verbindung mit Ergebnissen von Amado und Brueckner darauf hin, daß 
die Hartree-Fock-Näherung auf das starre Trägheitsmoment und damit 
nicht zur Übereinstimmung mit den experimentellen Werten führt. Das 
Bruecknersche Verfahren liefert Zusatzterme, über deren Größe jedoch noch 
nichts bekannt ist. are 


Die vorliegende Arbeit enthält weitere Untersuchungen zum Problem der 
Trägheitsmomente bei Kernrotationen, nachdem in einer vorangegangenen 
Arbeit?) bereits allgemeinere Überlegungen zu diesem Thema angestellt wurden. 
Hinweise auf einzelne Kapitel dieser Arbeit erhalten folgende kurze Form: 
14 bezeichnet z. B. das Kapitel 4 der Arbeit). Entsprechend werden Formel- 
zitate in der Form (14.6) angeführt. 


Nach einer Darlegung des allgemeinen Gedankenganges in Kapitel 1 ent- 
halten die Kapitel 2 und 3 die Theorie eines rotierenden Kerns in Hartree- 
Fock-Näherung. Dabei wird zunächst kein Gebrauch von der bekannten 
Tatsache gemacht, daß sich das Hartree-Fock-Verfahren direkt als Varia- 
tionsaufgabe formulieren läßt. Diese gleichwertige Möglichkeit wird in Kapitel 4 
dargestellt, die zusätzlich gewonnenen Erkenntnisse werden in Kapitel 5 bei 
der genäherten Auswertung der für das Trägheitsmoment abgeleiteten Formel 
verwendet. In Kapitel 6 wird die Theorie am Beispiel der Schwerpunkts- 
bewegung getestet, wo das Resultat von vornherein bekannt ist. Die Kapitel 7 
und 8 behandeln die bei Anwendung der erforder- 


lichen Modifikationen. 


1) Auszug aus der Dissertation des Verfassers, Dresden 1961. 
2) W.Macke u. R. Lenk, Ann. Physik 9, 1 (1961). fc 


dr 
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= 1. Gestörtes Modellproblem 
re i In (1) ist eine Theorie der Kernrotationen entwickelt worden, bei der das G 
Variationsprinzip 
—®3—AN] = 0 (1) 
zur Untersuchung eines rotierenden Kerns benutzt wurde. Die dabei gewon-| des I 
nenen Erkenntnisse über die Deutung dieses Verfahrens, die Einführung eines 
Trägheitstensors und dessen Entartung werden als bekannt vorausgesetzt. 
Außerdem wurde in 15 eine sehr allgemeine Lösung zu (1) berechnet, von der 
im folgenden jedoch kein Gebrauch gemacht wird. Statt dessen beschränken! aus ¢ 
wir uns hier von vornherein auf Näherungslösungen zu (1), die aber über die in! gura‘ 
13 verwendete grobe Näherung H = H® hinausgehen. 

Ein Näherungsverfahren kann immer so aufgefaßt werden, daß dem 
Hamiltonoperator H ein Modelloperator H® zugeordnet wird. Dadurch wird! aus 
eine Zerlegung 

H = H®+ U® (2) 
definiert, die Näherung besteht in der Vernachlässigung der Störung U". 


Wird die in (1) enthaltene Nebenbedingung vorgegebenen Drehimpulses $ | Dis 
bei der Durchführung des Näherungsverfahrens nicht gestellt, dann beschreibt | in (4 
H° in (2) einen— nach Voraussetzung deformierten — Modellkern ohne Rotation, | werd 


wie in 12 ausführlich diskutiert wurde. das | 
Sollen dagegen Näherungslösungen zu (1) aufgesucht werden, so ist H I 


durch H — @ $ zu ersetzen. Die Anwendung des gleichen Näherungsverfahrens, 
das ohne die Nebenbedingung zur Zerlegung (2) führt, liefert dann bei Anwen- 


dung auf H — @ entsprechend er defin 


H = + (3) | zusti 


eine durch den Zusatzterm — @ & geänderte Aufspaltung von H in Modell- == 
operator H® und Störung U®. Die Näherung besteht nun in der Vernachlässi- et ‘ 
gung von U®. Für ein Einteilchenmodell besteht der Unterschied zwischen | “°° 
H® und H® qualitativ in folgendem: In H® ist ein mittleres Einteilchenpotential 
V° enthalten. Dieses ist nicht willkürlich wählbar, sondern durch das gewählte 
Näherungsverfahren festgelegt. Nur bei diesem bestimmten Potential ist die } »egi; 


Vernachlässigung der Störung U® gerechtfertigt, ist das Verfahren selbst- | 7ust 
konsistent. Diese Konsistenz wird aber durch die Rotation des Kerns zerstört. | ge fj 
Das Potential V° stellt einen Ersatz für die Nukleonenwechselwirkungen im 
ruhenden Kern dar, nicht aber für die im rotierenden Kern. Denn das Einteil- 
chenpotential entsteht durch Mittelung der Zweiteilchenwechselwirkungen über 
die Zustände der einzelnen Nukleonen, diese Zustände werden aber durch die Dur: 


Rotation des Kerns, d. h. durch den Term — IR: in (3), verändert. Das Poten- v0; 
tial V® des rotierenden Kerns wird also vom Potential V° des ruhenden Kerns 
verschieden sein. Es ist daher nicht richtig, beim Aufsuchen von Näherungs- 
lösungen zu (1) einfach H durch H® zu ersetzen, sondern das jeweilige Nähe- folg 
rungsverfahren muB entsprechend (3) neu durchgefiihrt werden. Eine der 
Hauptaufgaben dieser Arbeit besteht deshalb darin, H® zu berechnen, wobei 
H® als bekannt vorausgesetzt wird. H® —@ & beschreibt dann einen rotierenden In ( 
selbstkonsistenten Modellkern, mit den zugehörigen Lösungen ®® können die 
erforderlichen Erwartungswerte von Drehimpuls und Energie gebildet werden. 
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2. Hartree-Fock -Néherung 
Gegeben sei ein von Einteilchenfunktionen 


+77 Mm = + e, Pr (Pn: Pm) (1) 


die zu einem bestimmten mittleren Potential V° gehéren. Die Modellfunktionen 
des Kerns sind Determinanten 


= (98 (m) | = (sj; . . . (2) 


aus den g?, N i die Gesamtheit der besetzten Zustände, die Konfi- 
guration des Kerns. Der Modelloperator hat die Form 
(3) 


aus (1) folgt dann sofort C cones 


rt asst das pte’ \, D9 Dy). 


D® ist nach (3) ein Operator, der in den ®\ nur Diagonalelemente besitzt und 

in (4) noch zur geeigneten Festlegung der Energieeigenwerte von H® benutzt 

werden kann. Im folgenden wird nur eine einzige Forderung, namlich (10), an 

das Spektrum von H® gestellt, für D® genügt dann einfach eine Konstante. 
Das Potential V° wird durch 


D = 06) 
m Las 


definiert. Hier und im folgenden sollen |k) und |/> stets beliebige Einteilchen- 
zustände, |n) und |m) die in einer bestimmten, ausgezeichneten Konfiguration 
(Grundzustand) besetzten Zustände sein. Die Wellenfunktion zu dieser Kon- 
figuration wird mit ®} bezeichnet. Die Elemente < |v| > der Wechselwirkung 
sollen den Austauschterm mit enthalten. Die Störung 


Hunde 
re (6) 


D 
ij 


besitzt auf Grund von (5) keine Übergangselemente, bei denen sich nur ein 
Zustand ändert, während die übrigen Zustände sämtlich zur bei der P Potential- 
definition in (5) ausgezeichneten Konfiguration gehören: N 


<N |U°| N’y = 0, falls 
= (7) 
N = (n,..1..%,). 


Durch geeignete Wahl von D® kann ferner das Diagonalelement der Störung 
U° im Grundzustand zum Verschwinden gebracht werden. Aus 


= = 0 | (8) 


1 
<n m |v| n m) (9) 
nm 


In (4) wir durch D® die in der Summe über alle Einteilchenenergien e} doppelt 
enthaltene Wechselwirkungsenergie einmal wieder abgezogen. uhr 
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: Aus (8) folgt, daß sich beim Hartree-Fock- Verfahren die Modellenergie Zur 
nach Störun 
Ey = = <H) (10) 


auch als Erwartungswert des exakten Hamiltonoperators H mit der Modell. 
funktion ©? schreiben läßt. Aus diesem Grund kann man ®0 auch als Variations. 
ansatz für die Wellenfunktion betrachten und das gesamte Verfahren aus der) Nunme 
Variationsforderung 


ableiten. Die Gleichungen (10) und (11) gelten nicht mehr in der Brueckner- kerechr 
schen Theorie. Daher wird die Variationseigenschaft (11) zunächst nicht ver- dann we; 
wendet, da dann das Verfahren einfacher verallgemeinert werden kann. die in ( 

Die Anwendung des Hartree-Fock-Verfahrens auf H — @  fiihrt zu 
einem neuen Modelloperator H®, einem geändertem Potential V® und den zuge- 
hérigen Wellenfunktionen ®¥ und g;,. Zwischen all diesen Größen bestehen 
genau die gleichen Beziehungen wie zwischen den entsprechenden ungestörten 
Größen. Dementsprechend gilt Das fü 
He — H°= (v} — D’). (12) 
v 
Nach der Bemerkung im Anschluß an (4) können für D® und D® einfach 
Konstanten gesetzt werden. Die Nichtdiagonalelemente von H® — H® werden Erst di 
dann durch die Potentialdifferenzen in (12) bestimmt. Offensichtlich existieren daß (2 
nur Übergangselemente von H® — H,, bei denen sich lediglich ein einziger ’ 
Einteilchenzustand ändert. Die erste Näherung hinsichtlich @ definiert nach 
ow abe tensor | 
eine Größe J, von der Dimension eines Drehimpulses, der in barf gilt fi 
H® —@ = (3 + (14) | #13) 


zum gewohnlichen Drehimpuls hinzutritt. Die in (13) und (14) vernachlassigten 

Terme sind zwar von höherer Ordnung in w, gestatten aber entsprechend (12) 
ebenso wie §, nur Einteilenenübergänge. In der zu (14) gehörigen Wellenfunk- | wobei 

tion im übı 

P | Marit (15) 

€, = — 3 =3+2 


fehlen daher Anteile der Ordnung w?, bei denen gegenüber ®® nur ein einziger 


Einteilchenzustand geändert ist. 
Zur Bildung des Drehimpuls-Erwartungswertes ist zunächst nur die erste | __ 
Näherung von (15) erforderlich. Das Ergebnis Die Sc 

I= (0°, JG") =<} + (16) |" 
führt sofort auf den Trägheitstensor 


Dabei treten gegenüber (13.14) die %,-Terme auf, die dc ie die Störung des | Mit de: 
Modellproblems zustande kommen. Der Index % bedeutet, daß 03 in (16) über 


i 
= 
| 
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Zur Berechnung der Energie wird ausgenutzt, daß das Diagonalelement der | 
Störung U® mit der gestörten Grundzustandsfunktion ®2 analog zu (8) ver- 
schwinden muß, es gilt also 

U? = (9, = 0. (18) 
vite 


Nunmehr kann die Energie des rotierenden Kerns nach 
(©), [H° + 
berechnet werden. Bei der Auswertung von (16) bis zur Ordnung w? verschwin- 
den wegen (7) und (8) eine Reihe von U°-Termen, wegen (7) liefern insbesondere 
die in (15) fehlenden Anteile der Wellenfunktion keinen Beitrag. Man erhält 
H® = 4+ [H®— <A] iy 
) + <@Y’ — U® + ©. ec. (20) 
Das führt auf das Trägheitsmoment past 
+ Yo! Ve! 
Erst die genauere Untersuchung dieses Ausdrucks im nächsten Abschnitt zeigt, 
daß (21) mit dem einfacheren Ausdruck (17) identisch ist. dott 5% 
st 


ind 
3. Trägheitsmoment und gestörtes Potential 


He = H = (19) 


In (2.17) und (2.21) wurden zwei Ausdrücke #3 und 6, für den Tragheits- 
tensor definiert, die im folgenden genauer untersucht werden sollen. Zunächst — 
gilt für die Übergangselemente des in beiden Formeln auftretenden $, nach 
(2.13) 

\ Cw \ | 
wobei |k> und |!) die beiden Einteilchenzustände sind, durch die sich die 
im übrigen beliebigen Konfigurationen N und N’ unterscheiden. Definiert — 
man noch analog zu (2.13) einen Drehimpuls j, durch 


. hör Ta 
so wird aus (1) einfach 


Die Schrédingergleichung fiir die gestérten Einteilchenfunktionen lautet 
nun 


[t+ V°— aj] =[t—aj'] = =itio (4) 
Dabei ist im ersten Gleichheitszeichen die Näherung (2) für V® eingesetzt. | 
Mit den oben eingeführten Größen lautet 63 ausführlich 


03 = & <n | k> |j| m> + ef, = ef — (5) 
nk Een 
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Bei der Auswertung von 6, ist zu beachten, daß U® nach (2.7) keine von 
ausgehenden Einteilchenübergänge gestattet. Der Term 


Ue 


vereinfacht sich daher wie in (6) angegeben, denn es sind zwei 3’-Ubergainge 
erforderlich, um die durch U® angeregten zwei Teilchen wieder in den Grund. 
zustand zurückzuführen. Ausführlich lautet (6) 

Ekisnm 
|i’ | m>——<k |i’ | m> im 
Een Ein 
= (ee + ef) — (en + em). (7) 


Die Summationen in (7) sind sämtlich voneinander unabhängig. Der Faktor 
1/4 berücksichtigt, daß durch Vertauschung von n und m sowie von k und | 
bei ununterscheidbaren Teilchen keine neuen Zustände des Gesamtsystem: 
erzeugt werden. Die Terme in der eckigen Klammer von (7) entstehen durch! 
die vier verschiedenen Möglichkeiten, den durch U® entstandenen angeregten 
Zustand mit dem Grundzustand zu verbinden. Der Vorzeichenwechsel in den 
beiden letzten Termen entsteht durch die Vertauschung von |m) und |n) in 
®0. Durch nochmalige Vertauschung der Summationsindizes m und n gehen 


sie bei Beachtung von (mn |v| kly = — <nm|v| kl in die ersten beiden Terme 
über. Benutzt man noch die Energienenner-Identität 
1 1 1 1 
| ) = 7777: seal 
Erıınm \Ekn Exn Eim ‘i 
so wird aus (7) i 
kn im 
Insgesamt erhält 6, nun die Form 
nk Een mi Elm 
+5 <km nl +e c. (10) 
ml 


Die Summationen in (5) und (10) laufen unabhängig über alle besetzten Zu- 
stände |m) und |n) und alle unbesetzten Zustände |k) und |l). Die letzte 
Einschränkung ist jedoch unnötig. Löst man nämlich (4) durch Störungsrech- 
nung, bildet die ®” aus den gestörten gy; und berechnet dann die Erwartungs- 
werte von Drehimpuls und Energie, so erhält man wieder (5) und (10). Da aber 
bei der Berechnung der gy» alte Zustände |k>, auch die besetzten, berücksichtigt 
werden, laufen die Summationen in (5) und (10) dann über alle |k) und |)). 
Man kann durch eine einfache Rechnung auch direkt zeigen, daß die Summa- 
tion über die besetzten Zustände |k) und |T) keine Beiträge zu 43 und 6, 
liefert. 
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2 


Eine Bestimmungsgleichung für die in (2) definierte Größe j, erhält man aus © 

der Definitionsgleichung wih 


anh tin oala = ck <r m |v| sm) <s 12 

m 

für die Elemente des gestörten Potentials V®. Diese unterscheiden sich also — 
von den entsprechenden Elementen des ungestörten Potentials V° nur dadurch, 
daß die Mittelung über die Wechselwirkungen v mit den gestörten Einteilchen- 


funktionen |m) ausgeführt wird. Zur Berechnung von jy ist V® bis zur ersten . 
Ordnung in w zu entwickeln, es genügt also die Näherung fos dia 
Em 
zur Auswertung von (11). Das Ergebnis 
fix Vit lass: Varck 1 
führt mit (2) auf die gesuchte Gleichung (14) 


, y 1 ] , \ , 1 \ 
= & <m |i <kl\olk m) + |v] | m>| 
tm Im 


für jp. Wegen j’ = j + jg ist (14) vom Typ einer Integralgleichung für jp. Löst 
man diese durch Iteration, so wird j, durch eine Reihenentwicklung nach 
Potenzen von v dargestellt. Dabei entsprechen die verschiedenen Potenzen von 
v genau den einzelnen Schritten bei der Lösung des gestörten Hartree-Fock- — 
Problems: In nullter Näherung hat man die Wellenfunktionen des ungestörten 


Kerns, es ist 99, = 9° und Vo, = V°. Die Wellenfunktionen erster Näherung © 
gehören also noch zum ungestörten Potential V®, durch die Störung — aj 
unterscheiden sie sich aber von den ungestörten Funktionen, gj) enthält 
also nur den Einfluß von —@ . Berechnet man mit gaynach (11) das gestörte 
Potential Vay, so unterscheidet es sich nur um Terme, die linear in v sind, von 
V°. Erst der nächste Schritt des Verfahrens berücksichtigt bei gs) die Störung 
des Potentials in der Näherung V(j), das mit gi) berechnete verbesserte 
Potential V(s) enthält die Rückwirkung der Potentialänderung auf das Poten- 
tial selbst durch die in v quadratischen Terme. Allgemein enthält das mit den 
Wellenfunktionen (m) berechnete Potential V(n) Terme bis einschließlich der 
n-ten Ordnung in v. Bei Konvergenz des skizzierten Iterationsverfahrens wird 
man nach einem der ersten Schritte abbrechen können, das entspricht genau 

der Lösung von (14) durch Iteration und Abbrechen der entstehenden Reihe. 


Die Bestimmungsgleichung (14) für j, gestattet sofort eine wesentliche Ver- 
einfachung von (10). Die geschweifte Klammer in (10) ist nämlich auf Grund 
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von (14) einfach (k |j’ — jo! >, der Vergleich von (14), (10) und (5) zeigt also 
die Gültigkeit von 
6, = 03 = 8. (15) 


Die Erwartungswerte von Energie und Drehimpuls fiihren also auf das gleiche 
Trägheitsmoment 9, falls die Störung des mittleren Potentials durch die 
Kernrotation richtig berücksichtigt wird. 


4. Variation des Trigheitsmoments 


Wie in Kapitel 2 2 näher ausgeführt wurde, ist beim Hartree-Fock-Ver- 
fahren wegen 


<U% =U°=0 (1) 


sowohl die ungestörte Energie <H®) wie die gestörte H® als Erwartungswert des 
Hamiltonoperators H darstellbar. Die Ergebnisse der letzten Kapitel 
müssen sich daher auch erzielen lassen, wenn das grundlegende Variations- 


prinzip (1.1) direkt mit dem Variationsansatz int ow ara 


Gow 


fiir die Wellenfunktion verwendet wird. Eine ahnliche Rechnung findet sich 


auch bei Thouless3). Nach (3.4) sind die 9, durch V” bzw. j, bestimmt. Last 
man j, zunächst offen, so entstehen durch Variation von j, verschiedene Lösun- 


gen (2). Unter diesen wird die richtige durch die Forderung ete 
ausgewählt. Die in (3) auftretenden Erwartungswerte sind in den Kapiteln 2 
und 3 berechnet worden, das Ergebnis ea 
= «H — > 
0= 20, — sid 


zeigt zunächst, daß die ia io- Abhängigkeit von H —@ im Trägheits- 


fiir 6, wobei 63 und 6g in (3.5) und (3.10) explizit angegeben wurden. 3 


h Zur Durchführung der in (6) geforderten Variation wird 263 in der Form 


1 

203 = |" <n |i'| <k |i] m> + <m|j|D + ee 
Een mi Erm 
( 


3) D. J. Thouless, Nuclear Phys. 21, 225 (1960). 
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(1) (3.14) für j, führt. Die Variationsaufgabe (6) ist dieser Gleichung äquivalent, > 
t des} der zu variierende Ausdruck 6 stimmt mit dem gesuchten Tragheitsmoment 5 
pitel | überein, falls die zu variierende Größe j,(t, p) die Variationsaufgabe löst, jg 
ions- | selbst hat dann die aus (3.2) ersichtliche Bedeutung. 

Die Ergebnisse dieses Kapitels sind nützlich für eine genäherte Auswertung 
des Verfahrens. Die Bestimmungsgleichung (3.14) für j,ist praktisch unlösbar, _ 
(2) | daher sind bei der Auswertung von ® noch weitere — über die von vornherein ee 


sun- | wert —— % für die richtige Lösung der Variationsaufgabe ein Minimum u 
besitzt, hat @ @ entsprechend (4) ein Maximum. der in (13.28) disku- 


(3) | Theorie erhalten bleibt, ist er 


= gle x o)?. teat (9) 

Näherungslösungen der Gleichung (3.14) fiir jg führen also immer zu einem 
größeren Trägheitsmoment g, als man für die exakte Lösung von (3.14) er- 
(4) | hielte. 
5. Genäherte Auswertung 

(5) In (4.5) wurde ein Ausdruck fiir den Tragheitstensor angegeben, der bei 
eits. | Tichtiger Berücksichtigung der Differenz zwischen gestörtem und ungestörtem 
Potential einen stationären Wert annimmt. In einfachster Näherung wird 


io =0 ji =j (1) = 
(6) gesetzt. Man erhält für das Trägheitsmoment um eine Achse senkrecht zur _ 


Figurenachse den Näherungswert 


= 


1 
| | k> <k | | 


|| 
Man erkennt ohne weiteres, daß 66 — 0 wieder zu der Bestimmungsgleichung i 
42 
NZ 
| 
= 
2 | 
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7, ist die Komponente von j in Richtung der Drehachse. Die Näherung (1) und 
(2) entspricht einer Entwicklung des Trägheitsmomentes nach Potenzen von 
v, die mit den in v linearen Termen abgebrochen wird. Denn der bei (1) gemachte 
Fehler ist von der Ordnung v, der bei g muß deshalb wegen der Extremaleigen- 
schaft (4.6) von @ von der Ordnung v? sein. Man erhält daher das Ergebnis (2) 
auch durch Reihenentwicklungen von 03 oder 6g bis zu den in v linearen Glie- 
dern. 


Der Ausdruck (2) für g ist bereits von Amado und Bruecknert) für den 
rechnerisch einfachsten Fall, nämlich für spinlose Teilchen im periodischen 
Volumen, untersucht worden. Die Einteilchenfunktionen 9% sind dann ebene 
Wellen, die Summationen über Zustände kénnen meist als Integrationen im 
Impulsraum ausgeführt werden. Außerdem wird j, = !,. Unter diesen Vor- 
aussetzungen wird in 4) die Beziehung 


1 
er > <n dı | k> „0 m dy | b k l |v | n m) Met +73 
it 


= 27% n|l,| k ie m |v) km) — <nm |v| nm>] <k |l, | n> 


nkm kn 


bewiesen. Die Rechnung wurde nachgeprüft und der Beweis erweitert auf den 
Fall eines Nukleonengases mit Neutronenüberschuß. Die Integrationen laufen 
dann über zwei Fermikugeln verschiedener Größe, (3) bleibt dabei gültig. 
Mit (3) erhält (2) die einfache Form 
Vel eS — 70 


Eon Een 


Bis auf Terme der ohehin vernachlässigten Ordnung »? ist (4) mit 


1 
= 22 <n |l, 5 
nk ( k | k> <n n>) | | » q t ( ) 


gleichwertig. In (5) stehen im Energienenner nur noch die Differenzen der 
kinetischen Energien, während die ursprünglichen Nenner ¢° durch die Impuls- 
abhängigkeit des Potentials V° geändert sind. Wie ebenfalls in *) bewiesen 
wurde, ist (5) das Trägheitsmoment g,,; für starre Rotation. Rechnungen mit 
realistischeren Wellenfunktionen g}, insbesondere Rechnungen für endliche 
axialsymmetrische Kerne, liegen nicht vor. 

Die grobe Näherung (1) läßt sich ohne erhöhten Rechenaufwand sofort ver- 
bessern durch den Ansatz 


j=u(lt+a)j (6) 
mit einem freien Variationsparameter x. Die in (6) enthaltene einfache Pro- 
portionalität zwischen j, und j muß nicht richtig sein, verglichen mit j, = 
bedeutet die in (6) zugelassene Variationsfreiheit aber jedenfalls eine Ver- 
besserung des Verfahrens. Beschränkt man sich auf eine quadratische Ab- 


4) R.D. Amado u. K. A. Brueckner, Physic. Rev. 115, 778 (1959). 
A> 
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) und § hängigkeit des Potentials V° vom Impuls, so lauten die Energienenner a 
( 9 2 N m 2 
achte Ehn (pi — Pn) (1+ n) == (1+ N) Ekin = m* Ekin (9) 
igen 
Es ist also 
> den <k | V°| k> — <n | n> = (8) 
“a Fir das zu a Trägheitsmoment erhält man nach (4.5) unter Ver- 
n im | wendung von ) bis (8) 
9¢ itn I1+n (l+n? “1l+y 1+n/' 
Für x = 0 erhält man aus (9) wieder ide sind rate 
g (0) ait 2n 1, (10) 


Ist (1+)? 


(3) | das ist das alte Ergebnis ¢4) mit der in (7) zusätzlich eingeführten Effektiv- 
massennäherung. Den besten Wert a, für & liefert die Variation 


> A 
dg /9 54 =2(1 =0 (%) _ (11) 
alte Ist 

den 1+7 

‚ufen 


Die Variation von « bestätigt also die in (10) bereits angedeutete Tendenz 
zu = Ist- 

Abschließend sollen Fälle diskutiert werden, in denen sich die in v linearen 
(4) | Terme nicht kompensieren. (Diese Annahme ist also nicht in Einklang mit den 
speziellen Rechnungen von ®). Mit einem neuen Parameter q vor den v-Termen 
von 9, also vor dem zu n proportionalen Anteil in (9), wird dann aus (9) 


(5) 2 9 
| 9(0) _ 1+ ——— 
9. [ 1+ N (13) 
liche | Der Extremwert liegt bei ‘ 
ver- 1+n 1+(1—q)7 
Falls g in der Nähe von 1 res unterscheidet sich g/gs¢ also relativ wenig von 
(6) | 1, für „= 2/3 (d.h. nach (7) m*/m = 0,6) gilt sire Kae Re 
Ver- gigst=1 3/4 3/5 3/7. (15) 
= Die Ubersicht (15) zeigt, daB erst ein Vorzeichenwechsel von g = 1 zug = — 1 


die experimentell gemessenen Trägheitsmomente erklären würde. 


a 


ilt 
| 
a) 
| 
; 


90 Annalen der Physik. 7. Folge. Band 9. 1961 


Wie es nach der allgemeinen Diskussion von Kapitel 4 sein muß, ist das also i 


Verhältnis 

stets kleiner oder höchstens gleich 1, die Näherung « x = 0 liefert also ein zu D 


kleines Trägheitsmoment, ein verbessertes Verfahren führt zu einem größeren | tial J 


ert für g. pa 
6. Schwerpunktsbewegung 


Es soll hier noch gezeigt werden, daß die bisher entwickelte Methode bei 
Anwendung auf die Schwerpunktsbewegung zum richtigen Ergebnis führt. 
Das wurde zwar in 16 bereits allgemein, ohne Benutzung von Störungsrech- 
nung, gezeigt. Dennoch ist es interessant zu untersuchen, auf welche Weise | betra 
eine zur Behandlung der Kernrotation analoge Entwicklung nach Potenzen | Einfi 
der Geschwindigkeit v auf die Gesamtmasse der Kerns führt. Um die Rechnung | nenn 
zu vereinfachen, soll nur der Fall gleicher Teilchenmassen, also 


mM (1) 


betrachtet werden. Die gestörten Wellenfunktionen’y, unterscheiden sich ent- | wird 
sprechend (16.3) von g, nur durch einen Phasenfaktor 


Gn =e" 
Sie geniigen der Gleichung E 
i exak 
vp+J (r,p—mo)| Pn =e + + V(r, (3 
aus der man sofort 
p) = p mv) (4) 
entnimmt. Die Entwicklung von V® nach Potenzen von » ist damit ebenfalls 
bekannt, in der erforderlichen Naherung ist a ae dars 
Ve = V°— pm V°—y 
Damit erhält der analog zu (3.5) gebildete Massentensor die Form 3 
0) 
Aus der te gebr 
lassen sich aber die Matrixelemente von ® + po) als edie IR ber 


berechnen. Klassisch wäre p+p,= mt, wie ebenfalls aus (7) abgelesen 
werden kann. Dagegen hat der kanonische Impuls bei impulsabhängigem 
Potential nicht diese einfache Bedeutung. Setzt man (8) in (6) ein, so folgt Die 


: 
=) 
3 
um d 
| 
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st das also insgesamt nach 
A = - nzahl - 
(16) u I = Einheitstensor 


das erwartete Ergebnis. 

in zu Dieses Resultat hätte man formal auch mit impulsunabhängigem Poten- 
Beren # tial V(r), also p, = 0, erhalten. Bei der obigen Rechnung ist V? + V°®, p ist 
dann einmal durch p + p, zu ersetzen. Dadurch wird gerade die von der Im- 
pulsabhängigkeit des Potentials verursachte Änderung der Energienenner 
je bei kompensiert. Um das zu sehen, wird wie in (5.7) der einfache Fall 

Neise | betrachtet, in dem sich die NEN des Potentials vollständig durch 


nzen Einführung einer effektiven Masse m/1 + n beschreiben läßt. Die Energie- 
nung | nenner unterscheiden sich entsprechend 


& = (1+ 7) & e() = Energienenner mit V (r) (12) 


2m 


A) um pm Faktor 1 + 7 von denen des lokalen Potentials V(r). Dieser Faktor 

ent- | wird durch den gleichseitigen Ubergang von p zu 
P+ Po = p(l + 9) 


im Ausdruck (6) für M gerade kompensiert. 


Es bleibt noch zu zeigen, daß die in (5) definierte Größe p, tatsächlich die 
exakte Lösung der zu (3.14) vollständig Bestimmungsgleichung 
3 
(3) <k = — > |p + N - -<kllv|k' m) 
mi 14 
ab 1 (14) 
(4 + <km |v| lb ——<l |p + p,| m) 
falls 
darstellt. Zunächst kann (14) unter mc, von (8)und 
eve 
in die Form 
k\(V®, k’> + = [<km |v| ly m> —<m <kl my] = 0 
(6) (16) 
gebracht werden. Im Anhang wird dieser Ausdruck untersucht. Die Forderung 
(16) stellt sich dabei also identisch mit uote 
(7) ~< m=0 (1%) 
heraus. p, ist also dann eine Lösung von (14), wenn die in (17) enthaltene Ver- 
(8) tauschung verschwindet. Das ist trivialerweise erfüllt bei Wechselwirkungen, 
die nur vom Ort abhängen, allgemeiner für 
sen i ö 
em E u + to] = + dp, = = 0 v = vlt, Ty, Py Po). (18) 
Igt Die letzte Form von v ist notwendig, falls der Hamiltonoperator H covariant 
P 
(9) gegen Galileitransformationen r = r — v ¢ sein soll. Eine genauere Betrach- 
A tung zeigt, daß auch nur dann der einfache Ansatz (2) für die Wellenfunktionen 


4. 


aa 
i 
i 
| 
Jeu 
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Y„ das Variationsprinzip (16.1) für die Schwerpunktsbewegung löst. Natürlich 
ist (18) in allen physikalisch sinnvollen Fällen gesichert. 


7. Trägheitsmoment bei starker Zweiteilchen-Korrelation 


Bisher wurde angenommen, daß die Wechselwirkungen v die Anwendung 
des Hartree-Fock-Verfahrens gestatten. Tatsächlich sind jedoch die Kern- 
kraftpotentiale nach den modernen Vorstellungen stark singulär bei kleinen 
Abständen, sie besitzen sogenannte hard cores. Daher müssen Korrelationen 
zwischen je zwei Nukleonen sorgfältig berücksichtigt werden. Das geschieht 
in dem von Brueckner u. a.*) ausgearbeiteten Verfahren, bei dem im Ver- 
gleich mit dem Hartree-Fock-Verfahren v bei allen (H®, AA 
E®) durch die sogenannte Streumatrix such 


ersetzt wird. Der Projektionsoperator Pf in (1) erlaubt nur unbesetzte Zwischen- 


zustände (Zustände außerhalb der Fermikugel), e° ist der Energienenner für 
die Anregung zweier Teilchen. Ausführlich lautet (1) 


(nm kl = «nm |v| kl (2) 
po 1 


= bedeutet dabei eine Summation über alle Paare von Zuständen |p) und 


a> auBerhalb der Fermikugel. Die Energie erhalt die Form 


0 


= = nm |K°| n m), (3) 


wobei H® ganz analog zu (2.3) und (2.5) definiert ist, nur daß überall » durch 
K® ersetzt werden muß. Nach (3) ist die Energie E® nicht als Erwartungswert 
von H gegeben, das Diagonalelement <U®> der Störung verschwindet nicht, 
und dem Verfahren kann kein äquivalentes Variationsprinzip zugeordnet 
werden. Die gleichzeitige Bestimmung von Wellenfunktionen g,, Streumatrix 
K° und Potential V° stellt ein außerordentlich kompliziertes Problem dar, da 
mit dem üblichen Hartree-Fock-Problem der Bestimmung von 9% und V° 
das Problem (1) gekoppelt ist. Das Verfahren ist speziell auf die Energie- 
berechnung nach (3) zugeschnitten. 

Zur Behandlung eines rotierenden Kerns ist das oben besprochene Verfahren 


für den Hamiltonoperator H —@ % durchzuführen. Aus 
3 = <P 000 (4) 


kann dann der Trägheitstensor 9 entnommen werden. Der in (4) auftretende 
Erwartungswert lautet ausfiihrlich 


+> Enm|K|nm) (5) 
n nm n 


5) K. A. Brueckner u. J.L. Gammel, Physic. Rev. 109, 1023 (1958) und dort 
zitierte Arbeiten. 


die g 
Nahe 
Erge' 
e 
wobe 
durch 
:; 
1 
nm 
wird 
gescl 
> 
| 
für q 
zu (i 
—(k 
für j, 
| 
’ Ande 


lich 


R. Lenk: Rotierende Modellkerne 93 


Bib Ke = K+ + 56 Gp (6) 
die gestörte Streumatrix nach Potenzen von w entwickelt, so liefert die nullte 
Näherung Kin (5) eingesetzt dort einen Ausdruck, der aus dem entsprechenden 
Ergebnis bei der Hartree- Fock -Methode hervorgeht, wenn man dort v durch 


K® ersetzt. Der Vergleich mit (4.4) und 4.5) zeigt also 


| 
K® \ 
+5 Ts nm \, 7 
/ (7) 


nm 
wobei sich also 9 von 6 aus NIMRE. nur dadurch unterscheidet, daß dort v 
durch K® ersetzt wird. Der öK®/&o-Term in (7) liefert unter Verwendung von 
(3.12) 
1 \= 13, /nm| m zs 


2 = \ = am \ (8) 
nkm Een \ ow | / 


einen in wlinearen und einen in quadratischen Term. Im! nachsten Kapitel 
wird gezeigt, daB der lineare Term verschwindet, (7) kann also in der Form (4) 
geschrieben werden. Aus (7) und (8) oo 


<n | kyo nm \+ ¢.¢. 
nkm Een \ 77777 
PAG 


— 


(9) 


fir das 
Der Ausdruck (9) läßt sich wieder erheblich vereinfachen, wenn man die 
€ 


tm 


fiir j, verwendet. (9) erhält dann die Form ice 3.4 
nkm n 
(11) 


(11) ist im Vergleich mit (3.5) kompliziert durch die Zusatzterme, die von der 
Anderung der effektiven Wechselwirkungen K° durch die Kernrotation her- 
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rühren. Das Trägheitsmoment des Kerns hängt also zumindest prinzipiell 
von diesem Effekt ab. Es ergeben sich Abweichungen von den Aussagen eines 
einfachen Einteilchenmodells dadurch, daß dem Einteilchenmodell nicht die 
wirklichen Wechselwirkungen v, sondern effektive Wechselwirkungen K® 
(Rechengrößen!) zugrunde liegen, die selbst vom Zustand des Kerns abhängen. 
Ein bekanntes Beispiel für Abweichungen dieser Art vom einfachen Einteil- 
chenmodell ist die sogenannte Rearrangement-Energie®). In diesem Fall sind 
die Korrekturen am Einteilchenmodell groß; die Frage, ob entsprechende 
Korrekturen bei 6 erheblich sind und die Übereinstimmung mit experimen- 
tellen Werten für 9 verbessern, bleibt zunächst offen. Um sie zu beantworten, 
müßten die K®-Terme in (11) abgeschätzt werden, das wird aber in dieser 
Arbeit nicht durchgeführt. 
” 8. Die gestörte Streumatrix 

_ Zur Auswertung der Formel (7.11) für das Trägheitsmoment muß die 
gestérte Streumatrix K”, bzw. nach (7.6) die Abweichung von K® gegenüber 
K® berechnet werden. Die Bestimmungsgleichungen für K® und K® lauten 


0 0 


Pr Pr Pr Pr 
und 

K* =v =. (2) 


Bei K® tritt also ein geänderter Projektionsoperator Py und geänderte Energie- 
nenner €” auf. Beide sind völlig analog zu P% und e° definiert, beziehen sich aber 
zum Unterschied von diesen auf die gestörten Einteilchenfunktionen g, und 
die zugehörigen Einteilchenenergien e). Durch Subtraktion erhält man aus 
(1) und (2) os. 
n0 


Po pe P Pp 
wobei R durch 
(4) 
€ € 


definiert ist. Durch Zusammenfassung entsteht aus (3) unter Verwendung von 
(1) und (2) 

K° = K°+ K°R K®°= K°+ K°RK°+ KIORK’RK’+... (6) 
Die gestörte Streumatrix läßt sich also durch die ungestörte ausdrücken. Die 
in (5) angegebenen Entwicklungsglieder genügen zur Auswertung der in @ 
linearen und quadratischen Glieder, also zur Berechnung der in (7.11) gefor- 
derten Ableitungen 0K°/@@ und ®K®/öwodw. 


6) K. A. Brueckner, Physic. Rev. 110, 597 (1958); P. Mittelstaedt, Nuclear Phys. 
9, 116 (1958/59). 
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Hier soll nur noch gezeigt werden, daß die Diagonalelemente von 0K°/éw 
verschwinden. Dazu wird die Differenz 


<n m | K® — K®| n my = (nm |K® R K®| n m) 
2 
bis zur Ordnung w ausgewertet. Zunächst gilt 
<nm | K® — K®|nm> = <nm | K®| <p q nm) 
wobei die Summation über die Zwischenzustände einmal über die gestörten 
und einmal über die ungestörten Zustände geführt wird. Die Projektionsopera- 
toren Pp und Pr sorgen dann dafür, daß nur über Zustände onthe der 
gestörten bzw. ungestörten Fermikugel summiert wird, das ist durch ’ an- 
gedeutet. Entwickelt man |?) und |) in (7) entsprechend (3.12) nach den 
ungestörten Funktionen |k) und |l), so erhält man 
| 
(8) 
Panm 


wobei drei ee a Terme in (8) nur durch ... angedeutet sind. 


Die (K°— @ j’ — K®)-Terme von (8) verschwinden aus folgendem Grund: 
Das Wechselwirkungspotential v ist drehinvariant, es besitzt daher nur Matrix- 
elemente zwischen. Zuständen <p q| und |n m) zum gleichen Gesamtwert des 
Drehimpulses in Richtung der Figurenachse e des Kerns: 


= entspr. für 94, vp, if 

+ 


Das gleiche gilt für v-Ketten beliebiger ‚„‚Länge‘‘ und damit auch für die Ele- 5 


mente von K®, Dagegen ändert sich bei dem @j’-Übergang in (8) der Dreh- 
impuls, so daß diese Terme, wie oben behauptet, keinen Beitrag zu (8) liefern. __ 
Der Energienenner e® lautet ausführlich 


ep + es — en — em 
= + <p (i + ig)| P> + — <n | | n> — <m | | my. 


Der Erw vartungswert von @ j reduziert sich out dem jeweiligen Eigenwert - 


Zt keinen Telteng zu (10). Ubrig bleibt 


+ o[<p [io | p> + <q lio | ©, lio | n> lio | m) ]. (11) 


Nun reduzieren sich (wieder wegen der unvereinbaren Auswahlregeln bei den 
Elementen von j’ und K®) die durch (7.10) definierten Diagonalelemente von j, 
nach 


/ 
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fahren schrittweise verbessert wird, so ergibt sich 


“ [77 In 
Dieses Ergebnis wurde in Kapitel 7 bereits ne tie TREE 
he 
Anhang 


steht 
= [<km |v| <I |r| k’y —<k |x| m |v| k’ my 
Die Elemente von v können nach 7 
(km |v| Im) = (km |v| Im) — (k m |v| ml) pe 
= (km |v| Im) — (m k |v| Im) 


tauscht werden. Dann gilt beispielsweise für den ersten Term in (1) 
= <km |v|Imy <I |r| k’> 
ie 


& = (km |v| Im) 


>— (mk |v| Im) |r} k’> 
= (km |v k’ m)— (mk|vr,|k’ m) = (km 
Ganz entsprechend beweist man die drei übrigen Gleichungen 


= (km |v| k’ m> = (km |v r,| k’ m) 
= <m |r| <kl|v| = <km | ty 0| k’ m). 
Insgesamt wird damit aus (1) die Gleichung (6.17). 


_ Anregungen zu dieser Arbeit. 


Od Bei der Redaktion eingegangen am 17. August 1961. 


nachlassigt, in dieser De @K°/é@ nach (8) berechnet und dann das Ver. 


(13) 


In (6.16) wird zunächst die Vertauschung in dem Element von [V®, r] aus. 
geschrieben und die Definition (2.5) des Potentials V° verwendet. Dann ent. 


(1) 


in normalen und Austauschterm aufgespalten werden. Dabei sind die beiden in 
(2) angegebenen Möglichkeiten völlig gleichwertig, in (1) wird immer diejenige 
benutzt, beider die an < |r| > unmittelbar anschließenden Zustände nicht ver- 


(3) 


(4) 


Herrn Prof. Dr. W. Macke danke ich sehr fiir wertvolle Diskussionen und 
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auf die gerade zu untersuchenden Diagonalelemente von @K°/éw. Denkt man 
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LEBENSDARSTELLUNGEN DEUTSCHER NATURFORSCHER 


Herausgegeben von der 
Deutschen Akademie der Naturforscher Leopoldina 


Lieferbare Titel 
Nummer 5: 
Wissenschaftliche Selbstbiographie 
Von MAX PLANCK 
3. Auflage 


34 Seiten mit einem Bildnis. 1955. 8°. DM 1,35 


Orion: Über siebzig Jahre Geschichte der Physik werden lebendig, und inmitten dieser 
Zeitspanne die Geburtsstunde der Quantentheorie. Hinter jedem der bescheidenen Sätze 
dieses in seiner kristallenen Klarheit meisterlichen Selbstzeugnisses spürt man den großen 
Forscher. Das Heft wird für alle Zeiten eines der klassischen Dokumente der Natur- 
forschung bleiben. 


Nummer 6: 


Johann Friedrich Naumann 
der Altmeister der deutschen Vogelkunde 
Sein Leben und seine Werke 
Von P. THOMSEN 
Nach des Verfassers Tode bearbeitet und ergänzt 
von Prof. Dr. E. STRESEMANN, Berlin. 

212 Seiten mit 12 Abbildungen im Text. 1957. 8°. Leinen DM 13,— 
Die Vogelwelt: Für die ältere Ornithologengeneration ist der Name Naumann noch heute 
ein Begriff. Wenn auch die wenigsten mit ihm durch die Lektüre seines originalen Werkes 
vertraut sein mögen, so kennen sie doch die kurz nach der Jahrhundertwende erschienene 


Neubearbeitung, den ‚Neuen Naumann“. Die Biographie ist ein würdiges Mal für einen 
großen deutschen Ornithologen, dessen Wirken weit über seine Zeit hinausreichte. 


2 Oscar Hertwig 

Leben und Werk eines deutschen Biologen 

Von Prof. Dr. R. WEISSENBERG, Philadelphia 

63 Seiten mit 8 Bildnissen u. 1 Abbildung im Text. 1959. 8°. DM 3,60 

Durch die Veröffentlichung der Biographie von Oscar Hertwig soll die Erinnerung an 

sein Lebensbild und sein Lebenswerk wachgerufen werden, verdankt ihm doch eine ganze 

Generation von Biologen und Ärzten, die entweder unmittelbar seine Hörer und Schüler 


waren oder aus seinen Büchern und Schriften lernten, nicht nureine Bereicherung ihres 
Wissens, sondern auch die Formung ihres biologischen Weltbildes. 


Nummer 8: 


Friedrich Traugott Kützing 
1807-1893 
ss Aufzeichnungen und Erinnerungen 


Herausgegeben 
von R. H. WALTHER MÜLLER und RUDOLPH ZAUNICK, Halle ‘Saale 


300 Seiten mit 2 Bildnissen, 4 Abbildungen und einer Karte. 1960. 8°. 
Leinen DM 16,60 


Die hier vorliegenden „Aufzeichnungen und Erinnerungen‘ des weiland Professors am 
Nordhäuser Realgymnasium, Friedrich Traugott Kützing, bringen Material in reicher 
Fülle und fesselnder Darstellung. Diese Selbstbiographie laßt uns den Lebenslauf eines 
Menschen miterleben, der, aus der kinderreichen Familie eines Müllers in Ritteburg an 
der Unstrut stammend, sich schon als blutjunger Apothekerlehrling und -gehilfe zum 
Naturforscher entwickelte, der schon als Student der Pharmazie und Naturwissenschaften 
in Halle sich die ersten Sporen als systematischer Botaniker verdiente, und der neben 
seinem Amt als Lehrer der Chemie und Naturgeschichte als Systematiker der Diatomeen 
und Algen grundiegende Werke schuf. 
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(Für Freunde des Steeahinmels 


KALENDER FÜR STERNFREUNDE 1962 


Herausgegeben von Dr. PAUL AHNERT, Sternwarte Sonneberg in Thüringen 
Unter Mitwirkung von Prof. Dr. R. Müller, Sonnenobservatorium Wendelstein, 
und K. Schöbel, Hirschfelde 
188 Seiten mit 37 Abbildungen im Text. 1962. 8°. DM 4,— 


Der Kalender enthält die für die Beobachtungen erforderlichen Unterlagen für das Jahr 
1962, mit Ausnahme aller immer wiederkehrender Angaben, die im nachstehenden Buch 
„Beobachtungsobjekte für Liebhaberastronomen* enthalten sind. 

Anschließend folgen: Neuere astronomische Arbeiten und Entdeckungen / Künstliche 
Erdsätelliten 1961 / Die Rückseite des Mondes / Über Finsternisse sowie Merkur- und 
Venusdurchgänge / Bedeckung von hellen Fixsternen durch Planeten / Neuere Forschungen 
über den Tunguskameteor / Die vermißten Objekte in Messiers Katalog / Beobachtung 
veränder!icher Sterne / Versuche am Lyot-Koronoskop / Nomogramm zur Ermittlung 


der Zenitdistanz. 


Beobachtungsobjekte 
fiir Liebhaberastronomen 
Mit Hilfstafeln und Sternkarten 
Von Dr. PAUL AHNERT 
Sternwarte Sonneberg in Thiiringen 
104 Seiten mit 19 Abbildungen im Text. 1961. 8°. Kartoniert DM 3,— 
Die vorliegende Sammlung von Beobachtungsobjekten ist für die Benutzung neben obigem 
„Kalender für Sternfreunde‘‘ gedacht. Nach einem einführenden Aufsatz über das Fernrohr 
des Liebhaberastronomen, in dem ein Überblick über Bauart und Leistung gegeben wird, 
folgen alle jährlich in gleicher Form erforderlichen Angaben. Einige weitere Tafeln. die für 
fortgeschrittene Beobachter und den Astronomieunterricht von Interesse sind (Präzession, 
Refraktion, Extinktion, Zonenzeiten) sind zusätzlich aufgenommen. Außerdem enthält das 


Buch die Sternkarten, die früher im Kalender erschienen. Eine Übersichtskarte für den Mond 
und ein Nomogramm für Gestirnshöhen bilden den Abschluß. 


a Himmelswunder im Feldstecher 


Von Dipl.-Optiker RUDOLF BRANDT 


Sternwarte Sonneberg in Thüringen 


Meran 6., verbesserte Auflage 
iD oa Seiten mit 83 Abbildungen im Text und 1 Tafel. 1961. 8 


> Kartoniert DM 4,70 
Astronomischer Nachrichtendienst der Archenhold-Sternwarte: 
Es liegt hier ein neu bearbeitetes kleines Werk vor, das jedem Sternfreund und jedem 
Lehrer warmstens empfohlen sei, weil es zweierlei meistert und bietet: eine ausführliche 
Behandlung der Optik des Feldstechers und Kleinfernrohrs sowie deren Benutzung bei 
astronomischen Beobachtungen. Das Buch ist dazu angetan, eine gute Anleitung für eine 
leicht zu errichtende ,,Feldstecher-Sternwarte zu bilden. D.W. 
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